2013IMA Journal of Numerical AnalysisRequires access

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

David Titley-Péloquin, Jennifer Pestana, Andy Wathen

Open publisher page 16 citations

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.

About this research paper

What this paper is about

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.

Why it matters

OpenAlex reports 16 citations for this work. Citation counts describe recorded attention and do not establish research quality.

Key contribution

A contribution statement is not available in the OpenAlex record.

Method / approach

Method details are not available in the OpenAlex metadata.

Main findings

Findings are not separately available in the OpenAlex metadata.

Limitations

Limitations are not available in the OpenAlex metadata.

Applications

Application details are not available in the OpenAlex metadata.

Available 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.

Key concepts: Generalized minimal residual method, Mathematics, Eigenvalues and eigenvectors, Invertible matrix, Diagonalizable matrix, Applied mathematics, Residual, Convergence (economics)

Related papers

Back to paper searchBrowse research topicsOriginal source
GMRES convergence bounds that depend on the right-hand-side vector — Research Paper | ScholarLens