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.

Improved upper bounds for the structured singular value

Structured Singular Value of a Repeated Complex Full-Block Uncertainty

Parallel Branch and Bound Algorithm for Computing Maximal Structured Singular Value

Structured singular value of a repeated complex full‐block uncertainty

Structured singular values and their application in computing eigenvalue backward errors of the Rosenbrock system matrix

Exact Solution for the Rank-One Structured Singular Value with Repeated Complex Full-Block Uncertainty

An Extended Ulm-like Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

Small singular values can increase in lower precision

On the Robust Stability Analysis with Real Block Structured Uncertainties.

Towards a more robust algorithm for computing the restricted singular value decomposition

A note on Oishi's lower bound for the smallest singular value of linearized Galerkin equations

A new lower bound for the smallest singular value

A parallel frequency sweeping algorithm for computing maximal structured singular value μ1

A lower bound for the smallest singular value of a matrix

Verified Error Bounds for Matrix Decompositions

Matrix perturbation analysis of methods for extracting singular values from approximate singular subspaces

New rigorous perturbation bounds for the Cholesky-like factorization of skew-symmetric matrix

Schur Stability Of Uncertain Matrices

An Enhanced Formulation Of Error Bound In Subspace Iteration Method

Estimation to structured distances to singularity for matrix pencils with symmetry structures: A linear algebra-based approach

On the Condition Number of the Total Least Squares Problem