Volume 46, pp. 162-189, 2017.

Stagnation of block GMRES and its relationship to block FOM

Kirk M. Soodhalter

Abstract

We analyze the convergence behavior of block GMRES and characterize the phenomenon of stagnation which is then related to the behavior of the block FOM method. We generalize the block FOM method to generate well-defined approximations in the case that block FOM would normally break down, and these generalized solutions are used in our analysis. This behavior is also related to the principal angles between the column-space of the previous block GMRES residual and the current minimum residual constraint space. At iteration j, it is shown that the proper generalization of GMRES stagnation to the block setting relates to the column space of the jth block Arnoldi vector. Our analysis covers both the cases of normal iterations as well as block Arnoldi breakdown wherein dependent basis vectors are replaced with random ones. Numerical examples are given to illustrate what we have proven, including one built from a small application problem to demonstrate the validity of the analysis in a less pathological case.

Full Text (PDF) [516 KB], BibTeX

Key words

block Krylov subspace methods, GMRES, FOM, stagnation

AMS subject classifications

65F10, 65F50, 65F08

Links to the cited ETNA articles

[12] Vol. 33 (2008-2009), pp. 207-220 L. Elbouyahyaoui, A. Messaoudi, and H. Sadok: Algebraic properties of the block GMRES and block Arnoldi methods
[25] Vol. 39 (2012), pp. 75-101 Gérard Meurant: The complete stagnation of GMRES for n≤4