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 Second-Order Biorthogonalization Procedure and Its Application to Quadratic Eigenvalue Problems

Second-order Krylov Subspace and Arnoldi Procedure

Restarted Generalized Second-Order Krylov Subspace Methods for Solving Quadratic Eigenvalue Problems

SOAR: A Second-order Arnoldi Method for the Solution of the Quadratic Eigenvalue Problem.

A Block Incomplete Orthogonalization Method for Large Nonsymmetric Eigenproblems

Triangularized Orthogonalization-Free Method for Solving Extreme Eigenvalue Problems

A Refined Second-order Arnoldi (RSOAR) Method for the Quadratic Eigenvalue Problem and Implicit Restarting

IMPLICITLY RESTARTED GENERALIZED SECOND-ORDER ARNOLDI TYPE ALGORITHMS FOR THE QUADRATIC EIGENVALUE PROBLEM

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

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

GPBiLQ and GPQMR: Two iterative methods for unsymmetric partitioned linear systems

Restarted Generalized Krylov Subspace Methods for Solving Large-Scale Polynomial Eigenvalue Problems

On the Inverse Symmetric Quadratic Eigenvalue Problem

Stability Analysis of the Two-level Orthogonal Arnoldi Procedure

A Refined Jacobi-Davidson Method for the Quadratic Eigenvalue Problem

Quadratic transformation and matrix biorthogonal polynomials: an $\mathcal{LU}$ factorization approach

Solutions of the Partially Described Inverse Quadratic Eigenvalue Problem

The eigenvalue decomposition of normal matrices by the decomposition of the skew-symmetric part with applications to orthogonal matrices

The Joint Bidiagonalization Process with Partial Reorthogonalization

Krylov subspace methods for large multidimensional eigenvalue computation

Spectral Decomposition of Real Symmetric Quadratic Λ-Matrices and Its Applications.