Article thumbnail
Location of Repository

How descriptive are GMRES convergence bounds?

By Mark Embree

Abstract

Eigenvalues with the eigenvector condition number, the field of values, and pseudospectra have all been suggested as the basis for convergence bounds for minimum residual Krylov subspace methods applied to non-normal coefficient matrices. This paper analyzes and compares these bounds, illustrating with six examples the success and failure of each one. Refined bounds based on eigenvalues and the field of values are suggested to handle low-dimensional non-normality. It is observed that pseudospectral bounds can capture multiple convergence stages. Unfortunately, computation of pseudospectra can be rather expensive. This motivates an adaptive technique for estimating GMRES convergence based on approximate pseudospectra taken from the Arnoldi process that is the basis for GMRES

Topics: Linear and multilinear algebra; matrix theory, Potential theory, Numerical analysis
Publisher: Unspecified
Year: 1999
OAI identifier: oai:generic.eprints.org:1290/core69

Suggested articles


To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.