Abstract
The convergence of GMRES for solving linear systems can be influenced heavily by the structure of the right hand side. Within the solution of eigenvalue problems via inverse iteration or subspace iteration, the right hand side is generally related to an approximate invariant subspace of the linear system. We give detailed and new bounds on (block) GMRES that take the special behavior of the right hand side into account and explain the initial sharp decrease of the GMRES residual. The bounds motivate the use of specific preconditioners for these eigenvalue problems, e.g. tuned and polynomial preconditioners, as we describe. The numerical results show that the new (block) GMRES bounds are much sharper than conventional bounds and that preconditioned subspace iteration with either a tuned or polynomial preconditioner should be used in practice.
Original language | English |
---|---|
Pages (from-to) | 203-222 |
Number of pages | 20 |
Journal | Computational Methods in Applied Mathematics |
Volume | 18 |
Issue number | 2 |
Early online date | 7 Jun 2017 |
DOIs | |
Publication status | Published - 30 Apr 2018 |
Keywords
- GMRES
- convergence analysis
- inexact inverse iteration
- inexact subspace iteration
- Krylov subspace methods
- block Krylov methods
- preconditioning