Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 45, Issue 1, Page 24-58, March 2024. The cross-product matrix-based CJ-FEAST SVDsolver proposed previously by the authors is shown to compute the left singular vector possibly much less accurately than the right singular vector and may be numerically backward unstable when a desired singular value is small. In this paper, an alternative augmented matrix-based CJ-FEAST SVDsolver is proposed to compute the singular triplets of a large matrix [math] with the singular values in an interval [math] contained in the singular spectrum. The new CJ-FEAST SVDsolver is a subspace iteration applied to an approximate spectral projector of the augmented matrix [math] associated with the eigenvalues in [math], and it constructs approximate left and right singular subspaces independently, onto which [math] is projected to obtain the Ritz approximations to the desired singular triplets. Compact estimates are given for the accuracy of the approximate spectral projector constructed by the Chebyshev–Jackson series expansion in terms of series degree, and a number of convergence results are established. The new solver is proved to be always numerically backward stable. A convergence comparison of the cross-product-based and augmented matrix-based CJ-FEAST SVDsolvers is made, and a general-purpose choice strategy between the two solvers is proposed for the robustness and overall efficiency. Numerical experiments confirm all the results and meanwhile demonstrate that the proposed solver is more robust and substantially more efficient than the corresponding contour integral-based versions that exploit the trapezoidal rule and the Gauss–Legendre quadrature to construct an approximate spectral projector.

An Augmented Matrix-Based CJ-FEAST SVDsolver for Computing a Partial Singular Value Decomposition with the Singular Values in a Given Interval

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

A CJ-FEAST GSVDsolver for computing a partial GSVD of a large matrix pair with the generalized singular values in a given interval

Using Cross-Product Matrices to Compute the SVD.

Extended Newton-type Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

An Extended Ulm-like Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

A Novel Fully Hardware-Implemented SVD Solver Based on Ultra-Parallel BCV Jacobi Algorithm

Accelerated singular value thresholding for matrix completion.

A Cross-Product Free Jacobi–Davidson Type Method for Computing a Partial Generalized Singular Value Decomposition of a Large Matrix Pair

A Contour Integral-Based Algorithm for Computing Generalized Singular Values

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

The Joint Bidiagonalization Method for Large GSVD Computations in Finite Precision.

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

Newton-type Methods for Inverse Singular Value Problems with Multiple Singular Values

High Order Numerical Methods To Approximate The Singular Value Decomposition

A mixed precision Jacobi SVD algorithm

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

PRIMME_SVDS: A High-Performance Preconditioned SVD Solver for Accurate Large-Scale Computations

Two Harmonic Jacobi–Davidson Methods for Computing a Partial Generalized Singular Value Decomposition of a Large Matrix Pair

On Inner Iterations of Jacobi--Davidson Type Methods for Large SVD Computations

An Efficient GIPM Algorithm for Computing the Smallest V-Singular Value of the Partially Symmetric Tensor