For solving the large scale quadratic eigenvalue problem L(λ)x: = (Aλ2 + Bλ + C)x = 0, a direct projection method based on the Krylov subspaces generated by a single matrix A-1B using the standard Arnoldi algorithm is considered. It is shown that, when iteratively combined with the shift-and-invert technique, it results in a fast converging algorithm. The important situations of inexact shift-and-invert are also discussed and numerical examples are presented to illustrate the new method.
|Number of pages||10|
|Journal||Applied Mathematics and Computation|
|Issue number||2 SPEC. ISS.|
|State||Published - Jan 15 2006|
ASJC Scopus subject areas
- Computational Mathematics
- Applied Mathematics