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.

On the global convergence of the block Jacobi method for the positive definite generalized eigenvalue problem

Convergence of the complex block Jacobi methods under the generalized serial pivot strategies

On the convergence of complex Jacobi methods

Convergence of the Cyclic and Quasi-cyclic Block Jacobi Methods

Generalized Block Lanczos Methods for Large Unsymmetric Eigenproblems

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Convergence of the alternating iterative method for hermitian positive definite linear systems

On the global convergence of the Jacobi method for symmetric matrices of order 4 under parallel strategies

Convergence of the Eberlein diagonalization method under the generalized serial pivot strategies

Convergence of generalized alternating iterative methods

Convergence of block iterative methods for linear systems with generalized H-matrices

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

Newton-type Methods for Inverse Singular Value Problems with Multiple Singular Values

A Two-Level Block Preconditioned Jacobi–Davidson Method for Multiple and Clustered Eigenvalues of Elliptic Operators

A Variation on the Block Arnoldi Method for Large Unsymmetric Matrix Eigenproblems

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

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

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

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

On parallel multisplitting methods for non-Hermitian positive definite linear systems

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