GMRES convergence bounds that depend on the right-hand-side vector
David Titley-Péloquin, Jennifer Pestana, Andy Wathen
Abstract
David Titley-Péloquin, Jennifer Pestana, Andy Wathen
Abstract
We consider the convergence of the algorithm GMRES of Saad and Schultz for solving linear equations Bx=b, where B ∈ ℂn × n is nonsingular and diagonalizable, and b ∈ ℂn. Our analysis explicitly includes the initial residual vector r0. 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.
OpenAlex reports 16 citations for this work. Citation counts describe recorded attention and do not establish research quality.
A contribution statement is not available in the OpenAlex record.
Method details are not available in the OpenAlex metadata.
Findings are not separately available in the OpenAlex metadata.
Limitations are not available in the OpenAlex metadata.
Application details are not available in the OpenAlex metadata.
We consider the convergence of the algorithm GMRES of Saad and Schultz for solving linear equations Bx=b, where B ∈ ℂn × n is nonsingular and diagonalizable, and b ∈ ℂn. Our analysis explicitly includes the initial residual vector r0. 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.
Key concepts: Generalized minimal residual method, Mathematics, Eigenvalues and eigenvectors, Invertible matrix, Diagonalizable matrix, Applied mathematics, Residual, Convergence (economics)