GMRES convergence bounds that depend on the right-hand-side vector

摘要

We consider the convergence of the algorithm GMRES of Saad and Schultz for solving linear equations Bx=b, where B ∈ C~(n×n) is nonsingular and diagonalizable, and b ∈ C~n. Our analysis explicitly includes the initial residual vector r_0. We show that the GMRES residual norm satisfies a weighted polynomial least-squares problem on the spectrum of B, and that GMRES convergence reduces to an ideal GMRES problem on a rank-1 modification of the diagonal matrix of eigenvalues of B. Numerical experiments show that the new bounds can accurately describe GMRES convergence.