Variable accuracy of matrix-vector products in projection methods for eigencomputation

摘要

We analyze the behavior of projection-type schemes, such as the Arnoldi and Lanczos methods, for the approximation of a few eigenvalues and eigenvectors of a matrix A, when A cannot be applied exactly but only with a possibly large perturbation. This occurs, for instance, in shift-and-invert procedures or when dealing with large generalized eigenvalue problems. We theoretically show that the accuracy with which A is applied at each iteration can be relaxed, as convergence to specific eigenpairs takes place. We show that the size of the perturbation is allowed to be inversely proportional to the current residual norm, in a way that also depends on the sensitivity of the matrix A. This result provides a complete understanding of reported experimental evidence in the recent literature. Moreover, we adapt our theoretical results to devise a practical relaxation criterion to achieve convergence of the inexact procedure. Numerical experiments validate our analysis.