Characterization of worst-case GMRES
Abstract
Given a matrix and iteration step , we study a best possible attainable upper bound on the GMRES residual norm that does not depend on the initial vector . This quantity is called the worst-case GMRES approximation. We show that the worst case behavior of GMRES for the matrices and is the same, and we analyze properties of initial vectors for which the worst-case residual norm is attained. In particular, we show that such vectors satisfy a certain "cross equality", and we characterize them as right singular vectors of the corresponding GMRES residual matrix. We show that the worst-case GMRES polynomial may not be uniquely determined, and we consider the relation between the worst-case and the ideal GMRES approximations, giving new examples in which the inequality between the two quantities is sharp at all iteration steps . Finally, we give a complete characterization of how the values of the approximation problems in the context of worst-case and ideal GMRES for a real matrix change, when one considers complex (rather than real) polynomials and initial vectors in these problems.
Cite
@article{arxiv.1302.5535,
title = {Characterization of worst-case GMRES},
author = {Vance Faber and Jörg Liesen and Petr Tichý},
journal= {arXiv preprint arXiv:1302.5535},
year = {2013}
}