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 Note on Wilkinson's Iterative Refinement of Solution with Automatic Step-Size Control

Wilkinson’s Iterative Refinement of Solution with Automatic Step-Size Control for Linear System of Equations

Extended version with the analysis of dynamic system for iterative refinement of solution

A New Iterative Refinement with Roundoff Error Analysis.

Iterative Refinement of Solution with Biparameter for Solving Ill-Conditioned Systems of Linear Algebraic Equations

Stable iterative refinement algorithms for solving linear systems

New Iterative Improvement of a Solution for an Ill-conditioned System of Linear Equations Based on a Linear Dynamic System

Iterative Refinement of Computing Inverse Matrix

A Ulm-like Method for Inverse Eigenvalue Problems

On iterative solvers of the matrix equation X+A/sup T/X/sup -1/A=I

Randomized Iterative Solver as Iterative Refinement: A Simple Fix Towards Backward Stability

Control-theoretic Design of Iterative Methods for Symmetric Linear Systems of Equations

Developing iterative algorithms to solve Sylvester tensor equations

Note on Predictor-Corrector Process for Ill-Conditioned Linear System of Equations

Note on the Improvement of Newton’s Method for System of Nonlinear Equations

On the Iterative Algorithm for Large Sparse Saddle Point Problems

Mixed precision iterative refinement for linear inverse problems

Mixed precision iterative refinement for least squares with linear equality constraints and generalized least squares problems

On iterative solvers of the matrix equation X+ATX-1 A=I

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

On Convergence of the Inexact Rayleigh Quotient Iteration with the Lanczos Method Used for Solving Linear Systems