PROPERTIES OF WORST-CASE GMRES

被引:6
|
作者
Faber, Vance [1 ]
Liesen, Joerg [2 ]
Tichy, Petr [3 ]
机构
[1] Vanco Res, Big Pine Key, FL 33043 USA
[2] Tech Univ Berlin, Inst Math, D-10623 Berlin, Germany
[3] Acad Sci Czech Republic, Inst Comp Sci, Prague 18207, Czech Republic
关键词
GMRES method; worst-case convergence; ideal GMRES; matrix approximation; problems; minmax; NONSYMMETRIC LINEAR-SYSTEMS; IDEAL GMRES; MATRIX; POLYNOMIALS; ALGORITHM;
D O I
10.1137/13091066X
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In the convergence analysis of the GMRES method for a given matrix A, one quantity of interest is the largest possible residual norm that can be attained, at a given iteration step k, over all unit norm initial vectors. This quantity is called the worst-case GMRES residual norm for A and k. We show that the worst-case behavior of GMRES for the matrices A and A(T) is the same, and we analyze properties of initial vectors for which the worst-case residual norm is attained. In particular, we prove that such vectors satisfy a certain "cross equality." 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 strict at all iteration steps k >= 3. Finally, we give a complete characterization of how the values of the approximation problems change in the context of worst-case and ideal GMRES for a real matrix, when one considers complex (rather than real) polynomials and initial vectors.
引用
收藏
页码:1500 / 1519
页数:20
相关论文
共 50 条