Abstract:In subspace iteration method (SIM), the relative difference of approximated eigenvalues between two consecutive iterations is usually employed as the convergence criterion. However, though it controls the convergence of eigenvalues well, it cannot guarantee the convergence of eigenvectors in all cases. In the case when there is no shifting, the best choice for the convergence criterion of eigenvalues may generally be the computable error bound proposed by Matthies, which is based on an estimation of Rayleigh quotient of approximated eigenvalues expressed in the subspace. Matthies' form guarantees the convergence of both eigenvalues and eigenvectors and can be computed with almost negligible operations. However, it is not as popular as expected in implementations, partly because it does not consider the popular shifting acceleration technique of subspace iterations. In this paper, we extend Matthies' form to the case of nonzero shifting and prove that this extended error bound form with nonzero shifting can be used generally as a convergence criterion for eigenpairs. Besides, this paper details the derivation to illustrate that the extended error bound can also be applied to the case of positive semi-definite mass matrix by only slightly modifying the subspace iteration procedure. Numerical tests are presented to illustrate the motivation and to demonstrate the better performance of the modified computable error bound. The studies in this paper indicate that the modified Matthies' form of error bound can be effectively used as a preferred convergence criterion in the SIM. Copyright (C) 2009 John Wiley & Sons, Ltd.

Molera, JM, see Dopico, FM 309 (2000) 3 Moro, J., see Dopico, FM 309 (2000) 3 Oliveira, S. and DE Stewart, Exponential splittings of products of matrices and accu-rately …

Perturbation of the Drazin Inverse for Matrices with Equal Eigenprojections at Zero

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

Mixed methods and lower eigenvalue bounds

Matrix‐less methods for the spectral approximation of large non‐Hermitian Toeplitz matrices: A concise theoretical analysis and a numerical study

Ju l 2 01 9 Eigenvalue estimates via Hömander ’ s L 2-method

Avoiding discretization issues for nonlinear eigenvalue problems

Spectral Multipliers II: Elliptic and Parabolic Operators and Bochner-Riesz Means

Error Bounds for Perturbation of the Drazin Inverse of Closed Operators with Equal Spectral Projections

A smoothing spline algorithm to interpolate and predict the eigenvalues of matrices extracted from the sequence of preconditioned banded symmetric Toeplitz matrices

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop theory

Numerical optimisation of Dirac eigenvalues

Wegner estimate and upper bound on the eigenvalue condition number of non-Hermitian random matrices

Relative perturbation bounds for the eigenvalues of diagonalizable and singular matrices – Application of perturbation theory for simple invariant subspaces

Fast computation of spectral densities for generalized eigenvalue problems

Eigenvalue superposition for Toeplitz matrix-sequences with matrix order dependent symbols

An Enhanced Formulation Of Error Bound In Subspace Iteration Method

On Uncertainty Quantification of Eigenvalues and Eigenspaces with Higher Multiplicity

Bounds for the Relative and Absolute Spectral Variations of Matrices