𝔖 Bobbio Scriptorium
✦   LIBER   ✦

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


CONJUGATE GRADIENT METHODS FOR SOLVING T
✍ Y. T. FENG; D. R. J. OWEN πŸ“‚ Article πŸ“… 1996 πŸ› John Wiley and Sons 🌐 English βš– 1000 KB

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

New results on the convergence of the co
✍ R. Bouyouli; G. Meurant; L. Smoch; H. Sadok πŸ“‚ Article πŸ“… 2009 πŸ› John Wiley and Sons 🌐 English βš– 120 KB

## 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 Multigrid Preconditioning of th
✍ Leesa Brieger; Giuditta Lecca πŸ“‚ Article πŸ“… 1998 πŸ› Elsevier Science 🌐 English βš– 186 KB

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

A Weak Form of the Conjugate Gradient FF
✍ George Pelekanos; Ralph E. Kleinman; Peter M. van den Berg πŸ“‚ Article πŸ“… 2000 πŸ› Elsevier Science 🌐 English βš– 144 KB

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