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.

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

A Refined Iterative Algorithm Based on the Block Arnoldi Process for Large Unsymmetric Eigenproblems

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

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

A Refined Variant of the Inverse-Free Krylov Subspace Method for Symmetric Generalized Eigenvalue Problems

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

A Refined Arnoldi Type Method for Large Scale Eigenvalue Problems

Improving eigenvectors in Arnoldi's method

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

Comparison of Subspace Iteration, Iterative Ritz Vector Method and Iterative Lanczos Method

Composite orthogonal projection methods for large matrix eigenproblems

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

Residuals of refined projection methods for large matrix eigenproblems

An Analysis of the Rayleigh-Ritz Method for Approximating Eigenspaces

Subspace Methods with Local Refinements for Eigenvalue Computation Using Low-Rank Tensor-Train Format.

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

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

A Semiorthogonal Generalized Arnoldi Method And Its Variations For Quadratic Eigenvalue Problems

Convergence of Q-Ritz Pairs, Refined Q-Ritz Vectors and Q-Rayleigh-ritz Method for Quadratic Eigenvalue Problems

Iterative Methods by Space Decomposition and Subspace Correction

Advances in Solution of Classical Generalized Eigenvalue Problem