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 Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

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

Triangularized Orthogonalization-Free Method for Solving Extreme Eigenvalue Problems

A new iterative method for simultaneous computation of several eigenpair derivatives of a large matrix

Accelerating Convergence by Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Fast parallelizable methods for the Hermitian eigenvalue problem

A Cubically Convergent Parallelizable Method for the Hermitian Eigenvalue Problem

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

A Ulm-like Method for Inverse Eigenvalue Problems

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

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

A Shifted Method for Solving Symmetric Eigenvalue Problems

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

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

A Sorted Partial Jacobi Method and Its Convergence Analysis

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

Block Algorithms with Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

On the convergence of orthogonalization-free conjugate gradient method for extreme eigenvalues of Hermitian matrices: a Riemannian optimization interpretation