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 Refined Arnoldi Type Method for Large Scale Eigenvalue Problems

A VARIANT OF THE REFINED ARNOLDI METHOD FOR LARGE SCALE MATRIX EIGENPROBLEMS

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

A Refined Iterative Algorithm Based on the Block Arnoldi Process for Large Unsymmetric Eigenproblems

Improving eigenvectors in Arnoldi's method

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

The Refined Harmonic Arnoldi Method and an Implicitly Restarted Refined Algorithm for Computing Interior Eigenpairs of Large Matrices

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

A Ulm-like Method for Inverse Eigenvalue Problems

A Global Harmonic Arnoldi Method for Large Non-Hermitian Eigenproblems with an Application to Multiple Eigenvalue Problems

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

ARNOLDI TYPE ALGORITHMS FOR LARGE UNSYMMETRIC MULTIPLE EIGENVALUE PROBLEMS

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Polynomial characterizations of the approximate eigenvectors by the refined Arnoldi method and an implicitly restarted refined Arnoldi algorithm

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

The Rayleigh–Ritz method, refinement and Arnoldi process for periodic matrix pairs

A new iterative method for simultaneous computation of several eigenpair derivatives of a large matrix

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

Arnoldi Reduction Algorithm for Large Scale Gyroscopic Eigenvalue Problem

An Analysis of the Rayleigh-Ritz Method for Approximating Eigenspaces

Restarted Generalized Second-Order Krylov Subspace Methods for Solving Quadratic Eigenvalue Problems