The convergence of Jacobi–Davidson iterations for Hermitian eigenproblems
✍ Scribed by Jasper van den Eshof
- Publisher
- John Wiley and Sons
- Year
- 2002
- Tongue
- English
- Weight
- 140 KB
- Volume
- 9
- Category
- Article
- ISSN
- 1070-5325
- DOI
- 10.1002/nla.266
No coin nor oath required. For personal study only.
✦ Synopsis
Abstract
Rayleigh quotient iteration is an iterative method with some attractive convergence properties for finding (interior) eigenvalues of large sparse Hermitian matrices. However, the method requires the accurate (and, hence, often expensive) solution of a linear system in every iteration step. Unfortunately, replacing the exact solution with a cheaper approximation may destroy the convergence. The (Jacobi‐) Davidson correction equation can be seen as a solution for this problem. In this paper we deduce quantitative results to support this viewpoint and we relate it to other methods. This should make some of the experimental observations in practice more quantitative in the Hermitian case. Asymptotic convergence bounds are given for fixed preconditioners and for the special case if the correction equation is solved with some fixed relative residual precision. A dynamic tolerance is proposed and some numerical illustration is presented. Copyright © 2002 John Wiley & Sons, Ltd.
📜 SIMILAR VOLUMES
## Abstract General stationary iterative methods with a singular matrix __M__ for solving range‐Hermitian singular linear systems are presented, some convergence conditions and the representation of the solution are also given. It can be verified that the general Ortega–Plemmons theorem and Keller
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
Ä 4 Ä 4 Ä 4 mapping. Given a sequence x in D and two real sequences t and s Ä 4 5 5 we prove that if x is bounded, then lim Tx y x s 0. The conditions on n n ª ϱn n D , X, and T are shown which guarantee the weak and strong convergence of the Ishikawa iteration process to a fixed point of T.