Abstract:Two innovations are presented to improve the orthogonal subspace iteration. Firstly, a refined strategy is proposed that replaces conventional Ritz vectors by refined Ritz vectors in subspace iteration. Each refined Ritz vector is chosen such that the norm of the residual formed with the Ritz value is minimized over the subspace involved, and it can be cheaply computed by solving a small problem. Secondly, a standard updating matrix, Ritz vector or Schur vector matrix, is replaced by a new one that is generated only by the refined Ritz vectors used to approximate the wanted eigenvectors. A combination of them yields a refined subspace iteration algorithm. A qualitative analysis is given for the refined algorithm, showing that it may be considerably more efficient than a conventional subspace iteration algorithm. When the subspace dimension is bigger than the number of the wanted eigenpairs, numerical examples show that the refined algorithm has a sharp superiority. Numerical comparisons are also drawn for the refined algorithm with the popular package ARPACK and an implicitly restarted version of the refined Arnoldi method. The algorithm is suited for computing a few dominant eigenpairs of large sparse matrices since the only action of a matrix A in question into the program is a subroutine to form the product A Q, where Q is a matrix having much less columns than rows. (C) 2000 IMACS. Published by Elsevier Science B.V. All rights reserved.

A Modified Inverse Iteration Method for Computing the Symmetric Tridiagonal Eigenvectors

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

A Novel Divisional Bisection Method for the Symmetric Tridiagonal Eigenvalue Problem

COMPUTING EIGENVECTOR OF SYMMETRIC TRIDIAGONAL MATRIX BY USING SIMPLE INVERSE ITERATION

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

Triangularized Orthogonalization-Free Method for Solving Extreme Eigenvalue Problems

A Ulm-like Method for Inverse Eigenvalue Problems

Fast parallelizable methods for the Hermitian eigenvalue problem

A new method for the computation of eigenvalues

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Computing Eigenvectors of Normal Matrices with Simple Inverse Iteration

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

An Improved Two‐Step Method for Inverse Eigenvalue Problems

Modified CRI iteration methods for complex symmetric indefinite linear systems

Modified Newton'S Algorithm for Computing the Group Inverses of Singular Toeplitz Matrices

A Class of Inverse Problem with Symmetric Skew Anti-symmetric Tridiagonal Matrices

A Novel Iterative Method for Computing Generalized Inverse

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

Fast Parallelizable Methods for Computing Invariant Subspaces of Hermitian Matrices