Abstract:Iterative algorithms for large-scale eigenpair computation of symmetric matrices are mostly based on subspace projections consisting of two main steps: a subspace update (SU) step that generates bases for approximate eigenspaces, followed by a Rayleigh-Ritz projection step that extracts approximate eigenpairs. A predominant methodology for the SU step makes use of Krylov subspaces and builds orthonormal bases piece by piece in a sequential manner. On the other hand, block methods such as the classic (simultaneous) subspace iteration, allow higher levels of concurrency than what is reachable by Krylov subspace methods, but may suffer from slow convergence. In this work, we analyze the rate of convergence for a simple block algorithmic framework that combines an augmented Rayleigh-Ritz (ARR) procedure with the subspace iteration. Our main results are Theorem 4.5 and its corollaries, which show that the ARR procedure can provide significant accelerations to convergence speed. Our analysis will offer useful guidelines for designing and implementing practical algorithms from this framework.

A Cubically Convergent Parallelizable Method for the Hermitian Eigenvalue Problem

Fast parallelizable methods for the Hermitian eigenvalue problem

Fast Parallelizable Methods for Computing Invariant Subspaces of Hermitian Matrices

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

A Quadratically Convergent QR-like Method Without Shifts for the Hermitian Eigenvalue Problem

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

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

Accelerating Parallel Jacobi Method for Matrix Eigenvalue Computation in DOA Estimation Algorithm

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

A Ulm-like Method for Inverse Eigenvalue Problems

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

Generalized Jacobi Method for Computing Eigenvalues of Dual Quaternion Hermitian Matrices

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

Convergence Analysis of Newton-Like Methods for Inverse Eigenvalue Problems with Multiple Eigenvalues

Accelerating Convergence by Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

Approximate Cayley Transform Methods for Inverse Eigenvalue Problems and Convergence Analysis

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Numerical Analysis of the Parallel Orbital-Updating Approach for Eigenvalue Problems

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems