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.

A Class of Inverse Eigenvalue Problems for Real Symmetric Banded Matrices with Odd Bandwidth

Solution of the symmetric band partial inverse eigenvalue problem for the damped mass spring system

Inverse eigenvalue problem for pseudo-symmetric Jacobi matrices with two spectra

Determinants and Inverses of banded Toeplitz Matrices over $\mathbb{F}_p$ Are Periodic

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

Solutions to a Quadratic Inverse Eigenvalue Problem

Isospectral Flow Method for Nonnegative Inverse Eigenvalue Problem with Prescribed Structure

An inversion algorithm for a banded matrix

Efficient Computation for Invertibility Sequence of Banded Toeplitz Matrices

Test for Bandedness of High-Dimensional Covariance Matrices and Bandwidth Estimation

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

Inverse Eigenvalue Problem for Bordered Diagonal Matrices

On the Inverse Symmetric Quadratic Eigenvalue Problem

On the discrete spectrum of complex banded matrices

A Communication-Avoiding Parallel Algorithm for the Symmetric Eigenvalue Problem

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

Banded matrices and their orthogonality

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

The Non-negative Inverse Eigenvalue Problem for Tridiagonal Matrix, Circulant Matrix and Symmetric Matrix

Application of the Schur-Cohn Theorem to the precise convergence domain for a p-cyclic SOR iteration matrix