Abstract:The FEAST eigensolver is extended to the computation of the singular triplets of a large matrix A with the singular values in a given interval. The resulting FEAST SVDsolver is subspace iteration applied to an approximate spectral projector of A^TA corresponding to the desired singular values in a given interval, and constructs approximate left and right singular subspaces corresponding to the desired singular values, onto which A is projected to obtain Ritz approximations. Differently from a commonly used contour integral-based FEAST solver, we propose a robust alternative that constructs approximate spectral projectors by using the Chebyshev–Jackson polynomial series, which are shown to be symmetric positive semi-definite with the eigenvalues in [0, 1]. We prove the pointwise convergence of this series and give compact estimates for pointwise errors of it and the step function that corresponds to the exact spectral projector of interest. We present error bounds for the approximate spectral projector and reliable estimates for the number of desired singular triplets, prove the convergence of the resulting FEAST SVDsolver, and propose practical selection strategies for determining the series degree and the subspace dimension. The solver and results on it are directly applicable or adaptable to the real symmetric and complex Hermitian eigenvalue problem. Numerical experiments illustrate that the FEAST SVDsolver is substantially more efficient than the contour integral-based FEAST SVDsolver, and it is also more robust and stable than the latter.

A Sequential Subspace Projection Method for Linear Symmetric Eigenvalue Problem.

Fast Calculation of Laplace-Beltrami Eigenproblems Via Subdivision Linear Subspace

A Refined Subspace Iteration Algorithm For Large Sparse Eigenproblems

Gradient-type subspace iteration methods for the symmetric eigenvalue problem

A Nesterov-style Accelerated Gradient Descent Algorithm for the Symmetric Eigenvalue Problem

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

Iterative Methods by Space Decomposition and Subspace Correction

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

A FEAST SVDsolver Based on Chebyshev–Jackson Series for Computing Partial Singular Triplets of Large Matrices

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

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

A Shifted Method for Solving Symmetric Eigenvalue Problems

Theoretical and Computable Optimal Subspace Expansions for Matrix Eigenvalue Problems

An Iterative Method for the Least-Squares Minimum-Norm Symmetric Solution

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

Subspace Methods with Local Refinements for Eigenvalue Computation Using Low-Rank Tensor-Train Format.

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

A Skew-Symmetric Lanczos Bidiagonalization Method for Computing Several Extremal Eigenpairs of a Large Skew-Symmetric Matrix

Fast and Accurate Randomized Algorithms for Linear Systems and Eigenvalue Problems

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

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