Asymptotic Convergence of Conjugate Gradient Methods for the Partial Symmetric Eigenproblem
β Scribed by Luca Bergamaschi; Giuseppe Gambolati; Giorgio Pini
- Publisher
- John Wiley and Sons
- Year
- 1997
- Tongue
- English
- Weight
- 149 KB
- Volume
- 4
- Category
- Article
- ISSN
- 1070-5325
No coin nor oath required. For personal study only.
β¦ Synopsis
Recently an efficient method (DACG) for the partial solution of the symmetric generalized eigenproblem Ax = Ξ»Bx has been developed, based on the conjugate gradient (CG) minimization of the Rayleigh quotient over successive deflated subspaces of decreasing size. The present paper provides a numerical analysis of the asymptotic convergence rate Ο j of DACG in the calculation of the eigenpair Ξ» j , u j , when the scheme is preconditioned with A -1 . It is shown that, when the search direction are A-conjugate, Ο j is well approximated by 4/ΞΎ j , where ΞΎ j is the Hessian condition number of a Rayleigh quotient defined in appropriate oblique complements of the space spanned by the leftmost eigenvectors u 1 , u 2 , . . ., u j-1 already calculated. It is also shown that 1/ΞΎ j is equal to the relative separation between the eigenvalue Ξ» j currently sought and the next higher one Ξ» j +1 . A modification of DACG (MDACG) is studied, which involves a new set of CG search directions which are made M-conjugate, with M a matrix approximating the Hessian. By distinction, MDACG has an asymptotic rate of convergence which appears to be inversely proportional to the square root of ΞΎ j , in complete agreement with the theoretical results known for the CG solution to linear systems.
π SIMILAR VOLUMES
In this paper, a detailed description of CG for evaluating eigenvalue problems by minimizing the Rayleigh quotient is presented from both theoretical and computational viewpoints. Three variants of CG together with their asymptotic behaviours and restarted schemes are discussed. In addition, it is s
## Abstract This paper is concerned with proving theoretical results related to the convergence of the conjugate gradient (CG) method for solving positive definite symmetric linear systems. Considering the inverse of the projection of the inverse of the matrix, new relations for ratios of the __A__
Parallel preconditioners are considered for improving the convergence rate of the conjugate gradient method for solving sparse symmetric positive definite systems generated by finite element models of subsurface flow. The difficulties of adapting effective sequential preconditioners to the parallel
The problem of two-dimensional scattering of elastic waves by an elastic inclusion can be formulated in terms of a domain integral equation, in which the grad-div operator acts on a vector potential. The vector potential is the spatial convolution of a Green's function with the product of the densit