The Shift-Invert Residual Arnoldi Method and the Jacobi-Davidson Method: Theory and Algorithms ∗

Zhongxiao Jia,Cen Li
2011-01-01
Abstract:Using a new analysis approach, we establish a general convergence theory of the ShiftInvert Residual Arnoldi (SIRA) method for computing a simple eigenvalue nearest to a target σ and/or the associated eigenvector. In SIRA, a subspace expansion vector at each step is obtained by solving a certain inner linear system. We prove that the inexact SIRA method mimic the exact SIRA well provided that the inner linear systems are iteratively solved with low or modest accuracy. The results demonstrate that SIRA is superior to the inexact Shift-Invert Arnoldi (SIA) method, where the inner linear system involved must be solved with very high accuracy whenever the approximate eigenpair is of poor accuracy and are only solved with decreasing accuracy after the approximate eigenpair starts converging. Based on the theory, we design a practical stopping criterion for inner solvers. For practical purpose, we propose a restarted SIRA algorithm. Our analysis approach applies to the JD method as well, a general convergence theory is obtained similarly for the standard Jacobi–Davidson (JD) method with a fixed target, and a practical JD algorithm is developed. Numerical experiments confirm our theory and the considerable superiority of the (non-restarted and restarted) inexact SIRA to the inexact SIA, and demonstrate that the inexact SIRA and JD are similarly effective and mimic the exact SIRA very well.
What problem does this paper attempt to address?