An augmented matrix-based CJ-FEAST SVDsolver for computing a partial singular value decomposition with the singular values in a given interval

by   Zhongxiao Jia, et al.

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 considered to compute the singular triplets of a large matrix A with the singular values in an interval [a,b] contained in the singular spectrum. The new CJ-FEAST SVDsolver is a subspace iteration applied to an approximate spectral projector of the augmented matrix [0, A^T; A, 0] associated with the eigenvalues in [a,b], and constructs approximate left and right singular subspaces with the desired singular values independently, onto which A is projected to obtain the Ritz approximations to the desired singular triplets. Compact estimates are given for the accuracy of the approximate spectral projector, 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 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.


page 1

page 2

page 3

page 4


A FEAST SVDsolver for the computation of singular value decompositions of large matrices based on the Chebyshev–Jackson series expansion

The FEAST eigensolver is extended to the computation of the singular tri...

A skew-symmetric Lanczos bidiagonalization method for computing several largest eigenpairs of a large skew-symmetric matrix

The spectral decomposition of a real skew-symmetric matrix A can be math...

Cross product-free matrix pencils for computing generalized singular values

It is well known that the generalized (or quotient) singular values of a...

Convergence directions of the randomized Gauss–Seidel method and its extension

The randomized Gauss–Seidel method and its extension have attracted much...

Towards a more robust algorithm for computing the restricted singular value decomposition

A new algorithm to compute the restricted singular value decomposition o...

Fast Interpolation-based Globality Certificates for Computing Kreiss Constants and the Distance to Uncontrollability

The Kreiss constant of a matrix and the distance to uncontrollability ca...

Please sign up or login with your details

Forgot password? Click here to reset