Relating computed and exact entities in methods based on Lanczos tridiagonalization

• ID: 2792, RIV: 10384797
• ISSN: not specified, ISBN: 978-3-319-97135-3
• source: Lecture Notes in Computer Science
• keywords: Krylov subspace; CGL; MINRES; finite-precision computations; loss of orthogonality; delay of convergence; Lanczos vectors
• authors: Tomáš Gergelits, Iveta Hnětynková, Marie Kubínová
• authors from KNM: Hnětynková Iveta

Abstract

Krylov subspace methods based on short recurrences such as CGL or MINRES represent an attractive way of solving large and sparse systems of linear algebraic equations. Loss of orthogonality in the underlying Lanczos process delays significantly their convergence in finite-precision computation, whose connection to exact computation is still not fully understood. In this paper, we exploit the idea of simultaneous comparison of finite-precision and exact computations for CGL and MINRES, by taking advantage of their relationship valid also in finite-precision arithmetic. In particular, we show that finite-precision CGL residuals and Lanczos vectors have to be aggregated over the intermediate iterations to form a counterpart to vectors from the exact computation. Influence of stagnation in exact MINRES computation is also discussed. Obtained results are supported by numerical experiments.