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.

A New Shift of the QL Algorithm for Irreducible Symmetric Tridiagonal Matrices

Global Convergence of Hessenberg Shifted QR I: Exact Arithmetic

An improved Shifted CholeskyQR based on columns

A Quadratically Convergent QR-like Method Without Shifts for the Hermitian Eigenvalue Problem

A Shifted Method for Solving Symmetric Eigenvalue Problems

Algorithm 1019: A Task-based Multi-shift QR/QZ Algorithm with Aggressive Early Deflation

A new real structure-preserving quaternion QR algorithm

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

QMRCGstab Algorithm for Families of Shifted Linear Systems

CholeskyQR for sparse matrices

A REFINED HARMONIC LANCZOS DIDIAGONALIZATION FOR COMPUTING SMALLEST SINGULAR TRIPLETS

ON THE PURE IMAGINARY QUATERNIONIC LEAST SQUARES SOLUTIONS OF MATRIX EQUATION

A SQMRCGstab Algorithm for Families of Shifted Linear Systems

Real dqds for the nonsymmetric tridiagonal eigenvalue problem

On Krylov subspace methods for skew-symmetric and shifted skew-symmetric linear systems

A Class of Inverse Problem with Symmetric Skew Anti-symmetric Tridiagonal Matrices

A Refined Shift-And-Invert Arnoldi Algorithm for Large Unsymmetric Generalized Eigenproblems

On rank‐revealing QR factorizations of quaternion matrices

A New Implementation of LSMR Algorithm for the Quaternionic Least Squares Problem

Algorithm 953

TR-2008007: Degeneration of Structured Integer Matrices Modulo an Integer