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.

Refined projection method for the quadratic eigenvalue problem and the corresponding Jacobi-Davidson method

A Refined Jacobi-Davidson Method for the Quadratic Eigenvalue Problem

A Refined Jacobi-Davidson Method and Its Correction Equation

Residuals of refined projection methods for large matrix eigenproblems

Composite orthogonal projection methods for large matrix eigenproblems

A Theoretical Comparison of Two Classes of Projection Methods

A Ulm-like Method for Inverse Eigenvalue Problems

On the Convergence of Ritz Pairs and Refined Ritz Vectors for Quadratic Eigenvalue Problems

Convergence of Q-Ritz Pairs, Refined Q-Ritz Vectors and Q-Rayleigh-ritz Method for Quadratic Eigenvalue Problems

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

Some Theoretical Comparisons of Refined Ritz Vectors and Ritz Vectors

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

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

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

A modification of the Jacobi-Davidson method

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

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

A Refined Second-order Arnoldi (RSOAR) Method for the Quadratic Eigenvalue Problem and Implicit Restarting

Generalized Jacobi Method for Computing Eigenvalues of Dual Quaternion Hermitian Matrices

A Projection Method Based on Extended Krylov Subspaces for Solving Sylvester Equations

A numerical projection technique for large-scale eigenvalue problems