A Lanczos-type method for multiple starting vectors
Open Access
- 20 May 1999
- journal article
- Published by American Mathematical Society (AMS) in Mathematics of Computation
- Vol. 69 (232) , 1577-1602
- https://doi.org/10.1090/s0025-5718-99-01163-1
Abstract
Given a square matrix and single right and left starting vectors, the classical nonsymmetric Lanczos process generates two sequences of biorthogonal basis vectors for the right and left Krylov subspaces induced by the given matrix and vectors. In this paper, we propose a Lanczos-type algorithm that extends the classical Lanczos process for single starting vectors to multiple starting vectors. Given a square matrix and two blocks of right and left starting vectors, the algorithm generates two sequences of biorthogonal basis vectors for the right and left block Krylov subspaces induced by the given data. The algorithm can handle the most general case of right and left starting blocks of arbitrary sizes, while all previously proposed extensions of the Lanczos process are restricted to right and left starting blocks of identical sizes. Other features of our algorithm include a built-in deflation procedure to detect and delete linearly dependent vectors in the block Krylov sequences, and the option to employ look-ahead to remedy the potential breakdowns that may occur in nonsymmetric Lanczos-type methods.Keywords
This publication has 25 references indexed in Scilit:
- Iterative solution of multiple radiation and scattering problems in structural acoustics using a block quasi-minimal residual algorithmComputer Methods in Applied Mechanics and Engineering, 1997
- QMR: a quasi-minimal residual method for non-Hermitian linear systemsNumerische Mathematik, 1991
- The nonsymmetric Lanczos algorithm and controllabilitySystems & Control Letters, 1991
- A Look-Ahead Lanczos Algorithm for Unsymmetric MatricesMathematics of Computation, 1985
- The determination of structural properties of a linear multivariable system by operations of system similarity 2. Non-proper systems in generalized state-space form†International Journal of Control, 1980
- Recursive Algorithms for the Matrix Pade ProblemMathematics of Computation, 1980
- Input-output invariants for linear multivariable systemsIEEE Transactions on Automatic Control, 1980
- Implementation Aspects of Band Lanczos Algorithms for Computation of Eigenvalues of Large Sparse Symmetric MatricesMathematics of Computation, 1979
- Matrix interpretations and applications of the continued fraction algorithmRocky Mountain Journal of Mathematics, 1974
- Algorithms for Triangular Decomposition of Block Hankel and Toeplitz Matrices with Application to Factoring Positive Matrix PolynomialsMathematics of Computation, 1973