# SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY

Save this PDF as:

Size: px
Start display at page:

Download "SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY"

## Transcription

1 SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY KLAUS NEYMEYR ABSTRACT. Multigrid techniques can successfully be applied to mesh eigenvalue problems for elliptic differential operators. They allow to compute a few of the smallest eigenvalues and the eigenfunctions within optimal or quasi-optimal computational complexity. A review is given on three multigrid solvers for elliptic eigenproblems. Numerical results are presented for a simple model problem. 1. INTRODUCTION Consider the eigenvalue problem for an elliptic partial differential operator. Numerical approximations of its eigenvalues and eigenfunctions can be gained by a finite element discretization of the variational problem associated with the differential operator. The discretization results in a generalized matrix eigenvalue problem of the form (1.1) Ax = λmx. The discretization matrix A and the mass matrix M are sparse with only a small number of nonzero elements per row. In typical applications, the dimension of these matrices may exceed and, at the same time, only a few of the smallest eigenvalues (in modulus) together with the eigenvectors are to be determined. These smallest eigenvalues characterize, for example, the base frequencies of a vibrating mechanical structure modeled by an eigenvalue problem for an elliptic partial differential operator. Most of the classical methods of numerical linear algebra for solving a generalized matrix eigenvalue problem cannot be applied to (1.1) because typically, the computer storage for a full representation of A and M or any factorization of these matrices is not available. The QR method, the Rayleigh quotient iteration or inverse iteration, for example, would all demand to generate matrix factorizations and to store the more dense matrix factors. In any case, none of these iterations would allow to determine an eigenpair within O(n) floating point operations for A and M being n n matrices. In contrast to this, boundary value problems for the same elliptic partial differential operators can be solved with optimal complexity or quasi-optimal computational complexity (i.e. O(n) or O(n log n)) by means of multigrid or domain decomposition methods. The aim of this paper is to rivet on eigensolvers for elliptic eigenvalue problems with multigrid efficiency, i.e. with optimal or quasi-optimal complexity. In Section 2 we first point out common elements of eigensolvers for mesh eigenvalue problems and then we briefly review and classify the following eigensolvers: (I) The Rayleigh quotient multigrid minimization, [4, 5, 16], (II) Direct multigrid eigensolvers, [3, 7, 8, 10], (III) Eigensolvers using multigrid as a linear solver, [14, 19]. Finally, Section 3 contains numerical results for the Laplacian in 2D. 1

2 2 KLAUS NEYMEYR 2. MULTIGRID EIGENSOLVERS In order to give a concise description of multigrid eigensolvers we first separate those parts which are common to any such iteration: (1) Subspace extension: Any eigensolver for computing only a single eigenvalue/eigenvector pair (typically the smallest eigenvalue), can simply be extended to a subspace iteration in order to determine a modest number of the smallest eigenvalues and eigenvectors. One can implement a subspace eigensolver by applying a given vector iteration to each of the Ritz vectors spanning an actual approximating subspace. Subsequent application of the Rayleigh-Ritz procedure serves to compute the new Ritz values and Ritz vectors. Alternatively, subspace computations can be realized by using the less stable deflation technique. (2) Nested iteration and adaptivity: The concepts of nested iteration and adaptivity, which are well known from boundary value problem solvers, can be transferred to the eigenvalue problem. Then nested iteration starts with the solution of a coarse grid eigenproblem. The following steps are successive prolongations of the actual eigenvector approximations to refined grids in combination with sufficiently many iterations of the eigensolver on each grid level. Nested iteration can be combined with the concept of adaptivity. By means of adaptively generated grids, numerical approximation of eigenvalues and eigenvectors within a prescribed tolerance can often be gained with only a small portion of the work necessary if uniform grid refinement is employed. An adaptive algorithm requires appropriate error estimators for the iteration error and for the discretization error. The first estimator is used to define a stopping criterion for the iterative eigensolver on the actual grid whereas the second estimator controls the mesh refinement. Realizations of such estimators are, for example, suggested in [5, 11, 20]. Let us turn to the multigrid eigensolvers I III, whose discussion because of (1) can be restricted to a single-vector iteration. I. Rayleigh quotient multigrid minimization (RQMG): We consider the eigenvalue problem for a self-adjoint elliptic partial differential operator. Then the discrete eigenproblem (1.1) can be treated as an optimization problem for the Rayleigh quotient (2.1) λ(x) = (x, Ax) (x, Mx). By the Courant-Fischer principle the minimum of (2.1) equals the smallest eigenvalue of (A, M) and is taken in the corresponding eigenvector. Hence iterative minimization of (2.1) can serve as an eigensolver. This minimization can be realized by means of a multigrid procedure. One applies a coordinate relaxation scheme, i.e. for each coordinate direction d k i (which is associated with the i-th finite element function on a certain grid level k) one computes the minimum (2.2) λ(x + ϑ d k (x + ϑd k i i ) = min, A(x + ϑdk i )) ϑ R (x + ϑd k i, M(x + ϑdk i )), which is at the same time the smallest Ritz value of (A, M) in the 2D space span{x, d k i }. The new iterate is x+ϑ d k i. A multigrid cycle of RQMG consists in a successive minimization of λ(x) for all finite element functions on all grid levels [4, 16, 18]. Grid independent convergence estimates for RQMG have been presented in [18]. For a generalization of RQMG to non-selfadjoint elliptic operators see [5].

3 SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY 3 II. Direct multigrid eigensolvers: Direct multigrid eigensolvers and the third class of eigensolvers (see III) are related to approximate variants of classical iterative eigensolvers, namely inverse iteration and the Rayleigh quotient iteration [23]. Inverse iteration applied to (1.1) generates a sequence of iterates x (k) by solving the linear systems (2.3) (A σm)x (k+1) = Mx (k), k = 0, 1, 2,..., where σ denotes an eigenvalue approximation. In practice, the iterates are normalized after each step. The main hurdle emerging from (2.3) is that if σ is near to an eigenvalue, then A σm is almost singular, cf. [17]. For instance, for the Rayleigh quotient iteration it holds σ = λ(x (k) ). Solving such equations with multigrid methods is a hard task [21, 25, 26]. To overcome this difficulty, one can either bound away σ from the eigenvalues of (A, M), see [1]. Alternatively, one can solve (instead of (2.3)) a non-singular coarse grid correction equation within the orthogonal complement of the actual eigenvector approximation. The latter approach provides the basis for the direct multigrid eigensolver introduced in [7, 8]. In order to realize such a coarse grid projection on the orthogonal complement, the algorithm is closely related to nested iteration which serves to provide sufficiently accurate eigenfunctions for building reliable coarse grid projections. The resulting two-grid method maps a given iterate x having the Rayleigh quotient λ(x) to the new eigenvector approximations x. It is given by (2.4) x = Sx, (Smoothing step), d c = R(A λ(x)m) x, (Coarse grid projection of the residual), d c = Q c d c, (M-orthogonal projection), v c = (A c λ(x)m c ) 1 d c, (Solution of correction equation), x = x P Q c v c, (Prolongation and correction). Therein, the index c denotes coarse grid quantities. R is a restriction operator, P is a prolongation and Q c is the orthogonal projection operator to the M-orthogonal complement of the actual eigenvector approximation. The coarse grid problem (2.4) in the orthogonal complement is quite different from (2.3) and can be solved by a multigrid scheme for singular equations [8]. III. Eigensolvers using multigrid as a linear solver: Let us consider an alternative way of solving (2.3) without getting into trouble due to singularity of A σm. We can do this by ignoring the shift parameter, i.e. σ = 0. Additionally, we assume an eigenproblem for a coercive and self-adjoint partial differential operator so that A and M are symmetric positive definite matrices. Now we apply multigrid preconditioning for A in order to determine an approximate solution of the linear system (2.5) Ax (k+1) = λ(x (k) )Mx (k), k = 0, 1, 2,.... In comparison to (2.3), we have introduced the scaling constant λ(x (k) ) on the right hand side of (2.5) in order to achieve its stationarity in eigenvectors. Such a scaling prepares the approximate solution of (2.5) by means of preconditioning. The multigrid preconditioner B 1 is an approximate inverse of A which is assumed to be a symmetric positive definite operator and which is characterized by a spectral equivalence γ 0 (x, Bx) (x, Ax) γ 1 (x, Bx), x R n,

4 4 KLAUS NEYMEYR for positive constants γ 0 and γ 1. It is often more convenient to assume a properly scaled preconditioner in such a way that (2.6) I B 1 A A γ, for a constant γ [0, 1); therein A denotes the operator norm induced by A. This assumption is typical of multigrid and domain decomposition preconditioners. The best preconditioners satisfy (2.6) with γ bounded away from 1 independently of the mesh size and hence the number of unknowns [2]. The approximate solution of (2.5) for x = x (k) and by using B 1 as a preconditioner yields a new iterate x approximating x (k+1) (2.7) x = x B 1 (Ax λ(x)mx). The iteration (2.7) can be considered as the most simple eigensolver embodying the idea of multigrid as a linear solver. It can be interpreted as a (multigrid) preconditioned variant of inverse iteration (PINVIT) [19]. In principle, its behavior can be understood by the associated error propagation equation (2.8) x λ(x)a 1 Mx = (I B 1 A)(x λ(x)a 1 Mx), which maps the error of the initial iterate x λ(x)a 1 Mx to the error of the new iterate x λ(x)a 1 Mx by premultiplying the error propagation matrix I B 1 A. We can also derive from (2.8) that x λ(x)a 1 Mx r A 1 B 1, which shows the Lipschitz continuity in A 1 of the mapping from the set of admissible preconditioners, as characterized by (2.6), to the set of possible new iterates. The Lipschitz constant r equals the norm of the residual vector Ax λ(x)mx. Keeping these interpretations in mind, it is not surprising that the convergence estimates for (2.7) are closely related to inverse iteration. Sharp convergence estimates in terms of the Rayleigh quotients of the iterates are available. If the Rayleigh quotient λ(x) is located between the smallest eigenvalue λ 1 and the next larger eigenvalue λ 2 then it holds (2.9) λ(x ) λ 1 λ 2 λ(x ) ( γ + (1 γ) λ 1 λ 2 ) 2 λ(x) λ 1 λ 2 λ(x), see [14, 19]. This estimate can be applied recursively and guarantees linear convergence of the eigenvalue approximations λ(x (k) ) to the eigenvalue λ 1. The case γ = 0 corresponding to an exact solution of (2.5) leads to the convergence factor (λ 1 /λ 2 ) 2 of inverse iteration. For the other extreme case of poorest preconditioning, i.e. γ 1, the convergence factor tends to 1. We have assumed in (2.7) a symmetric positive definite preconditioner B 1. But there is no consent in the literature whether positive definite or indefinite preconditioners lead to more efficient eigensolvers, cf. [12]. This has also opened the question of how accurately (2.5) is to be solved in order to achieve reasonable convergence estimates [6, 15, 22, 24]. In any way the cubic convergence of the Rayleigh quotient iteration cannot be transferred to the preconditioned multigrid case. A significant acceleration of (2.7) can be achieved by minimizing the Rayleigh quotient with respect to a certain search subspace. The idea of the Locally Optimal Block Preconditioned Conjugate Gradient method (LOBPCG, [13]) is to determine the new eigenvalue/vector approximation as the smallest Ritz value/vector with respect to the threedimensional space (2.10) span{x, B 1 (Ax λ(x)mx), x },

5 SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY 5 where x denotes the iterate preceding x. To implement the latter iteration only one additional vector (compared to (2.7)) is to be stored (i.e. the old iterate) and few additional inner products are to be computed to determine the Rayleigh-Ritz projections. It has been observed that these additional costs are more than compensated by a characteristic speedup of convergence. Practically, the LOBPCG scheme behaves like a preconditioned conjugate gradient iteration for boundary value problems [13]. 3. NUMERICAL EXPERIMENTS All numerical experiments covered in this section treat the eigenproblem for the Laplacian on [0, 1] 2 with homogeneous Dirichlet boundary conditions as a simple model problem. See [5, 20] for adaptive multigrid eigensolvers on non-convex domains and with discontinuous coefficient functions. The eigenproblem is discretized by using linear finite elements on uniform triangle meshes with mesh sizes h l = 2 l. The initial grid consists of 9 interior nodes for l = 2. Computations are presented up to l = 12 corresponding to a grid with interior nodes. Table 1 lists the convergence history of RQMG, PINVIT and LOBPCG for the smallest eigenvalue and the corresponding residuals for l = 6 (3969 inner nodes). We use a unigrid version of RQMG, which is easy to code but which is computationally expensive, see [16] for the details. The coordinate relaxation scheme is applied to a black-white ordering of the finite element functions. The preconditioned eigensolvers are equipped with a multigrid V (2, 2)-cycle preconditioner working on the sequence of grids for l = 2,..., 6. Pre- and postsmooting is done by each 2 steps of the weighted Jacobi method. For all computations the initial function in each node (x 1, x 2 ) is given by x x 2 2. Table 1 contains the eigenvalue approximations, i.e. the Rayleigh quotients of the iterates, and the corresponding Euclidean norms of the residuals Ax λ(x)mx with x M = 1. RQMG and LOBPCG converge within 5 iterations to a six-digit accuracy, whereas more than 10 iterations of PINVIT are required to reach the same accuracy. All these iterations are guaranteed to generate a sequence of monotone decreasing Rayleigh quotients. Rayleigh quotients λ(x (k) ) Residuals Ax (k) λ(x (k) )Mx (k) k RQMG PINVIT LOBPCG RQMG PINVIT LOBPCG TABLE 1. RQMG, PINVIT and LOBPCG convergence history for l = 6. Mesh independence of the convergence of (2.7) is illustrated by Table 2. On each grid with h = 1/2 l, l = 4,..., 12, the iteration is started with the restriction of x x 2 2 to the nodes. Then 25 iterations of (2.7) with V (2, 2)-Jacobi preconditioning have been applied resulting in a final iterate x h. The final eigenvalue approximations λ(x h ) are listed in Table 2. They reflect that λ(x h ) λ 1 = O(h 2 ), whereas λ 1 denotes the smallest eigenvalue of the continuous problem. Moreover, the norms of the residuals (defined as above) are shown, too. Note that the M 1 norm of the residual is an upper estimate for the error of the eigenvalue approximations. As the condition number of M is bounded independently of h, the last column of Table 2 indicates the mesh independent convergence, since the final

6 6 KLAUS NEYMEYR Level l Grid Inner nodes λ(x h ) Ax h λ(x h )Mx h TABLE 2. Eigenvalues and residuals after 25 iterations of (2.7) for h l = 2 l, l = 4,..., 12. residuals do not deteriorate as the mesh size gets smaller. Let us finally mention that it is a nontrivial task to give a fair comparison (in terms of total computational costs and efficiency) of the direct multigrid eigensolver [7, 8] with those eigensolvers discussed above. This is the case since the direct multigrid scheme is inseparably connected with nested iteration in order to generate the required projection operators. Therefore, we only refer to [9], where the direct multigrid eigensolver is applied to the biharmonic plate equation. 4. SUMMARY A review has been given on three quite different multigrid techniques for solving the eigenvalue problem for a (coercive and self-adjoint) elliptic partial differential operator. Each of these solvers achieves optimal or quasi-optimal complexity for the partial eigenproblem, i.e. to compute only a small number of eigenpairs. These multigrid eigensolvers can be considered as the eigenproblem-counterparts of the celebrated, highly efficient multigrid solvers for boundary value problems. Eigensolvers using multigrid as a linear solver are conceptually very simple, and recent results allow a new theoretical understanding. They can be implemented easily: For the multigrid part one can adopt any existing program code realizing multigrid preconditioning. The remaining task is only to implement the linear algebra of (2.7) or (2.10). REFERENCES [1] R.E. Bank. Analysis of a multilevel inverse iteration procedure for eigenvalue problems. SIAM J. Numer. Anal., 19: , [2] J.H. Bramble, J.E. Pasciak, and J. Xu. Parallel multilevel preconditioners. Math. Comp., 55(191):1 22, [3] A. Brandt, S. McCormick, and J. Ruge. Multigrid methods for differential eigenproblems. SIAM J. Sci. Statist. Comput., 4(2): , [4] Z. Cai, J. Mandel, and S. McCormick. Multigrid methods for nearly singular linear equations and eigenvalue problems. SIAM J. Numer. Anal., 34: , [5] T. Friese, P. Deuflhard, and F. Schmidt. A multigrid method for the complex Helmholtz eigenvalue problem. In Eleventh International Conference on Domain Decomposition Methods (London, 1998), pages (electronic). DDM.org, Augsburg, [6] G.H. Golub and Q. Ye. Inexact inverse iteration for generalized eigenvalue problems. BIT, 40(4): , [7] W. Hackbusch. On the computation of approximate eigenvalues and eigenfunctions of elliptic operators by means of a multi-grid method. SIAM J. Numer. Anal., 16: , 1979.

7 SOLVING MESH EIGENPROBLEMS WITH MULTIGRID EFFICIENCY 7 [8] W. Hackbusch. Multi-grid methods and applications. Springer series in computational mathematics 4. Springer, Berlin, [9] W. Hackbusch and G. Hofmann. Results of the eigenvalue problem for the plate equation. Z. Angew. Math. Phys., 31(6): , [10] V. Heuveline and C. Bertsch. On multigrid methods for the eigenvalue computation of nonselfadjoint elliptic operators. East-West J. Numer. Math., 8(4): , [11] V. Heuveline and R. Rannacher. A posteriori error control for finite element approximations of elliptic eigenvalue problems. Adv. Comput. Math., 15(1-4): , [12] A.V. Knyazev. Preconditioned eigensolvers an oxymoron? Electron. Trans. Numer. Anal., 7: , [13] A.V. Knyazev. Toward the optimal preconditioned eigensolver: Locally optimal block preconditioned conjugate gradient method. SIAM J. Sci. Comp., 23: , [14] A.V. Knyazev and K. Neymeyr. A geometric theory for preconditioned inverse iteration. III: A short and sharp convergence estimate for generalized eigenvalue problems. Accepted for Linear Algebra Appl., [15] Y. Lai, K. Lin, and W. Lin. An inexact inverse iteration for large sparse eigenvalue problems. Numer. Linear Algebra Appl., 4: , [16] J. Mandel and S.F. McCormick. A multilevel variational method for Au = λbu on composite grids. J. Comput. Phys., 80: , [17] S.F. McCormick. A mesh refinement method for Ax = λbx. Math. Comp., 36(154): , [18] S.F. McCormick. Multilevel adaptive methods for elliptic eigenproblems: a two-level convergence theory. SIAM J. Numer. Anal., 31: , [19] K. Neymeyr. A geometric theory for preconditioned inverse iteration. I: Extrema of the Rayleigh quotient. Linear Algebra Appl., 322:61 85, [20] K. Neymeyr. A posteriori error estimation for elliptic eigenproblems. Numer. Linear Algebra Appl., 9: , [21] M.G. Neytcheva and P.S. Vassilevski. Preconditioning of indefinite and almost singular finite element elliptic equations. SIAM J. Sci. Comput., 19(5): , [22] Y. Notay. Convergence analysis of inexact Rayleigh quotient iterations. Technical report, Service de Métrologie Nucléaire, Université Libre de Bruxelles, [23] B.N. Parlett. The symmetric eigenvalue problem. Prentice Hall, Englewood Cliffs New Jersey, [24] P. Smit and M. Paardekooper. The effects of inexact solvers in algorithms for symmetric eigenvalue problems. Linear Algebra Appl., 287: , [25] J. Xu. A new class of iterative methods for nonselfadjoint or indefinite problems. SIAM J. Numer. Anal., 29: , [26] H. Yserentant. Preconditioning indefinite discretization matrices. Numer. Math., 54(6): , MATHEMATISCHES INSTITUT DER UNIVERSITÄT TÜBINGEN, AUF DER MORGENSTELLE 10, TÜBINGEN, GERMANY. address:

### EIGIFP: A MATLAB Program for Solving Large Symmetric Generalized Eigenvalue Problems

EIGIFP: A MATLAB Program for Solving Large Symmetric Generalized Eigenvalue Problems JAMES H. MONEY and QIANG YE UNIVERSITY OF KENTUCKY eigifp is a MATLAB program for computing a few extreme eigenvalues

### Using the Karush-Kuhn-Tucker Conditions to Analyze the Convergence Rate of Preconditioned Eigenvalue Solvers

Using the Karush-Kuhn-Tucker Conditions to Analyze the Convergence Rate of Preconditioned Eigenvalue Solvers Merico Argentati University of Colorado Denver Joint work with Andrew V. Knyazev, Klaus Neymeyr

### The Algebraic Multigrid Projection for Eigenvalue Problems; Backrotations and Multigrid Fixed Points. Sorin Costiner and Shlomo Ta'asan

The Algebraic Multigrid Projection for Eigenvalue Problems; Backrotations and Multigrid Fixed Points Sorin Costiner and Shlomo Ta'asan Department of Applied Mathematics and Computer Science The Weizmann

### On the Superlinear Convergence of MINRES. Valeria Simoncini and Daniel B. Szyld. Report January 2012

On the Superlinear Convergence of MINRES Valeria Simoncini and Daniel B. Szyld Report 12-01-11 January 2012 This report is available in the World Wide Web at http://www.math.temple.edu/~szyld 0 Chapter

### A Review of Preconditioning Techniques for Steady Incompressible Flow

Zeist 2009 p. 1/43 A Review of Preconditioning Techniques for Steady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/43 PDEs Review : 1984 2005 Update

### A Balancing Algorithm for Mortar Methods

A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz

### A STOKES INTERFACE PROBLEM: STABILITY, FINITE ELEMENT ANALYSIS AND A ROBUST SOLVER

European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,

### Numerical behavior of inexact linear solvers

Numerical behavior of inexact linear solvers Miro Rozložník joint results with Zhong-zhi Bai and Pavel Jiránek Institute of Computer Science, Czech Academy of Sciences, Prague, Czech Republic The fourth

### AN INEXACT INVERSE ITERATION FOR COMPUTING THE SMALLEST EIGENVALUE OF AN IRREDUCIBLE M-MATRIX. 1. Introduction. We consider the eigenvalue problem

AN INEXACT INVERSE ITERATION FOR COMPUTING THE SMALLEST EIGENVALUE OF AN IRREDUCIBLE M-MATRIX MICHIEL E. HOCHSTENBACH, WEN-WEI LIN, AND CHING-SUNG LIU Abstract. In this paper, we present an inexact inverse

### Multigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu

Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu Outline Introduction Standard Cascadic Multigrid method(cmg) Economical

### Conjugate gradient method. Descent method. Conjugate search direction. Conjugate Gradient Algorithm (294)

Conjugate gradient method Descent method Hestenes, Stiefel 1952 For A N N SPD In exact arithmetic, solves in N steps In real arithmetic No guaranteed stopping Often converges in many fewer than N steps

### the sum of two projections. Finally, in Section 5, we apply the theory of Section 4 to the case of nite element spaces. 2. Additive Algorithms for Par

ON THE SPECTRA OF SUMS OF ORTHOGONAL PROJECTIONS WITH APPLICATIONS TO PARALLEL COMPUTING PETTER E. BJRSTAD y AND JAN MANDEL z Abstract. Many parallel iterative algorithms for solving symmetric, positive

### 2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9

Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers

### PARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS

PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method

### Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations

Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg

### ANONSINGULAR tridiagonal linear system of the form

Generalized Diagonal Pivoting Methods for Tridiagonal Systems without Interchanges Jennifer B. Erway, Roummel F. Marcia, and Joseph A. Tyson Abstract It has been shown that a nonsingular symmetric tridiagonal

### On complex shifted Laplace preconditioners for the vector Helmholtz equation

On complex shifted Laplace preconditioners for the vector Helmholtz equation C. Vuik, Y.A. Erlangga, M.B. van Gijzen, C.W. Oosterlee, D. van der Heul Delft Institute of Applied Mathematics c.vuik@tudelft.nl

### A Posteriori Estimates for Cost Functionals of Optimal Control Problems

A Posteriori Estimates for Cost Functionals of Optimal Control Problems Alexandra Gaevskaya, Ronald H.W. Hoppe,2 and Sergey Repin 3 Institute of Mathematics, Universität Augsburg, D-8659 Augsburg, Germany

### A Balancing Algorithm for Mortar Methods

A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz

### EXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS

EXACT DE RHAM SEQUENCES OF SPACES DEFINED ON MACRO-ELEMENTS IN TWO AND THREE SPATIAL DIMENSIONS JOSEPH E. PASCIAK AND PANAYOT S. VASSILEVSKI Abstract. This paper proposes new finite element spaces that

### The parallel computation of the smallest eigenpair of an. acoustic problem with damping. Martin B. van Gijzen and Femke A. Raeven.

The parallel computation of the smallest eigenpair of an acoustic problem with damping. Martin B. van Gijzen and Femke A. Raeven Abstract Acoustic problems with damping may give rise to large quadratic

### Domain decomposition methods via boundary integral equations

Domain decomposition methods via boundary integral equations G. C. Hsiao a O. Steinbach b W. L. Wendland b a Department of Mathematical Sciences, University of Delaware, Newark, Delaware 19716, USA. E

### Aggregation Algorithms for K-cycle Aggregation Multigrid for Markov Chains

Aggregation Algorithms for K-cycle Aggregation Multigrid for Markov Chains by Manda Winlaw A research paper presented to the University of Waterloo in partial fulfillment of the requirements for the degree

### FINITE-DIMENSIONAL LINEAR ALGEBRA

DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup

### A fast method for the solution of the Helmholtz equation

A fast method for the solution of the Helmholtz equation Eldad Haber and Scott MacLachlan September 15, 2010 Abstract In this paper, we consider the numerical solution of the Helmholtz equation, arising

### Computational Linear Algebra

Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 3: Iterative Methods PD

### Convergence of The Multigrid Method With A Wavelet. Abstract. This new coarse grid operator is constructed using the wavelet

Convergence of The Multigrid Method With A Wavelet Coarse Grid Operator Bjorn Engquist Erding Luo y Abstract The convergence of the two-level multigrid method with a new coarse grid operator is studied.

### Adaptive Time Space Discretization for Combustion Problems

Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk

### ALGEBRAIC MULTILEVEL METHODS FOR GRAPH LAPLACIANS

The Pennsylvania State University The Graduate School Department of Mathematics ALGEBRAIC MULTILEVEL METHODS FOR GRAPH LAPLACIANS A Dissertation in Mathematics by Yao Chen Submitted in Partial Fulfillment

### Lecture 3: Review of Linear Algebra

ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak, scribe: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters,

### Conjugate Gradient (CG) Method

Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous

### ECS231: Spectral Partitioning. Based on Berkeley s CS267 lecture on graph partition

ECS231: Spectral Partitioning Based on Berkeley s CS267 lecture on graph partition 1 Definition of graph partitioning Given a graph G = (N, E, W N, W E ) N = nodes (or vertices), E = edges W N = node weights

### Linear Algebra Final Exam Study Guide Solutions Fall 2012

. Let A = Given that v = 7 7 67 5 75 78 Linear Algebra Final Exam Study Guide Solutions Fall 5 explain why it is not possible to diagonalize A. is an eigenvector for A and λ = is an eigenvalue for A diagonalize

### Largest Bratu solution, lambda=4

Largest Bratu solution, lambda=4 * Universiteit Utrecht 5 4 3 2 Department of Mathematics 1 0 30 25 20 15 10 5 5 10 15 20 25 30 Accelerated Inexact Newton Schemes for Large Systems of Nonlinear Equations

### Theoretical Bounds for Algebraic Multigrid Performance: Review and Analysis

NUMERICL LINER LGEBR WITH PPLICTIONS Numer. Linear lgebra ppl. 2013; 00:1 27 Published online in Wiley InterScience (www.interscience.wiley.com). Theoretical Bounds for lgebraic Multigrid Performance:

### arxiv: v1 [math.na] 16 Dec 2015

ANALYSIS OF A NEW HARMONICALLY ENRICHED MULTISCALE COARSE SPACE FOR DOMAIN DECOMPOSITION METHODS MARTIN J. GANDER, ATLE LONELAND, AND TALAL RAHMAN arxiv:52.05285v [math.na] 6 Dec 205 Abstract. We propose

### A Robust Approach to Minimizing H(div)-Dominated Functionals in an H 1 -Conforming Finite Element Space

NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 23 [Version: 2002/09/18 v1.02] A Robust Approach to Minimizing H(div)-Dominated Functionals in an H 1 -Conforming Finite

### Quasi-Newton Methods

Quasi-Newton Methods Werner C. Rheinboldt These are excerpts of material relating to the boos [OR00 and [Rhe98 and of write-ups prepared for courses held at the University of Pittsburgh. Some further references

### Numerische Mathematik

umer. Math. 73: 149 167 (1996) umerische Mathematik c Springer-Verlag 1996 Overlapping Schwarz methods on unstructured meshes using non-matching coarse grids Tony F. Chan 1, Barry F. Smith 2, Jun Zou 3

### Math 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.

Math 443 Differential Geometry Spring 2013 Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Endomorphisms of a Vector Space This handout discusses

### MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.

MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If

### ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS

J. OPERATOR THEORY 44(2000), 243 254 c Copyright by Theta, 2000 ADJOINTS, ABSOLUTE VALUES AND POLAR DECOMPOSITIONS DOUGLAS BRIDGES, FRED RICHMAN and PETER SCHUSTER Communicated by William B. Arveson Abstract.

### The Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment

he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State

### FOR many computational fluid dynamics (CFD) applications, engineers require information about multiple solution

AIAA Aviation 13-17 June 2016, Washington, D.C. 46th AIAA Fluid Dynamics Conference AIAA 2016-3809 Improved Functional-Based Error Estimation and Adaptation without Adjoints William C. Tyson *, Katarzyna

### 1. Introduction. In this work, we consider the solution of finite-dimensional constrained optimization problems of the form

SIAM J. SCI. COMPUT. Vol. 31, No. 6, pp. 4152 4175 c 2009 Society for Industrial and Applied Mathematics MULTILEVEL ALGORITHMS FOR LARGE-SCALE INTERIOR POINT METHODS MICHELE BENZI, ELDAD HABER, AND LAUREN

### A NEW METHOD FOR SOLVING ILL-CONDITIONED LINEAR SYSTEMS. Fazlollah Soleymani

Opuscula Math. 33, no. 2 (2013), 337 344 http://dx.doi.org/10.7494/opmath.2013.33.2.337 Opuscula Mathematica A NEW METHOD FOR SOLVING ILL-CONDITIONED LINEAR SYSTEMS Fazlollah Soleymani Communicated by

### Iterative Methods for Linear Systems

Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the

### RANA03-02 January Jacobi-Davidson methods and preconditioning with applications in pole-zero analysis

RANA03-02 January 2003 Jacobi-Davidson methods and preconditioning with applications in pole-zero analysis by J.Rommes, H.A. van der Vorst, EJ.W. ter Maten Reports on Applied and Numerical Analysis Department

### (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? Solution: dim N(A) 1, since rank(a) 3. Ax =

. (5 points) (a) If A is a 3 by 4 matrix, what does this tell us about its nullspace? dim N(A), since rank(a) 3. (b) If we also know that Ax = has no solution, what do we know about the rank of A? C(A)

### An Algorithmist s Toolkit September 15, Lecture 2

18.409 An Algorithmist s Toolkit September 15, 007 Lecture Lecturer: Jonathan Kelner Scribe: Mergen Nachin 009 1 Administrative Details Signup online for scribing. Review of Lecture 1 All of the following

### ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS

ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis

### Conjugate Gradient Tutorial

Conjugate Gradient Tutorial Prof. Chung-Kuan Cheng Computer Science and Engineering Department University of California, San Diego ckcheng@ucsd.edu December 1, 2015 Prof. Chung-Kuan Cheng (UC San Diego)

### Line Search Methods for Unconstrained Optimisation

Line Search Methods for Unconstrained Optimisation Lecture 8, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Generic

### Numerical Methods for PDE-Constrained Optimization

Numerical Methods for PDE-Constrained Optimization Richard H. Byrd 1 Frank E. Curtis 2 Jorge Nocedal 2 1 University of Colorado at Boulder 2 Northwestern University Courant Institute of Mathematical Sciences,

### AN HELMHOLTZ ITERATIVE SOLVER FOR THE THREE-DIMENSIONAL SEISMIC IMAGING PROBLEMS?

European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate, J. Périaux (Eds) c TU Delft, Delft The Netherlands, 2006 AN HELMHOLTZ ITERATIVE SOLVER FOR THE THREE-DIMENSIONAL

### Overlapping Schwarz Preconditioners for Spectral. Problem in H(curl)

Overlapping Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Technical Report TR2002-83 November 22, 2002 Department of Computer Science Courant Institute of Mathematical

### Wavelet-Based Numerical Homogenization for Scaled Solutions of Linear Matrix Equations

International Journal of Discrete Mathematics 2017; 2(1: 10-16 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.20170201.13 Wavelet-Based Numerical Homogenization for Scaled Solutions

### EE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2

EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 2 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 5, 2012 Andre Tkacenko

### HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)

HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe

### On the mathematical background of Google PageRank algorithm

Working Paper Series Department of Economics University of Verona On the mathematical background of Google PageRank algorithm Alberto Peretti, Alberto Roveda WP Number: 25 December 2014 ISSN: 2036-2919

### Homework set 4 - Solutions

Homework set 4 - Solutions Math 407 Renato Feres 1. Exercise 4.1, page 49 of notes. Let W := T0 m V and denote by GLW the general linear group of W, defined as the group of all linear isomorphisms of W

### Orthogonal iteration to QR

Week 10: Wednesday and Friday, Oct 24 and 26 Orthogonal iteration to QR On Monday, we went through a somewhat roundabout algbraic path from orthogonal subspace iteration to the QR iteration. Let me start

### Conjugate Gradient Method

Conjugate Gradient Method direct and indirect methods positive definite linear systems Krylov sequence spectral analysis of Krylov sequence preconditioning Prof. S. Boyd, EE364b, Stanford University Three

### Singular Value Decomposition

Singular Value Decomposition Motivatation The diagonalization theorem play a part in many interesting applications. Unfortunately not all matrices can be factored as A = PDP However a factorization A =

### A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators

A Framework for Analyzing and Constructing Hierarchical-Type A Posteriori Error Estimators Jeff Ovall University of Kentucky Mathematics www.math.uky.edu/ jovall jovall@ms.uky.edu Kentucky Applied and

### 1 Introduction and preliminaries

Proximal Methods for a Class of Relaxed Nonlinear Variational Inclusions Abdellatif Moudafi Université des Antilles et de la Guyane, Grimaag B.P. 7209, 97275 Schoelcher, Martinique abdellatif.moudafi@martinique.univ-ag.fr

### Properties of Matrices and Operations on Matrices

Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,

### Ph.D. Katarína Bellová Page 1 Mathematics 2 (10-PHY-BIPMA2) EXAM - Solutions, 20 July 2017, 10:00 12:00 All answers to be justified.

PhD Katarína Bellová Page 1 Mathematics 2 (10-PHY-BIPMA2 EXAM - Solutions, 20 July 2017, 10:00 12:00 All answers to be justified Problem 1 [ points]: For which parameters λ R does the following system

### Solutions of Selected Problems

1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative

### Variational Formulations

Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that

### 2. Review of Linear Algebra

2. Review of Linear Algebra ECE 83, Spring 217 In this course we will represent signals as vectors and operators (e.g., filters, transforms, etc) as matrices. This lecture reviews basic concepts from linear

### Key words. linear equations, deflation, eigenvalues, Lanczos, conjugate gradient, QCD, multiple right-hand sides, symmetric, Hermitian

DEFLATED AND RESTARTED SYMMETRIC LANCZOS METHODS FOR EIGENVALUES AND LINEAR EQUATIONS WITH MULTIPLE RIGHT-HAND SIDES ABDOU M. ABDEL-REHIM, RONALD B. MORGAN, DYWAYNE A. NICELY, AND WALTER WILCOX Abstract.

### THE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE FOURTH ORDER PROBLEM *

Journal of Computational Mathematics Vol.30, No.5, 2012, 449 460. http://www.global-sci.org/jcm doi:10.4208/jcm.1203-m3855 THE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE

### Maximum-weighted matching strategies and the application to symmetric indefinite systems

Maximum-weighted matching strategies and the application to symmetric indefinite systems by Stefan Röllin, and Olaf Schenk 2 Technical Report CS-24-7 Department of Computer Science, University of Basel

### arxiv: v1 [math.na] 1 May 2013

arxiv:3050089v [mathna] May 03 Approximation Properties of a Gradient Recovery Operator Using a Biorthogonal System Bishnu P Lamichhane and Adam McNeilly May, 03 Abstract A gradient recovery operator based

### PRECONDITIONING ORBITAL MINIMIZATION METHOD FOR PLANEWAVE DISCRETIZATION 1. INTRODUCTION

PRECONDITIONING ORBITAL MINIMIZATION METHOD FOR PLANEWAVE DISCRETIZATION JIANFENG LU AND HAIZHAO YANG ABSTRACT. We present an efficient preconditioner for the orbital minimization method when the Hamiltonian

### Orthogonal iteration revisited

Week 10 11: Friday, Oct 30 and Monday, Nov 2 Orthogonal iteration revisited Last time, we described a generalization of the power methods to compute invariant subspaces. That is, starting from some initial

### Enforcing Passivity for Admittance Matrices Approximated by Rational Functions

IEEE TRANSACTIONS ON POWER SYSTEMS, VOL. 16, NO. 1, FEBRUARY 2001 97 Enforcing Passivity for Admittance Matrices Approximated by Rational Functions Bjørn Gustavsen, Member, IEEE and Adam Semlyen, Life

### Multigrid Solution of the Debye-Hückel Equation

Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections 5-17-2013 Multigrid Solution of the Debye-Hückel Equation Benjamin Quanah Parker Follow this and additional works

### Math 25a Practice Final #1 Solutions

Math 25a Practice Final #1 Solutions Problem 1. Suppose U and W are subspaces of V such that V = U W. Suppose also that u 1,..., u m is a basis of U and w 1,..., w n is a basis of W. Prove that is a basis

### Introduction to Optimization

Introduction to Optimization Gradient-based Methods Marc Toussaint U Stuttgart Gradient descent methods Plain gradient descent (with adaptive stepsize) Steepest descent (w.r.t. a known metric) Conjugate

### SELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY

Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS

### An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems

An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A. We present an efficient implicit

### Large automorphism groups of 16-dimensional planes are Lie groups

Journal of Lie Theory Volume 8 (1998) 83 93 C 1998 Heldermann Verlag Large automorphism groups of 16-dimensional planes are Lie groups Barbara Priwitzer, Helmut Salzmann Communicated by Karl H. Hofmann

### Technische Universität Graz

Technische Universität Graz Schur complement preconditioners for the biharmonic Dirichlet boundary value problem L. John, O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2013/4

### Math 18, Linear Algebra, Lecture C00, Spring 2017 Review and Practice Problems for Final Exam

Math 8, Linear Algebra, Lecture C, Spring 7 Review and Practice Problems for Final Exam. The augmentedmatrix of a linear system has been transformed by row operations into 5 4 8. Determine if the system

### An error estimate for matrix equations

Applied Numerical Mathematics 50 (2004) 395 407 www.elsevier.com/locate/apnum An error estimate for matrix equations Yang Cao, Linda Petzold 1 Department of Computer Science, University of California,

### MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T KEMP) Definition 01 If T (x) = Ax is a linear transformation from R n to R m then Nul (T ) = {x R n : T (x) = 0} = Nul (A) Ran (T ) = {Ax R m : x R n } = {b R m

### Singular Value Decomposition

Chapter 6 Singular Value Decomposition In Chapter 5, we derived a number of algorithms for computing the eigenvalues and eigenvectors of matrices A R n n. Having developed this machinery, we complete our

### Math 2331 Linear Algebra

5. Eigenvectors & Eigenvalues Math 233 Linear Algebra 5. Eigenvectors & Eigenvalues Shang-Huan Chiu Department of Mathematics, University of Houston schiu@math.uh.edu math.uh.edu/ schiu/ Shang-Huan Chiu,

### Today: eigenvalue sensitivity, eigenvalue algorithms Reminder: midterm starts today

AM 205: lecture 22 Today: eigenvalue sensitivity, eigenvalue algorithms Reminder: midterm starts today Posted online at 5 PM on Thursday 13th Deadline at 5 PM on Friday 14th Covers material up to and including

### Stability of Krylov Subspace Spectral Methods

Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot

### FINDING PARALLELISM IN GENERAL-PURPOSE LINEAR PROGRAMMING

FINDING PARALLELISM IN GENERAL-PURPOSE LINEAR PROGRAMMING Daniel Thuerck 1,2 (advisors Michael Goesele 1,2 and Marc Pfetsch 1 ) Maxim Naumov 3 1 Graduate School of Computational Engineering, TU Darmstadt

### 5 Quiver Representations

5 Quiver Representations 5. Problems Problem 5.. Field embeddings. Recall that k(y,..., y m ) denotes the field of rational functions of y,..., y m over a field k. Let f : k[x,..., x n ] k(y,..., y m )

### Algorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method

Algorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method Ilya B. Labutin A.A. Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 3, acad. Koptyug Ave., Novosibirsk

### Math 5520 Homework 2 Solutions

Math 552 Homework 2 Solutions March, 26. Consider the function fx) = 2x ) 3 if x, 3x ) 2 if < x 2. Determine for which k there holds f H k, 2). Find D α f for α k. Solution. We show that k = 2. The formulas

### IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each

### SOR as a Preconditioner. A Dissertation. Presented to. University of Virginia. In Partial Fulllment. of the Requirements for the Degree

SOR as a Preconditioner A Dissertation Presented to The Faculty of the School of Engineering and Applied Science University of Virginia In Partial Fulllment of the Reuirements for the Degree Doctor of