ETNA Kent State University
|
|
- Oscar Simmons
- 5 years ago
- Views:
Transcription
1 Eectronic Transactions on Numerica Anaysis. Voume 7, 1998, pp Copyright 1998,. ISSN ETNA A THEORETICAL COMPARISON BETWEEN INNER PRODUCTS IN THE SHIFT-INVERT ARNOLDI METHOD AND THE SPECTRAL TRANSFORMATION LANCZOS METHOD KARL MEERBERGEN y Abstract. The spectra transformation Lanczos method and the shift-invert Arnodi method are probaby the most popuar methods for the soution of inear generaized eigenvaue probems originating from engineering appications, incuding structura and acoustic anayses and fuid dynamics. The orthogonaization of the Kryov vectors requires inner products. Often, one empoys the standard inner product, but in many engineering appications one uses the inner product using the mass matrix. In this paper, we make a theoretica comparison between these inner products in the framework of the shift-invert Arnodi method. The concusion is that when the square-root of the condition number of the mass matrix is sma, the convergence behavior does not strongy depend on the choice of inner product. The theory is iustrated by numerica exampes arising from structura and acoustic anayses. The theory is extended to the discretized Navier-Stokes equations. Key words. Lanczos method, Arnodi s method, generaized eigenvaue probem, shift-invert. AMS subect cassifications. 65F Introduction. This paper is concerned with the soution of generaized eigenvaue probems of the form (1.1) Ax = Bx; A; B 2 R nn ; x 6= 0; where A may be symmetric or non-symmetric, and B is symmetric positive (semi) definite, by the spectra transformation Lanczos method [6, 18] and the shift-invert Arnodi method [17]. Appications incude the moda anaysis of structures without damping, which eads to (1.2) Ku =! 2 Mu; where K and M are symmetric matrices and often positive definite [9]. Typicay, the number of wanted eigenmodes for representing the structura properties for ow and mid frequencies ranges from a few tens to a few thousands. The moda extraction of acoustic finite eement modes aso eads to a probem of the form (1.2). The required number of eigenmodes is often sma, since the modes are usuay empoyed for a ow frequency anaysis. For the (Navier) Stokes probem, we have K (1.3) C C T 0 u p M 0 = 0 0 u p where M is symmetric positive definite, C is of fu rank and K is symmetric (Stokes [14]) or nonsymmetric (Navier-Stokes). This eigenvaue probem arises in the determination of the stabiity of a steady state soution. Here ony the rightmost eigenvaue is wanted [14, 3]. This paper concentrates on the soution of (1.2), but (1.3) wi aso be touched on. In both appications, M is a discretization of the continuous identity operator, i.e., the continuous inner product hx; yi is repaced by the discrete x T My. As a resut, the condition number of M is usuay sma. We study this specific case. The theory is iustrated by numerica exampes arising from rea appications. ; Received November 11, Accepted for pubication August 4, Recommended by R. Lehoucq. y LMS Internationa, Intereuvenaan 70, 3001 Leuven, Begium. Current address: Rutherford Appeton Laboratory, Chiton, Didcot, OX11 0QX, UK. (K.Meerbergen@r.ac.uk) 90
2 Kar Meerbergen 91 One approach to the soution of generaized eigenvaue probems is the shift-invert Arnodi method [17, 20, 15]. Instead of soving (1.1) directy, one soves the shifted and inverted probem (1.4) (A, B),1 Bx = x by the Arnodi method. The scaar is caed the shift, which expains the name shiftinvert. If (; x) is an eigenpair of (A, B),1 B,then( +,1 ;x) is an eigenpair of Ax = Bx. This reation demonstrates that s can be computed from s. Without oss of generaity, we assume a shift = 0 is used. In genera, A,1 B is a nonsymmetric matrix, even when A and B are symmetric, and this is the reason why the Arnodi method is used. However, when A is symmetric and B is symmetric positive definite, A,1 B is sef-adoint with respect to the B-inner product. This impies that the Lanczos method can be used, when the B-inner product x T By is empoyed instead of the standard inner product x T y. This idea was proposed by Ericsson [5] and Nour-Omid, Parett, Ericsson and Jensen [18]. A bock version was proposed by Grimes, Lewis and Simon [10]. In the case where B is positive semi-definite, which, e.g. arises in appications of the form (1.3), the B-semi-inner product can be used in the Lanczos method or the Arnodi method. This is suggested by Ericsson [5], Nour-Omid, Parett, Ericsson and Jensen [18], and Meerbergen and Spence [16] and appied to inearized and discretized Navier-Stokes equations by Lehoucq and Scott [12]. In this paper, we show by both anaysis and numerica exampes that if the square root of the condition number of the mass matrix B is sma, the choice of inner product does not infuence the convergence speed. The choice of inner product shoud be based on other criteria than rate of convergence. We iustrate this for two casses of appications. When A is symmetric, the use of the B-inner product reduces the Arnodi method to the Lanczos method. The Lanczos method has two advantages over the Arnodi method. First, the eigenvaues have quadratic error bounds and their convergence is we understood [19, 20]. Second, the cost per iteration consists of the action of A,1 B on a vector and the orthogonaization of the new iteration vector against the previous ones. The cost for the construction of the Kryov basis is smaer than for the Arnodi method, since ony the ast two basis vectors are used in the orthogonaization process. The Arnodi method uses a vectors. The Lanczos method uses the B-inner product which can be quite expensive compared to the standard inner product. The overa orthogonaization cost, however, can be much smaer than for the Arnodi method with standard inner product, when the number of iteration vectors is arge. This is often the case for a structura anaysis for ow and mid frequencies since a arge number of eigenmodes is wanted. The use of the standard inner product instead of the B-inner product may be preferred when A is nonsymmetric and the Arnodi method needs to be used anyway, so fu orthogonaization against a previous basis vectors cannot be avoided. Lehoucq and Scott [12] demonstrate for discretized Navier-Stokes appications that the B-inner product is more expensive than the standard inner product, but eads to a more reiabe Arnodi method. A side effect of the use of B-orthogonaization is that the approximate eigenvectors are B-orthogona, when A is symmetric. This is very natura since the exact eigenvectors corresponding to different eigenvaues are B-orthogona. In finite eement appications, it is assumed that the computed eigenmodes satisfy this property. This is automaticay satisfied by the Lanczos method with B-orthogonaization, but not by the Arnodi method. The pan of this paper is as foows. In x2, a theoretica comparison between standard and B orthogonaization is estabished. In x3, we present an easy way of obtaining B-orthogona eigenvectors from the Arnodi method when A is symmetric. In x4, we iustrate the theory by numerica exampes. Section 5 generaizes the ideas from x2 to the Navier-Stokes probem. Finay, we summarize the main concusions in x6. We assume computations in exact arithmetic.
3 92 Theoretica comparison between inner products 2. A reation between standard and B-orthogonaization. In this section, a theoretica study of the Arnodi method with standard orthogonaization and B-orthogonaization is estabished for the eigenvaue probem (1.2). The goa is to reate the residua norms as we as the Hessenberg matrix (which is tridiagona for the Lanczos method) and the computed eigenvaues for both types of inner product. The anaysis assumes exact arithmetic. First, in x2.1, some preiminaries and notation are presented. Second, x2.2 puts both types of orthogonaization into a singe theoretica framework: we present the agorithm and some properties. The reation between standard and B-inner products in the Arnodi method wi be formuated and derived in x Notation and preiminaries. This section is devoted to some notation and matrix properties. In genera, we use the Eucidean norm for vectors and matrices, denoted by kk 2 or kk. The matrix Frobenius norm is denoted by kk F. Let (C) denote the condition number of the matrix C. First, since B is a positive definite matrix, there exists L 2 R nn such that B = L T L. LEMMA 2.1. Consider V;W 2 R nk.letv T V = I, W T BW = I and et the coumns of V span the same space as the coumns of W. Then there is an S such that V = WS. Moreover, (S) =(W ) p (B). Proof. It is cear that there is an S such that V = WS. Hence V T BV = S T S and W T W = S,T S,1. Since kv k = 1, wehaveksk 2 = kv T BV k kbk, and since W T BW =(LW ) T (LW )=I, wehaveks,1 k 2 = kw T W k = k(lw ) T B,1 (LW )k kb,1 k. This competes the proof. We wi compare two agorithms that differ primariy in their choice of inner product or norm. We wi use the notation hx; yi to stand for a generic inner product, such as x T y or x T By. The notation is aso generaized to matrices V =[v 1 ;:::;v k ] and W =[w 1 ;:::;w ], as foows : hv;xi =[hv ;xi] k =1 2 Rk hw;vi =[hw i ;v i] (;k) (i;)=(1;1) 2 Rk : Since h; i is an inner product, it foows that hws;vzi = S T hw;viz. The B norm of a vector x is defined by kxk B = p xt Bx. TheB norm of a matrix C is defined by kck B = klck 2,whereB = L T L. Obviousy, for two matrices, V 2 R nk and W 2 R n,wehavekv T BWk kv k B kw k B. For a matrix C, the Kryov space K k (C; v 1 ) of order k with starting vector v 1 is defined by K k (C; v 1 ) = spanfv 1 ;Cv 1 ;C 2 v 1 ;:::;C k,1 v 1 g : We assume that a Kryov spaces of order k have dimension k. In practice, a space is represented by a basis. The foowing emma gives a reation between two different bases. LEMMA 2.2. Let V k ;W k 2 R nk be such that the first coumns of V k and the first coumns of W k form two bases for K (C; v 1 ) for =1; ::; k. Then there is a fu rank upper trianguar matrix S k 2 R kk such that V k = W k S k. The foowing emma shows the uniqueness of a normaized Kryov basis. LEMMA 2.3 (Impicit Q Theorem). ([8, Theorem 7.4.2]) Let the first coumns of V k and W k 2 R nk form two bases for K (C; v 1 ) for =1;:::;k and hv k ;V k i = I = hw k ;W k i. Then v i = w i i with i = 1 for i =1;:::;k.
4 Kar Meerbergen A genera (theoretica) framework for Kryov methods. The Arnodi and Lanczos methods for the soution of (1.2) are Kryov subspace methods, i.e., the eigenvaues and eigenvectors are computed from the proection of A,1 B on a Kryov space. The foowing agorithm covers both methods. Reca that hx; yi denotes the inner product, e.g. the standard inner product hx; yi = x T y or the B-inner product hx; yi = x T By. ALGORITHM 1. Genera framework for the Lanczos and Arnodi methods. 0. Given v 1 with hv 1 ;v 1 i =1. 1. For =1to k do 1.1. Form p = A,1 Bv Compute the Gram Schmidt coefficients h i = hv i ;p i;i =1;:::; Update q = p, P v i=1 ih i Compute norm h +1; = hq ;q i 1= Normaize : v +1 = q =h +1;. 2. Let H k =[h i ] (k+1;k) (i;)=(1;1) 2 Rk+1k where h i =0whenever i>+1. Let H k be the first k rows of H k. Let V k =[v 1 ;:::;v k ]. 3. Compute eigenpairs (; z) of H k, with z 2 R k, by the QR method. 4. Compute the Ritz vector x = V k z 2 R n. 5. Compute the residua norm = h k+1;k e T k z. Step 1.1 is performed by a matrix vector mutipication with B and the soution of a inear system with A. The soution of the inear system is usuay performed by a direct method, since one can take advantage of the fact that A needs to be factored ony once. Moreover, direct methods usuay give a soution with a sma backward error. This is required for the Arnodi method to find the eigenpairs of (1.1) [15]. Steps 1.1 to 1.5 compute a basis V k+1 =[v 1 ;:::;v k+1 ] of the Kryov space K k+1 (A,1 B;v 1 ) = spanfv 1 ;A,1 Bv 1 ;:::;(A,1 B) k v 1 g ; normaized such that hv k+1 ;V k+1 i = I. The normaization is performed by Gram-Schmidt orthogonaization. For reasons of numerica stabiity, practica impementations use reorthogonaization [4] in the Arnodi method and modified Gram-Schmidt with partia reorthogonaization in the Lanczos method [10]. The Gram-Schmidt coefficients are coected in the upper Hessenberg matrix H k. Note that H k is a k +1by k matrix. We denote the k k upper submatrix of H k by H k. By the eimination of p and q from Steps 1.1 to 1.5, it foows that (2.1a) (2.1b) A,1 BV k = V k+1 H k = V k H k + v k+1 h k+1;k e T : k This is the we known recurrence reation for the Arnodi and Lanczos methods. Usuay, k n, sothath k has much smaer dimensions than A and B. In Steps 3-4, an eigenpair (; x) is computed by the Gaerkin proection of A,1 B on K k, i.e., x 2K k and the residua r = A,1 Bx, x is orthogona to the Kryov space : (2.2) x = V k z with H k z = z : The s are the eigenvaues of H k and are caed Ritz vaues and the x s are the corresponding Ritz vectors. They form an approximate eigenpair of A,1 B. Reca that H k is an upper Hessenberg matrix, so the eigenpairs (; z) are efficienty computed by the QR method [8, x7.5]. The residua r = A,1 Bx, x aso foows from (2.1b) : (2.3) r = A,1 Bx, x = v k+1 h k+1;k e T k z
5 94 Theoretica comparison between inner products and the induced norm is very cheapy computed as (2.4) = hr;ri 1=2 = h k+1;k e T k z : This expains Step 5 in the agorithm. The residua norm is a measure of the accuracy of the eigenvaues. (See eigenvaue perturbation theory, e.g., [20, Ch. III] and [1, Ch. 2].) In the foowing, we aso use a bock formuation of the recurrence reation for the Ritz vectors. Let X = [x 1 ;:::;x ] = V k Z and D = diag( 1 ;:::; ) represent k Ritz pairs and R =[r 1 ;:::;r ] the corresponding residua terms of the form (2.3). Then it foows that (2.5) A,1 BX = X D + R ; R = h k+1;k v k+1 e T k Z : In the rest of this paper, we use the foowing notation to distinguish between the standard and B-orthogonaization. For standard orthogonaization, the Kryov vectors are denoted by V k+1 and the Hessenberg matrix is H k. They satisfy A,1 BV k = V k+1 H k ; H k = V T k+1a,1 BV k ; V T k+1v k+1 = I: For B-orthogonaization, the Kryov vectors are denoted by W k+1 and the Hessenberg matrix by T k. They satisfy A,1 BW k = W k+1 T k ; T k = W T k+1ba,1 BW k ; W T k+1bw k+1 = I: It is cear that if A is symmetric, the Hessenberg matrix T k = V T k BA,1 BV k is symmetric, hence tridiagona Comparison between both orthogonaization schemes. In this section, a reationship between standard and B-orthogonaization is estabished. Ceary, the computed Kryov spaces are the same, but the proections differ and this may ead to different eigenvaue approximations. The theoretica resut in this section uses Lemma 2.2, which makes the ink between two Kryov bases by use of an upper trianguar matrix. First, in Theorem 2.4, we reate h k+1;k and t k+1;k and H k and T k, and in Theorem 2.5, we reate the eigenvaues of T k and H k. THEOREM 2.4. Let v 1 = w 1 =kw 1 k. Then there is a matrix S k 2 R kk, such that (2.6) T k = S k H k S,1 k + E kek 2 h k+1;k (W k+1 ) (W k+1 ),1 t k+1;k h k+1;k (W k+1 ) ; and (S k );(W k+1 ) p (B) : Proof. Reca that V k+1 and H k satisfy the Arnodi recurrence reation (2.1a). Let S k+1 2 R k+1k+1 be the upper trianguar Choesky factor of V T k+1 BV k+1, i.e., V T k+1bv k+1 = S T k+1s k+1 : As a consequence, V k+1 S,1 k+1 forms a B-orthogona Arnodi basis for the Kryov space. Since Kryov bases are unique (see Lemma 2.3), W k+1 = V k+1 S,1 k+1. (Eventuay, the coumns
6 Kar Meerbergen 95 S k+1 must be mutipied by,1, see i in Lemma 2.3.) The Arnodi reation (2.1a) can be rewritten as A,1 B(V k S,1 k )=(V k+1s,1 k+1 )(S k+1h k S,1 k ); where S k is the k k principe submatrix of S k+1. Consequenty, T k = S k+1 H k S,1. k Decompose S k+1 H k S,1 Sk s H k = S,1 k 0 s k+1;k+1 h k+1;k e T k k T T k k = : t k+1;k e T k Then, we have where (2.7) and T k = S k H k S,1 k + E E = h k+1;k s,1 k;k set k t k+1;k = h k+1;k s k+1;k+1 s,1 k;k : The proof now foows from Lemma 2.1. This theorem says that when the square root of the condition number of B is sma, the residua terms are of the same order. It aso says that T k is a rank-one update of H k of the order of h k+1;k. The foowing theorem estabishes a reationship between eigenpairs of T k and H k. THEOREM 2.5. Let (; V k z) be a Ritz pair of the Arnodi method and et be the corresponding residua norm. Assume that v 1 = w 1 =kw 1 k. Then there exists an eigenvaue of T k such that, (Y k )(W k+1 ); where Y k is the matrix whose coumns contain the eigenvectors of T k. Proof. From (2.6) and the fact that H k z = z,wehave S k H k z, T k S k z =,ES k z T k S k z, S k z = ES k z: With y = S k z=ks k zk and = h k+1;k e T k z, it foows from (2.7) and et k S k = s k;k e T k that (2.8) kt k y, yk h k+1;ks,1 k;k ksket k S kz ks k zk h k+1;kkske T z ksk k = ks k zk ks k zk h k+1;ke T k z (S k+1 ) = (W k+1 ): Foowing the Bauer-Fike Theorem [20, Theorem 3.6], it foows that T k has an eigenvaue for which, (Y k )kt k y, yk ;
7 96 Theoretica comparison between inner products from which the proof foows. The eigenvaues of H k can simiary be reated to the eigenvaues of T k. When A is symmetric, so is T k. Therefore (Y k )=1, and the error bound becomes sharper. The most important concusion is that the Ritz vaues with a sma residua norm are amost the same for both the standard and the B-inner products. 3. Computation of B-orthogona eigenvectors in the Arnodi method. When the B- inner product is used and A and B are symmetric, the tridiagona matrix T k is symmetric, so T k has an orthogona set of eigenvectors z ; =1;:::;k. Therefore, the Ritz vectors W k z, =1;:::;kform a B-orthogona set of k vectors. When the standard inner product is used, the Ritz vectors are not necessariy B- orthogona. The foowing theorem shows the dependence of the B-orthogonaity on the separation between the eigenvaues and on the norms of the residuas. THEOREM 3.1. Let ( ;x ), =1;:::; < k be Ritz pairs and r the corresponding residuas computed by the Arnodi method with the standard inner product. Then, for 1 i;, we have x T i Bx kx i k B kx k B kr k B =kx k B + kr i k B =kx i k B i, provided that i 6=. Proof. Denote by D the diagona matrix containing 1 ;:::;, and et the coumns of X be the corresponding Ritz vectors. By mutipying (2.5) on the eft by X T B,wehave ; (3.1) X T BA,1 BX = X T BX D + X T BR : Since the eft-hand side is symmetric, (3.2) X T BX D + X T BR = D X T BX + R T BX and, so X T BX D, D X T BX = R T BX, X T BR : The (i; ) eement of this matrix eads to the resut of the theorem. This theorem shows that eigenvectors corresponding to different distinct eigenvaues are amost B-orthogona, but the eigenvectors corresponding to custered eigenvaues can ose B-orthogonaity. So, an additiona B-orthogonaization of the eigenvectors is desirabe. In the foowing, we use the notations of (2.5), i.e., D denotes the diagona matrix containing k Ritz vaues 1 ;:::; and the coumns of X denote the corresponding Ritz vectors. The most robust way to obtain B-orthogona Ritz vectors, is to compute the B- orthogona proection of A,1 B onto the range of X and compute new eigenpairs from this proection. This is estabished as foows. Let S be the upper trianguar Choesky factor of X T BX = S T S. Then, with Y = X S,1,wehaveY T BY = I. A Ritz pair obtained by B-orthogonaization is of the form (; y = Y z) with (; z) satisfying (3.3) Y T BA,1 BY z = z : This additiona proection eads to quadratic error bounds for the eigenvaues. (See e.g. the Kato-Tempe theorem [20, Theorem 3.8].) The matrix on the eft-hand side of (3.3) can be computed expicity, but can as we easiy be computed without the action of A,1 by the use of (3.1) : (3.4) Y T BA,1 BY = S D S,1 + Y T BR S,1
8 Kar Meerbergen 97 Box car interior exhaust pipe FIG.4.1.Meshes for the three eigenvaue probems in x4 (by courtesy of Bosa) In practice, we coud compute (; z) from (3.5) (S D S,1 + Y T BR S,1 )z = z instead of (3.3). However, there is an aternative to performing an additiona proection. The coumns of Y can be taken as the B-orthogona Ritz vectors. From the foowing theorem, it foows that if the residua norms of ( ;x ) for = 1;:::; are sma and p (B) is sma, the residuas for ( ;Y e ) for = 1;:::; are aso sma. Of course, the Ritz vaues do not satisfy a quadratic error bound. THEOREM 3.2. Reca the definition of D and X from Eq. (2.5). Let S be the upper trianguar Choesky factor of X T BX and et Y = X S,1.Then ka,1 BY, Y D k B (p +1) p (B)kR k : Proof. From (3.4), we have Y T BA,1 BY = S D S,1 + E with E = Y T BR S,1.SinceS is upper trianguar, and D diagona, S D S,1 = D + U where U is stricty upper trianguar. Note that S D S,1 + E is symmetric and so is U + E : U + E = U T + E T : This impies U, U T = E, E T 2kUk 2 F 2kEk 2 F and so kuk 2 kek F p kek 2. From (2.5), we have from which the proof foows. A,1 BY, Y D = Y (S D S,1, D )+R S,1 = Y U + R S,1 ka,1 BY, Y D k B kuk 2 + p (B)kR k 2 (p +1) p (B)kR k 2 4. Numerica exampes. A exampes have been generated using SYSNOISE, a software too for vibro-acoustic simuation [22]. The matrices A and B arise from acoustic and structura finite eement modes and are symmetric. We present resuts for the foowing probems.
9 98 Theoretica comparison between inner products PROBLEM 1. The first appication concerns the structura moda anaysis of a square box with Young moduus 2: Pa, Poisson coefficient 0:3 and voume density 7800 kg=m 3. The thickness of the materia is 0:01m. One of the six faces of the box is fixed. For the finite eement anaysis, we used 150 she eements. This anaysis eads to an eigenvaue probem (1.2). For this probem, A = K and B = M. The dimension of the probem is n = 852 and p (B) 1: We used a shift =,10 in the shift-invert transformation of (1.4). PROBLEM 2. The acoustic moda anaysis of a Vovo 480 3D car interior discretized with 744 cubic eements and 68 pentaedra finite eements eads to a probem of the form (1.2) with K positive semi-definite. The fuid has voume density 1:225 kg=m 3 and sound has a speed of 340 m/s (air). For this probem, A = K and B = M. The dimension of the probem is n = 1176 and p (B) 21:2. We used a shift =,10 in the shift-invert transformation of (1.4). PROBLEM 3. This concerns the couping of a structura and an acoustic probem. We consider a square box of 1 m 3 with Young moduus 2: Pa, Poisson coefficient 0:3 and voume density 7800 kg=m 3. The thickness of the materia is 0:01 m. One of the six faces of the box is fixed. The box is fied with air (speed of sound 340 m/s and voume density 1:29 kg=m 3 ). The box is discretized by 294 she eements, which represents the structura mode. The acoustic behavior of the fuid is modeed by a boundary eement mesh of 294 quadranguar eements. The matrix A is the structura stiffness matrix, whie B is the sum of the structura mass matrix and the added mass matrix coming from the acoustic mode. The dimension of the probem is n = 1692 and p (B) 1: Ashift =0was used. PROBLEM 4. The acoustic moda anaysis of a Jaguar X100MS exhaust pipe discretized with cubic finite eements eads to a probem of the form (1.2) with K positive semidefinite. The fuid has voume density 1:225kg=m 3 and the sound has a speed of 340 m/s (air). The dimension of the probem is n = and p (B) 42:8. For this probem, A = K and B = M. Theshift is chosen such that the 11 dominant eigenvaues of (A, B),1 B correspond to the eigenfrequencies between 0 and 200Hz Iustration for Theorem 3.1. We compare the Ritz vaues of (A +10B),1 B for Probems 1 and 2 and their residua norms after k = 10 Arnodi/Lanczos steps, started with a random initia vector. Theorem 3.1 states that the Ritz vaues obtained by B- orthogonaization and standard orthogonaization ie within a distance of p (B) where is the residua norm (2.4). For Probem 1, we have h k+1;k = 3: ,7, t k+1;k = 4: ,6 and (W k+1 )=24:755 and for Probem 2, we have h k+1;k =4: ,7, t k+1;k =3: ,7 and (W k+1 )=2:7. These vaues are consistent with Theorem 2.4. The matching significant digits of the respective Ritz vaues and the residua norms are shown in Tabes 4.1 and 4.2, except for the mutipe eigenvaue at 0:1 ( 1 ;:::; 4 ) for Probem 1, which is not dispayed. Ritz vaues that do not share any digits are not shown. A Ritz vaues satisfy Theorem Iustration of the B-orthogonaization of the Ritz vectors. When Ritz vectors are computed by Arnodi s method with standard orthogonaization, they are not B- orthogona. We coud perform an expicit B-orthogona proection by soving (3.5). Instead, we use the coumns of Y as Ritz vectors, as suggested in x3. We iustrate Theorem 3.2, which states that the coumns of Y are sufficienty accurate Ritz vectors. The are unchanged. The resuts reported come from the Arnodi computations from x4.1. The Ritz vectors before and after B-orthogonaization are denoted by x and y respectivey. Tabes 4.3 and 4.4 show the expicity computed residua norms before and after B-orthogonaization. (Note that the (2) in Tabes 4.3 and 4.4 corresponds we to the in Tabes 4.1 and 4.2 respectivey. The difference is that is computed from (2.4) and (2) is the expicity computed residua norm.)
10 Kar Meerbergen 99 TABLE 4.1 Matching the eigenvaues of H k and T k after 10 steps for Probem 1. The vaue is defined by (2.4). hx; yi = x T y hx; yi = x T By 5 5:94110,6 710,9 410,9 6 3: ,6 510,8 410, ,6 910,7 610, ,7 210,7 210, ,7 310, ,7 210,8 TABLE 4.2 Matching the eigenvaues of H k and T k after 10 steps for Probem 2. The vaue is defined by (2.4). hx; yi = x T y hx; yi = x T By 1 9: ,2 0: 0: 2 6: ,6 210,13 210,13 3 2:20810,6 210,9 210,9 4 1:810,6 210,7 210,7 5 1:36410,6 310,8 310, ,7 210,7 310, ,7 310,7 210, ,7 110,7 110, ,8 810, ,8 210,8 The resuts are consistent with Theorem 3.2. In Tabe 4.3, (B) 4 10,6, which is much arger than (2). This is possibe since y 4 4 is a inear combination of x 1 ;:::;x 4 and (B) 4 depends on kr 4 k (2) 1 10, Iustration for the impicity restarted Arnodi/Lanczos methods. In practica cacuations, if sma residua norms need to be obtained for the desired eigenvaues, the Kryov subspaces often become very arge. In the iterature, some remedies against the growth of this subspace are proposed. One idea is to restart the Kryov method with a new poe in (1.4), as suggested in [10, 7]. In this paper, we use the impicity restarted Lanczos and Arnodi methods with exact shifts [21, 11]. Roughy speaking, the impicity restarted Arnodi method is mathematicay equivaent to the foowing scheme.
11 100 Theoretica comparison between inner products TABLE 4.3 Residua norms for Probem 1 before and after B-orthogonaization. We denote (2) = k(a+10b),1 Bx, x k and (B) = k(a +10B),1 By, y k B. (2) (B) 1 310,5 510, ,9 810, ,8 210, ,9 210, ,9 110, ,8 110, ,7 310, ,7 810, ,7 310, ,7 410,8 TABLE 4.4 Residua norms for Probem 2 before and after B-orthogonaization. We denote (2) = k(a+10b),1 Bx, x k and (B) = k(a +10B),1 By, y k B. (2) (B) 1 610,17 410, ,13 410, ,9 510, ,8 310, ,7 110, ,7 210, ,7 210, ,7 110, ,8 110, ,8 910,8 ALGORITHM 2. (impicity) restarted Arnodi method 0. Given is an initia vector v Iterate : 1.1. Form V k+1 and H k by k steps of Arnodi Compute Ritz pairs (; x) and residua norm Get a new initia vector v 1 = V k z. Unti TOL The impicity restarted Arnodi method is an efficient and reiabe impementation of this agorithm. For the probems that are soved in this paper, the new initia vector v 1 is a inear combination of the wanted Ritz vectors, so that unwanted eigenvaues do not show up in foowing iterations. (For a precise definition and practica agorithms, see [21, 11].) In genera, the impicity restarted Arnodi and Lanczos methods do not have the same subspaces after two iterations, since the Ritz vectors are different and so is the inear combination. This can ead to different subspaces and therefore to different rates of convergence. The numerica resuts are obtained using the ARPACK routines dsaupd (Lanczos) and dnaupd (Arnodi) [13]. Note that ARPACK uses the Arnodi process, i.e., fu (re)orthogonaization,
12 Kar Meerbergen 101 TABLE 4.5 Comparison between impicity restarted Lanczos (x T By) and Arnodi (x T y) for Probem 3 hx; yi inear soves matrix vector iterations time products with B (sec.) x T By : x T y : TABLE 4.6 Comparison between impicity restarted Lanczos (x T By) and Arnodi (x T y) for Probem 4 hx; yi inear soves matrix vector iterations time products with B (sec.) x T By x T y for symmetric probems, rather than the more efficient Lanczos process with partia reorthogonaization. This impies that timings reductions are possibe for the resuts of the Lanczos method with B-orthogonaization. The matrixb was representedby a compressed sparse row matrix storage format, which aows for efficient storage and matrix vector products. The inear system sovers where soved by a bock skyine (or profie) sover. A computations were carried out within the SYSNOISE [22] environment on an HP PA 7100 C110 workstation. Tabe 4.5 contains the numerica resuts for Probem 3. The 20 eigenmodes nearest 0 are required. The Impicity Restarted Arnodi/Lanczos methods were run with the same initia vector, Kryov subspace dimension k = 30, and reative residua toerance TOL = 10,5. The number of iterations are equa for both methods. The Lanczos method with B-orthogonaization is more expensive due to the additiona matrix vector products in the reorthogonaization of the B-orthogonaization. For the Arnodi method with standard orthogonaization, the number of matrix vector products with B is equa to the number of soves with A, B pus the number needed for the B-orthogonaization of the Ritz vectors (x3). Tabe 4.6 contains the numerica resuts for Probem 4. The eigenfrequencies between 0 and 250 Hz are sought. The impicity restarted Arnodi/Lanczos methods were run with the same initia vector, Kryov subspace dimension k = 21, and reative residua toerance TOL =10,5. A matrix vector product by B is cheap compared to a back transformation. The concusions are very simiar to Probem 3, but the differences between the CPU times are ess pronounced. 5. Extension to (Navier) Stokes probems. Reca the (Navier) Stokes probem (1.3). For this probem, B is typicay positive semi-definite. The use of the B-semi-inner product in the Lanczos method was suggested by Ericsson [5], Nour-Omid, Parett, Ericsson and Jensen [5]. The use of this semi-inner product was ustified by Ericsson [5] and Meerbergen and Spence [16]. Moreover, it purifies the Ritz vaues from possibe contamination arising from the singuarity of B. For the detais, see [5, 18, 16]. Ciffe, Garratt, Goding and Spence [2] suggested the use of the matrix P = I instead of M in the inner product. This P -inner product has the same purification property as the B-inner product (see [16, Theorem 2] with M = I). Thus, an aternative to the Arnodi method with
13 102 Theoretica comparison between inner products B-orthogonaization is the Arnodi method with P -orthogonaization. In fact, the situation is very simiar to acoustic and structura eigenvaue probems. Indeed, decompose Uk+1 U ~ V k+1 = and W k+1 k+1 = ~Q k+1 Q k+1 and normaize V k+1 such that V T PV k+1 k+1 = U T U k+1 k+1 = I and W k+1 such that W T BW k+1 k+1 = U ~ T M ~ k+1 U k+1 = I. Consequenty, Theorem 2.4 is sti vaid with the ony difference that (B) shoud be repaced by (M ), V k+1 by U k+1 and W k+1 by U ~ k Concusions and fina remarks. In this paper, we estabished a theoretica comparison between standard orthogonaization and B-orthogonaization for the soution of eigenvaue probems, for which the square root of the condition number of the mass matrix is sma. This is often the case in acoustic and structura eigenvaue probems and for the discretized (Navier) Stokes equations. Roughy speaking, the orthogonaization seems to pay a minor roe in the convergence of the eigenvectors in the Kryov space. The theoretica resuts are extended to the Navier-Stokes probem as we, where the Arnodi method with B-semiorthogonaization and the Arnodi method with P -semi-orthogonaization behave simiary. Foowing the theory and our numerica experiments the type of inner product shoud not be chosen on the basis of rates of convergence. The decision shoud be based on cost per iteration, or reiabiity. For exampe, it foows from the numerica resuts for Probem 4 that there is no advantage in performance for the Arnodi method with standard orthogonaization with respect to the Lanczos method with B-orthogonaization, though the inner product is more expensive for the Lanczos method. Moreover, when the number of iteration vectors is high, it is expected that the Gram-Schmidt process with reorthogonaization becomes prohibitive, so that the Lanczos method becomes much cheaper than the Arnodi method. We aso noticed that the Lanczos method with B-orthogonaization is more reiabe than the Arnodi method with standard orthogonaization: sometimes, sought-after eigenvaues are missed by the Arnodi method. For Navier-Stokes appications, where the Arnodi method shoud be used anyway since A is nonsymmetric, the use of the P -inner product rather than the B-inner product may ead to a more efficient code. Lehoucq and Scott [12] report that the B-inner product seems to be more reiabe than the standard inner product. They do not report resuts for the P -inner product. 7. Acknowedgment. The author is gratefu to John Lewis who made usefu comments on the impementation and appication of the spectra transformation Lanczos method and suggested some other materia for further discussion, and to Jean-Pierre Coyette and Christophe Lecomte for hepfu comments on the appications in the numerica exampes. The author aso thanks the referees for the many hepfu suggestions that improved the presentation of the paper. REFERENCES [1] F. CHATELIN, Eigenvaues of matrices, John Wiey and Sons, [2] K. CLIFFE, T.GARRATT, J.GOLDING, AND A. SPENCE, On the equivaence of bock matrices under Arnodi s method, Tech. Rep. Preprint 94/08, Schoo of Mathematica Sciences, University of Bath, Bath, UK, [3] K. CLIFFE, T. GARRATT, AND A. SPENCE, Eigenvaues of bock matrices arising from probems in fuid mechanics, SIAM J. Matrix Ana. App., 15 (1994), pp [4] J. DANIEL, W. GRAGG, L. KAUFMAN, AND G. STEWART, Reorthogonaization and stabe agorithms for updating the Gram-Schmidt QR factorization, Math. Comp., 30 (1976), pp [5] T. ERICSSON, A generaised eigenvaue probem and the Lanczos agorithm, in Large Scae Eigenvaue Probems, J. Cuum and R. Wioughby, eds., Esevier Science Pubishers BV, 1986, pp
14 Kar Meerbergen 103 [6] T. ERICSSON AND A. RUHE, The spectra transformation Lanczos method for the numerica soution of arge sparse generaized symmetric eigenvaue probems, Math. Comp., 35 (1980), pp [7] T. GARRATT, G. MOORE, AND A. SPENCE, Two methods for the numerica detection of Hopf bifurcations, in Bifurcation and Chaos: Anaysis, Agorithms, Appications, R. Seyde, F. Schneider, and H. Troger, eds., Birkhäuser, 1991, pp [8] G. GOLUB AND C. VAN LOAN, Matrix computations, The Johns Hopkins University Press, 3rd ed., [9] R. GRIMES, J. LEWIS, AND H. SIMON, Eigenvaue probems and agorithms in structura engineering, in Large Scae Eigenvaue Probems, J. Cuum and R. Wioughby, eds., Esevier Science Pubishers B.V., 1986, pp [10], A shifted bock Lanczos agorithm for soving sparse symmetric generaized eigenprobems, SIAMJ. Matrix Ana. App., 15 (1994), pp [11] R. B. LEHOUCQ, Restarting an Arnodi reduction, Preprint MCS-P , Argonne Nationa Laboratory, Argonne, IL 60439, [12] R. B. LEHOUCQ AND J. A. SCOTT, Impicity restarted Arnodi methods and eigenvaues of the discretized Navier-Stokes equations, Technica Report SAND J, Sandia Nationa Laboratory, Abuquerque, New Mexico, [13] R. B. LEHOUCQ, D.C.SORENSEN, AND C. YANG, ARPACK USERS GUIDE: Soution of Large Scae Eigenvaue Probems with Impicity Restarted Arnodi Methods, SIAM, Phiadephia, PA, [14] D. MALKUS, Eigenprobems associated with the discrete LBB condition for incompressibe finite eements, Internat. J. Engrg. Sci., 19 (1981), pp [15] K. MEERBERGEN AND D. ROOSE, Matrix transformations for computing rightmost eigenvaues of rea nonsymmetric matrices, IMA J. Numer. Ana., 16 (1996), pp [16] K. MEERBERGEN AND A. SPENCE, Impicity restarted Arnodi and purification for the shift-invert transformation, Math. Comp., 66 (1997). [17] R. NATARAJAN, An Arnodi-based iterative scheme for nonsymmetric matrix pencis arising in finite eement stabiity probems, J. Comput. Phys., 100 (1992), pp [18] B. NOUR-OMID, B. PARLETT, T. ERICSSON, AND P. JENSEN, How to impement the spectra transformation, Math. Comp., 48 (1987), pp [19] B. PARLETT, The symmetric eigenvaue probem, Prentice Ha series in computationa mathematics, Prentice- Ha, [20] Y. SAAD, Numerica methods for arge eigenvaue probems, Agorithms and Architectures for Advanced Scientific Computing, Manchester University Press, Manchester, UK, [21] D. C. SORENSEN, Impicit appication of poynomia fiters in a k-step Arnodi method, SIAM J. Matrix Ana. App., 13 (1992), pp [22] SYSNOISE REV 5.4, A software too for vibro-acoustic simuation. LMS Internationa, Intereuvenaan 70, 3001 Leuven, Begium, 1999.
An implicit Jacobi-like method for computing generalized hyperbolic SVD
Linear Agebra and its Appications 358 (2003) 293 307 wwweseviercom/ocate/aa An impicit Jacobi-ike method for computing generaized hyperboic SVD Adam W Bojanczyk Schoo of Eectrica and Computer Engineering
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationBALANCING REGULAR MATRIX PENCILS
BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity
More informationPartial permutation decoding for MacDonald codes
Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics
More informationTHE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES
THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationPreconditioned Locally Harmonic Residual Method for Computing Interior Eigenpairs of Certain Classes of Hermitian Matrices
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Preconditioned Locay Harmonic Residua Method for Computing Interior Eigenpairs of Certain Casses of Hermitian Matrices Vecharynski, E.; Knyazev,
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More information18-660: Numerical Methods for Engineering Design and Optimization
8-660: Numerica Methods for Engineering esign and Optimization in i epartment of ECE Carnegie Meon University Pittsburgh, PA 523 Side Overview Conjugate Gradient Method (Part 4) Pre-conditioning Noninear
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationDistributed average consensus: Beyond the realm of linearity
Distributed average consensus: Beyond the ream of inearity Usman A. Khan, Soummya Kar, and José M. F. Moura Department of Eectrica and Computer Engineering Carnegie Meon University 5 Forbes Ave, Pittsburgh,
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vo. 2, No. 5, pp. 909 926 c 2000 Society for Industria and Appied Mathematics A DEFLATED VERSION OF THE CONJUGATE GRADIENT ALGORITHM Y. SAAD, M. YEUNG, J. ERHEL, AND F. GUYOMARC H
More informationAbsolute Value Preconditioning for Symmetric Indefinite Linear Systems
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.mer.com Absoute Vaue Preconditioning for Symmetric Indefinite Linear Systems Vecharynski, E.; Knyazev, A.V. TR2013-016 March 2013 Abstract We introduce
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More informationCryptanalysis of PKP: A New Approach
Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in
More informationORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION
J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationA SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS
A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY
More informationAlgorithms to solve massively under-defined systems of multivariate quadratic equations
Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
More informationResearch Article Deflated BiCG with an Application to Model Reduction
Hindawi Pubishing Corporation Mathematica Probems in Engineering Voume 25, Artice ID 3729, 8 pages http://dx.doi.org/.55/25/3729 Research Artice Defated BiCG with an Appication to Mode Reduction Jing Meng,
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationT.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA
ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network
More informationThe EM Algorithm applied to determining new limit points of Mahler measures
Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationHomotopy Perturbation Method for Solving Partial Differential Equations of Fractional Order
Int Journa of Math Anaysis, Vo 6, 2012, no 49, 2431-2448 Homotopy Perturbation Method for Soving Partia Differentia Equations of Fractiona Order A A Hemeda Department of Mathematics, Facuty of Science
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationSASIMI: Sparsity-Aware Simulation of Interconnect-Dominated Circuits with Non-Linear Devices
SASIMI: Sparsity-Aware Simuation of Interconnect-Dominated Circuits with Non-Linear Devices Jitesh Jain, Stephen Cauey, Cheng-Kok Koh, and Venkataramanan Baakrishnan Schoo of Eectrica and Computer Engineering
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More informationResearch Article Numerical Range of Two Operators in Semi-Inner Product Spaces
Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda
More informationEquilibrium of Heterogeneous Congestion Control Protocols
Equiibrium of Heterogeneous Congestion Contro Protocos Ao Tang Jiantao Wang Steven H. Low EAS Division, Caifornia Institute of Technoogy Mung Chiang EE Department, Princeton University Abstract When heterogeneous
More informationUsing constraint preconditioners with regularized saddle-point problems
Comput Optim App (27) 36: 249 27 DOI.7/s589-6-94-x Using constraint preconditioners with reguarized sadde-point probems H.S. Doar N.I.M. Goud W.H.A. Schiders A.J. Wathen Pubished onine: 22 February 27
More informationPrimal and dual active-set methods for convex quadratic programming
Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:
More informationAn Algorithm for Pruning Redundant Modules in Min-Max Modular Network
An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationApproximated MLC shape matrix decomposition with interleaf collision constraint
Approximated MLC shape matrix decomposition with intereaf coision constraint Thomas Kainowski Antje Kiese Abstract Shape matrix decomposition is a subprobem in radiation therapy panning. A given fuence
More informationResearch of Data Fusion Method of Multi-Sensor Based on Correlation Coefficient of Confidence Distance
Send Orders for Reprints to reprints@benthamscience.ae 340 The Open Cybernetics & Systemics Journa, 015, 9, 340-344 Open Access Research of Data Fusion Method of Muti-Sensor Based on Correation Coefficient
More informationIntroduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationOn Complex Preconditioning for the solution of the Kohn-Sham Equation for Molecular Transport
On Compex Preconditioning for the soution of the Kohn-Sham Equation for Moecuar Transport Danie Osei-Kuffuor Lingzhu Kong Yousef Saad James R. Cheikowsky December 23, 2007 Abstract This paper anayzes the
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More informationNumerical linear algebra and some problems in computational statistics
IASC28: December 5-8, Yokohama, Japan Numerica inear agebra and some probems in computationa statistics Zdeněk Strakoš 1 1 Institute of Computer Science, Academy of Sciences, Pod Vod. věží 2, 182 7 Prague
More informationInvestigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l
Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationInternational Journal of Mass Spectrometry
Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This artice appeared in a journa pubished by Esevier. The attached copy is furnished to the author for interna non-commercia research and education use, incuding for instruction at the authors institution
More informationGlobal Optimality Principles for Polynomial Optimization Problems over Box or Bivalent Constraints by Separable Polynomial Approximations
Goba Optimaity Principes for Poynomia Optimization Probems over Box or Bivaent Constraints by Separabe Poynomia Approximations V. Jeyakumar, G. Li and S. Srisatkunarajah Revised Version II: December 23,
More information17 Lecture 17: Recombination and Dark Matter Production
PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was
More informationMode in Output Participation Factors for Linear Systems
2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.
More informationMath 220B - Summer 2003 Homework 1 Solutions
Math 0B - Summer 003 Homework Soutions Consider the eigenvaue probem { X = λx 0 < x < X satisfies symmetric BCs x = 0, Suppose f(x)f (x) x=b x=a 0 for a rea-vaued functions f(x) which satisfy the boundary
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationMaejo International Journal of Science and Technology
Fu Paper Maejo Internationa Journa of Science and Technoogy ISSN 1905-7873 Avaiabe onine at www.mijst.mju.ac.th A study on Lucas difference sequence spaces (, ) (, ) and Murat Karakas * and Ayse Metin
More informationSubstructuring Preconditioners for the Bidomain Extracellular Potential Problem
Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationExpectation-Maximization for Estimating Parameters for a Mixture of Poissons
Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating
More informationHaar Decomposition and Reconstruction Algorithms
Jim Lambers MAT 773 Fa Semester 018-19 Lecture 15 and 16 Notes These notes correspond to Sections 4.3 and 4.4 in the text. Haar Decomposition and Reconstruction Agorithms Decomposition Suppose we approximate
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationStatistical Learning Theory: a Primer
??,??, 1 6 (??) c?? Kuwer Academic Pubishers, Boston. Manufactured in The Netherands. Statistica Learning Theory: a Primer THEODOROS EVGENIOU AND MASSIMILIANO PONTIL Center for Bioogica and Computationa
More informationModal analysis of a multi-blade system undergoing rotational motion
Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
SIAM J. NUMER. ANAL. Vo. 0, No. 0, pp. 000 000 c 200X Society for Industria and Appied Mathematics VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW
More informationResearch on liquid sloshing performance in vane type tank under microgravity
IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the
More informationSequential Decoding of Polar Codes with Arbitrary Binary Kernel
Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient
More informationComponentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems
Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationTorsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg
Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Deft University of Technoogy Marijn Drienburg October 2017 Contents 1 Introduction 2 1.1 Hand Cacuation....................................
More informationSolution methods for linear discrete ill-posed problems for color image restoration
BIT manuscript No. (wi be inserted by the editor) Soution methods for inear discrete i-posed probems for coor image restoration A. H. Bentbib M. E Guide K. Jbiou E. Onunwor L. Reiche Received: date / Accepted:
More informationAvailable online at ScienceDirect. Procedia Computer Science 96 (2016 )
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia Computer Science 96 (206 92 99 20th Internationa Conference on Knowedge Based and Inteigent Information and Engineering Systems Connected categorica
More informationAvailable online at ScienceDirect. IFAC PapersOnLine 50-1 (2017)
Avaiabe onine at www.sciencedirect.com ScienceDirect IFAC PapersOnLine 50-1 (2017 3412 3417 Stabiization of discrete-time switched inear systems: Lyapunov-Metzer inequaities versus S-procedure characterizations
More informationCONJUGATE GRADIENT WITH SUBSPACE OPTIMIZATION
CONJUGATE GRADIENT WITH SUBSPACE OPTIMIZATION SAHAR KARIMI AND STEPHEN VAVASIS Abstract. In this paper we present a variant of the conjugate gradient (CG) agorithm in which we invoke a subspace minimization
More informationGeneral Decay of Solutions in a Viscoelastic Equation with Nonlinear Localized Damping
Journa of Mathematica Research with Appications Jan.,, Vo. 3, No., pp. 53 6 DOI:.377/j.issn:95-65...7 Http://jmre.dut.edu.cn Genera Decay of Soutions in a Viscoeastic Equation with Noninear Locaized Damping
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationThe Group Structure on a Smooth Tropical Cubic
The Group Structure on a Smooth Tropica Cubic Ethan Lake Apri 20, 2015 Abstract Just as in in cassica agebraic geometry, it is possibe to define a group aw on a smooth tropica cubic curve. In this note,
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationApproximated MLC shape matrix decomposition with interleaf collision constraint
Agorithmic Operations Research Vo.4 (29) 49 57 Approximated MLC shape matrix decomposition with intereaf coision constraint Antje Kiese and Thomas Kainowski Institut für Mathematik, Universität Rostock,
More informationBayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?
Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine
More informationJENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF SEVERAL VARIABLES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Voume 128, Number 7, Pages 2075 2084 S 0002-99390005371-5 Artice eectronicay pubished on February 16, 2000 JENSEN S OPERATOR INEQUALITY FOR FUNCTIONS OF
More informationNumerical simulation of a high viscosity bubble column
20th Internationa Congress on Modeing and Simuation, Adeaide, Austraia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Numerica simuation of a high viscosity bubbe coumn Danio Carvajaa, Victor Meendez-Vejara,
More informationDo Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix
VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationReichenbachian Common Cause Systems
Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,
More informationAALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009
AALBORG UNIVERSITY The distribution of communication cost for a mobie service scenario by Jesper Møer and Man Lung Yiu R-29-11 June 29 Department of Mathematica Sciences Aaborg University Fredrik Bajers
More information