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.

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

The joint bidiagonalization of a matrix pair with inaccurate inner iterations

The Joint Bidiagonalization Process with Partial Reorthogonalization

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

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

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

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

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

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

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

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

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

A Joint Bidiagonalization Based Algorithm for Large Scale Linear Discrete Ill-posed Problems in General-Form Regularization

A Joint Bidiagonalization Based Iterative Algorithm for Large Scale General-Form Tikhonov Regularization

A mixed precision Jacobi SVD algorithm

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

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

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

A Matrix Polynomial Spectral Approach for General Joint Block Diagonalization.