Papers
Topics
Authors
Recent
Gemini 2.5 Flash
Gemini 2.5 Flash
144 tokens/sec
GPT-4o
7 tokens/sec
Gemini 2.5 Pro Pro
46 tokens/sec
o3 Pro
4 tokens/sec
GPT-4.1 Pro
38 tokens/sec
DeepSeek R1 via Azure Pro
28 tokens/sec
2000 character limit reached

Extending Elman's Bound for GMRES (2312.15022v2)

Published 22 Dec 2023 in math.NA and cs.NA

Abstract: If the numerical range of a matrix is contained in the right half of the complex plane, the GMRES algorithm for solving linear systems will reduce the norm of the residual at every iteration. In his Ph.D. dissertation, Howard Elman derived a bound that guarantees convergence. When the numerical range contains the origin, GMRES need not make progress at every step and Elman's bound does not apply, even if all the eigenvalues are located in the right half-plane. However by solving a Lyapunov equation, one can construct an inner product in which the numerical range is contained in the right half-plane. One can then bound GMRES (run in the standard Euclidean norm) by applying Elman's bound in this new inner product, at the cost of a multiplicative constant that characterizes the distortion caused by the change of inner product. Using Lyapunov inverse iteration, one can build a family of such inner products, trading off the location of the numerical range with the size of constant. This approach complements techniques that Greenbaum and colleagues have recently proposed for excising the origin from the numerical range to gain greater insight into the convergence of GMRES for nonnormal matrices.

Citations (2)

Summary

We haven't generated a summary for this paper yet.

X Twitter Logo Streamline Icon: https://streamlinehq.com