Abstract:The harmonic Lanczos bidiagonalization method can be used to compute the smallest singular triplets of a large matrix $A$. We prove that for good enough projection subspaces harmonic Ritz values converge if the columns of $A$ are strongly linearly independent. On the other hand, harmonic Ritz values may miss some desired singular values when the columns of $A$ are almost linearly dependent. Furthermore, harmonic Ritz vectors may converge irregularly and may even fail to converge. Based on the refined projection principle for large matrix eigenproblems due to the first author, we propose a refined harmonic Lanczos bidiagonalization method that takes the Rayleigh quotients of the harmonic Ritz vectors as approximate singular values and extracts the best approximate singular vectors, called the refined harmonic Ritz approximations, from the given subspaces in the sense of residual minimizations. The refined approximations are shown to converge to the desired singular vectors once the subspaces are sufficiently good and the Rayleigh quotients converge. An implicitly restarted refined harmonic Lanczos bidiagonalization algorithm (IRRHLB) is developed. We study how to select the best possible shifts, and suggest refined harmonic shifts that are theoretically better than the harmonic shifts used within the implicitly restarted harmonic Lanczos bidiagonalization algorithm (IRHLB). We propose a novel procedure that can numerically compute the refined harmonic shifts efficiently and accurately. Numerical experiments are reported that compare IRRHLB with five other algorithms based on the Lanczos bidiagonalization process. It appears that IRRHLB is at least competitive with them and can be considerably more efficient when computing the smallest singular triplets.

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

Refined and Refined Harmonic Jacobi–Davidson Methods for Computing Several GSVD Components of a Large Regular Matrix Pair

A Cross-Product Free Jacobi–Davidson Type Method 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

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

Preconditioning correction equations in Jacobi--Davidson type methods for computing partial singular value decompositions of large matrices

On Choices of Formulations of Computing the Generalized Singular Value Decomposition of a Large Matrix Pair.

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

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

Harmonic and Refined Harmonic Shift-Invert Residual Arnoldi and Jacobi-Davidson Methods for Interior Eigenvalue Problems

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

Characterizing GSVD by singular value expansion of linear operators and its computation

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

The Joint Bidiagonalization of a Matrix Pair with Inaccurate Inner Iterations

A Higher-Order Generalized Singular Value Decomposition for Rank Deficient Matrices

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

A Matrix Polynomial Spectral Approach for General Joint Block Diagonalization.

The Joint Bidiagonalization Process with Partial Reorthogonalization

A Refined Harmonic Lanczos Bidiagonalization Method and an Implicitly Restarted Algorithm for Computing the Smallest Singular Triplets of Large Matrices

A REFINED HARMONIC LANCZOS DIDIAGONALIZATION FOR COMPUTING SMALLEST SINGULAR TRIPLETS