Inner iterations in the shift-invert residual Arnoldi method and the Jacobi-Davidson method

From MaRDI portal
Publication:477157

DOI10.1007/S11425-014-4791-5zbMATH Open1312.65054arXiv1109.5455OpenAlexW2118497999MaRDI QIDQ477157FDOQ477157

Zhongxiao Jia, Cen Li

Publication date: 2 December 2014

Published in: Science China. Mathematics (Search for Journal in Brave)

Abstract: Using a new analysis approach, we establish a general convergence theory of the Shift-Invert Residual Arnoldi (SIRA) method for computing a simple eigenvalue nearest to a given target sigma and 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 mimics the exact SIRA well, that is, the former uses almost the same outer iterations to achieve the convergence as the latter does if all the inner linear systems are iteratively solved with {em low} or {em modest} accuracy during outer iterations. Based on the theory, we design practical stopping criteria for inner solves. Our analysis is on one step expansion of subspace and the approach applies to the Jacobi--Davidson (JD) method with the fixed target sigma as well, and a similar general convergence theory is obtained for it. Numerical experiments confirm our theory and demonstrate that the inexact SIRA and JD are similarly effective and are considerably superior to the inexact SIA.


Full work available at URL: https://arxiv.org/abs/1109.5455




Recommendations




Cites Work


Cited In (19)

Uses Software





This page was built for publication: Inner iterations in the shift-invert residual Arnoldi method and the Jacobi-Davidson method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q477157)