Abstract:We give a quantitative analysis of the Shift-Invert Residual Arnoldi (SIRA) method and the Jacobi–Davidson (JD) method for computing a simple eigenvalue nearest to a target σ and/or the associated eigenvector. In SIRA and JD, subspace expansion vectors at each step are obtained by solving certain (different) inner linear systems, respectively. We show that (i) SIRA and the JD method with the fixed target σ are mathematically equivalent when the inner linear systems are solved exactly and (ii) the inexact SIRA is asymptotically equivalent to the JD method when the inner linear systems in them are solved with the same accuracy. Remarkably, we prove that the inexact SIRA and JD methods mimic the exact SIRA well provided that the inner linear systems are iteratively solved with a fixed low or modest accuracy. It is opposed 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 is only solved with decreasing accuracy after the approximate eigenpair starts converging. We also show that SIRA and JD expand subspaces in a computationally optimal way. We propose restarted SIRA and JD algorithms and design practical stopping criteria for inner solvers. Numerical experiments confirm our theory and the considerable superiority of the (nonrestarted and restarted) inexact SIRA and JD to the inexact SIA, and demonstrate that the inexact SIRA and JD are similarly effective and mimic the exact SIRA very well.

IMPLICITLY RESTARTED GENERALIZED SECOND-ORDER ARNOLDI TYPE ALGORITHMS FOR THE QUADRATIC EIGENVALUE PROBLEM

A Refined Second-order Arnoldi (RSOAR) Method for the Quadratic Eigenvalue Problem and Implicit Restarting

SOAR: A Second-order Arnoldi Method for the Solution of the Quadratic Eigenvalue Problem.

Restarted Generalized Second-Order Krylov Subspace Methods for Solving Quadratic Eigenvalue Problems

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

Polynomial characterizations of the approximate eigenvectors by the refined Arnoldi method and an implicitly restarted refined Arnoldi algorithm

Stability Analysis of the Two-level Orthogonal Arnoldi Procedure

Second-order Krylov Subspace and Arnoldi Procedure

Randomized Implicitly Restarted Arnoldi method for the non-symmetric eigenvalue problem

The Second-Order Biorthogonalization Procedure and Its Application to Quadratic Eigenvalue Problems

A Refined Harmonic Rayleigh-Ritz Procedure and an Explicitly Restarted Refined Harmonic Arnoldi Algorithm

On IGMRES: an Incomplete Generalized Minimal Residual Method for Large Unsymmetric Linear Systems

Improving eigenvectors in Arnoldi's method

The Refined Harmonic Arnoldi Method and an Implicitly Restarted Refined Algorithm for Computing Interior Eigenpairs of Large Matrices

A VARIANT OF THE REFINED ARNOLDI METHOD FOR LARGE SCALE MATRIX EIGENPROBLEMS

Subspace Expansion in the Shift-Invert Residual Arnoldi Method and the Jacobi--Davidson Method: Theory and Algorithms

Block GMRES Method with Inexact Breakdowns and Deflated Restarting

A Refined Shift-And-Invert Arnoldi Algorithm for Large Unsymmetric Generalized Eigenproblems

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

Krylov Subspace Methods for the Generalized Sylvester Equation.

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems