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.

Smoothing algorithm for the maximal eigenvalue of non-defective positive matrices

An Efficient Algorithm for Finding the Maximal Eigenvalue of Zero Symmetric Nonnegative Matrices

On a minimization problem of the maximum generalized eigenvalue: properties and algorithms

An Improved Diagonal Transformation Algorithm for the Maximum Eigenvalue of Zero Symmetric Nonnegative Matrices

A Smoothing-Type Algorithm for Solving System of Inequalities

Linear Convergence of an Algorithm for Computing the Largest Eigenvalue of a Nonnegative Tensor.

Analysis Of Max-Min Eigenvalue Of Constrained Linear Combinations Of Symmetric Matrices

Optimal correcting of absolute value equations by using smoothing techniques

An Enhanced Formulation Of Error Bound In Subspace Iteration Method

A High Accuracy Post-processing Algorithm for the Eigenvalues of Elliptic Operators

A Note on Improvement on Bounds for Nonmaximal Eigenvalues of Symmetric Positive Matrices

Towards positive definite solutions of a class of nonlinear matrix equations.

FAST ALGORITHMS FOR DETERMINING POSITIVE DEFINITENESS OF LARGE SCALE MATRICES IN SOLID STRUCTURE ANALYSIS1)

A New Preconditioner for Solving Weighted Toeplitz Least Squares Problems

Accurate computation of the smallest eigenvalue of a diagonally dominant M-matrix.

Linear System Based Approach for Solving Some Related Problems of M-matrices

A Smoothing Approximation Algorithm for Nonlinear Complementarity Problems

On Hermitian positive definite solution of nonlinear matrix equation X plus A* X-2 A = Q

Smoothing Algorithm with Edge-preserving by Extrema Constraints

Efficient Algorithms for Solving Condition Number-Constrained Matrix Minimization Problems

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