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.

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

Second-order Krylov Subspace and Arnoldi Procedure

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

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

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

Stability Analysis of the Two-level Orthogonal Arnoldi Procedure

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

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

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

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

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

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

A Ulm-like Method for Inverse Eigenvalue Problems

An Arnoldi method for unsymmetric eigenproblems and its parallel implementation

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

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

Krylov Subspace Methods for the Generalized Sylvester Equation.

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

ARNOLDI TYPE ALGORITHMS FOR LARGE UNSYMMETRIC MULTIPLE EIGENVALUE PROBLEMS

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