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.

Second-order Krylov Subspace and Arnoldi Procedure

SOAR: A Second-order Arnoldi Method for the Solution of the Quadratic Eigenvalue Problem.

The Second-Order Biorthogonalization Procedure and Its Application to Quadratic Eigenvalue Problems

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

Model-order Reduction of Large-Scale Kth-Order Linear Dynamical Systems Via a Kth-Order Arnoldi Method

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

The Arnoldi Process, Short Recursions, and Displacement Ranks

A Two-Directional Arnoldi Process and Its Application to Parametric Model Order Reduction

Krylov Subspace Methods for the Generalized Sylvester Equation.

Dimension Reduction of Large-Scale Second-Order Dynamical Systems Via a Second-Order Arnoldi Method

Stability Analysis of the Two-level Orthogonal Arnoldi Procedure

An Arnoldi-based approach to polynomial and rational least squares problems

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

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

IMPLICITLY RESTARTED GENERALIZED SECOND-ORDER ARNOLDI TYPE ALGORITHMS FOR THE QUADRATIC EIGENVALUE PROBLEM

A sketch-and-select Arnoldi process

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

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

Improving eigenvectors in Arnoldi's method

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

Sketched and Truncated Polynomial Krylov Subspace Methods: Matrix Sylvester Equations