Abstract:Most iterative algorithms for eigenpair computation consist of two main steps: a subspace update (SU) step that generates bases for approximate eigenspaces, followed by a Rayleigh-Ritz (RR) projection step that extracts approximate eigenpairs. So far the predominant methodology for the SU step is based on Krylov subspaces that builds orthonormal bases piece by piece in a sequential manner. In this work, we investigate block methods in the SU step that allow a higher level of concurrency than what is reachable by Krylov subspace methods. To achieve a competitive speed, we propose an augmented Rayleigh-Ritz (ARR) procedure and analyze its rate of convergence under realistic conditions. Combining this ARR procedure with a set of polynomial accelerators, as well as utilizing a few other techniques such as continuation and deflation, we construct a block algorithm designed to reduce the number of RR steps and elevate concurrency in the SU steps. Extensive computational experiments are conducted in Matlab on a representative set of test problems to evaluate the performance of two variants of our algorithm in comparison to two well-established, high-quality eigensolvers ARPACK and FEAST. Numerical results, obtained on a many-core computer without explicit code parallelization, show that when computing a relatively large number of eigenpairs, the performance of our algorithms is competitive with, and frequently superior to, that of the two state-of-the-art eigensolvers.

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

Online Dominant Generalized Eigenvectors Extraction Via a Randomized Algorithm

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

A Two-Level Block Preconditioned Jacobi–Davidson Method for Multiple and Clustered Eigenvalues of Elliptic Operators

Preconditioned iterative methods for eigenvalue counts

Efficient Block Preconditioned Eigensolvers for Linear Response Time-Dependent Density Functional Theory.

Fast parallelizable methods for the Hermitian eigenvalue problem

On the convergence of orthogonalization-free conjugate gradient method for extreme eigenvalues of Hermitian matrices: a Riemannian optimization interpretation

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

Preconditioned Bi-conjugate gradient method of large-scale sparse complex linear equation group

Block Algorithms with Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

A greedy algorithm for computing eigenvalues of a symmetric matrix

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

A Greedy Algorithm for Computing Eigenvalues of a Symmetric Matrix with Localized Eigenvectors

Efficient Preconditioner and Iterative Method for Large Complex Symmetric Linear Algebraic Systems

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

A numerical projection technique for large-scale eigenvalue problems

Parallel Eigenvalue Calculation Based On Multiple Shift-Invert Lanczos And Contour Integral Based Spectral Projection Method