Abstract
We investigate the convergence of the weighted GMRES method for solving linear systems. Two different weighting variants are compared with unweighted GMRES for three model problems, giving a phenomenological explanation of cases where weighting improves convergence, and a case where weighting has no effect on the convergence. We also present a new alternative implementation of the weighted Arnoldi algorithm which under known circumstances will be favourable in terms of computational complexity. These implementations of weighted GMRES are compared for a large number of examples. We find that weighted GMRES may outperform unweighted GMRES for some problems, but more often this method is not competitive with other Krylov subspace methods like GMRES with deflated restarting or BICGSTAB, in particular when a preconditioner is used.
Original language | English |
---|---|
Pages (from-to) | 733-752 |
Number of pages | 21 |
Journal | Numerical Algorithms |
Volume | 67 |
Issue number | 4 |
Early online date | 10 Jan 2014 |
DOIs | |
Publication status | Published - Dec 2014 |
Keywords
- weighted GMRES
- linear systems
- Krylov subspace method
- harmonic Ritz values
- algorithms
- algebra