205 research outputs found
Some observations on weighted GMRES
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 new alternative implementations of the weighted Arnoldi algorithm which may be favorable in terms of computational complexity, and examine stability issues connected with these implementations. Two 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
A flexible and adaptive Simpler GMRES with deflated restarting for shifted linear systems
In this paper, two efficient iterative algorithms based on the simpler GMRES
method are proposed for solving shifted linear systems. To make full use of the
shifted structure, the proposed algorithms utilizing the deflated restarting
strategy and flexible preconditioning can significantly reduce the number of
matrix-vector products and the elapsed CPU time. Numerical experiments are
reported to illustrate the performance and effectiveness of the proposed
algorithms.Comment: 17 pages. 9 Tables, 1 figure; Newly update: add some new numerical
results and correct some typos and syntax error
Deflated Iterative Methods for Linear Equations with Multiple Right-Hand Sides
A new approach is discussed for solving large nonsymmetric systems of linear
equations with multiple right-hand sides. The first system is solved with a
deflated GMRES method that generates eigenvector information at the same time
that the linear equations are solved. Subsequent systems are solved by
combining restarted GMRES with a projection over the previously determined
eigenvectors. This approach offers an alternative to block methods, and it can
also be combined with a block method. It is useful when there are a limited
number of small eigenvalues that slow the convergence. An example is given
showing significant improvement for a problem from quantum chromodynamics. The
second and subsequent right-hand sides are solved much quicker than without the
deflation. This new approach is relatively simple to implement and is very
efficient compared to other deflation methods.Comment: 13 pages, 5 figure
- …