Abstract:A generalized skew-symmetric Lanczos bidiagonalization (GSSLBD) method is proposed to compute several extreme eigenpairs of a large matrix pair $(A,B)$, where $A$ is skew-symmetric and $B$ is symmetric positive definite. The underlying GSSLBD process produces two sets of $B$-orthonormal generalized Lanczos basis vectors that are also $B$-biorthogonal and a series of bidiagonal matrices whose singular values are taken as the approximations to the imaginary parts of the eigenvalues of $(A,B)$ and the corresponding left and right singular vectors premultiplied with the left and right generalized Lanczos basis matrices form the real and imaginary parts of the associated approximate eigenvectors. A rigorous convergence analysis is made on the desired eigenspaces approaching the Krylov subspaces generated by the GSSLBD process and accuracy estimates are made for the approximate eigenpairs. In finite precision arithmetic, it is shown that the semi-$B$-orthogonality and semi-$B$-biorthogonality of the computed left and right generalized Lanczos vectors suffice to compute the eigenvalues accurately. An efficient partial reorthogonalization strategy is adapted to GSSLBD in order to maintain the desired semi-$B$-orthogonality and semi-$B$-biorthogonality. To be practical, an implicitly restarted GSSLBD algorithm, abbreviated as IRGSSLBD, is developed with partial $B$-reorthogonalizations. Numerical experiments illustrate the robustness and overall efficiency of the IRGSSLBD algorithm.

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

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

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

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

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

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

A Contour Integral-Based Algorithm for Computing Generalized Singular Values

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

A generalized skew-symmetric Lanczos bidiagonalization method for computing several extreme eigenpairs of a large skew-symmetric/symmetric positive definite matrix pair

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

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

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

Generalizing the SVD of a matrix under non-standard inner product and its applications to linear ill-posed problems

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

Modifying the Generalized Singular Value Decomposition with Application with Application in Direction-of-arrival Finding

High Order Numerical Methods To Approximate The Singular Value Decomposition

A Power Method for Computing Singular Value Decomposition

Limited Memory Block Krylov Subspace Optimization for Computing Dominant Singular Value Decompositions.

A mixed precision Jacobi SVD algorithm

The Singular Value Decomposition, Applications and Beyond.

Two-sided uniformly randomized GSVD for large-scale discrete ill-posed problems with Tikhonov regularizations