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.

Generalized Golub-Kahan bidiagonalization for nonsymmetric saddle point systems

Preconditioned iterative method for nonsymmetric saddle point linear systems

On Semi-Convergence and Inexact Iteration of the GSS Iteration Method for Nonsymmetric Singular Saddle Point Problems

Deflation for the off-diagonal block in symmetric saddle point systems

A null-space approach for large-scale symmetric saddle point systems with a small and non zero (2, 2) block

On the invertibility of matrices with a double saddle-point structure

The Alternate Direction Iterative Methods for Generalized Saddle Point Systems

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

A generalized relaxed block positive-semidefinite splitting preconditioner for generalized saddle point linear system

A Block Positive-Semidefinite Splitting Preconditioner for Generalized Saddle Point Linear Systems.

A general alternating direction implicit iteration method for solving complex symmetric linear systems

A Class of Preconditioners Based on Symmetric-Triangular Decomposition and Matrix Splitting for Generalized Saddle Point Problems

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

A Note on the Block Alternating Splitting Implicit Iteration Method for Complex Saddle‐point Problems

Ideal Preconditioners for Saddle Point Systems with a Rank-Deficient Leading Block

An Enhancement of the Accuracy of the BiCGStab Method for Solving Linear Systems with Single or Multiple Right-Hand Sides

Regularized DPSS Preconditioners for Generalized Saddle Point Linear Systems

A Condition for the Nonsymmetric Saddle Point Matrix Being Diagonalizable and Having Real and Positive Eigenvalues

A Note on Constraint Preconditioners for Nonsymmetric Saddle Point Problems

The Joint Bidiagonalization of a Matrix Pair with Inaccurate Inner Iterations