Relating computed and exact entities in methods based on Lanczos tridiagonalization

T Gergelits, I Hnětynková, M Kubínová - High Performance Computing in …, 2018 - Springer
T Gergelits, I Hnětynková, M Kubínová
High Performance Computing in Science and Engineering: Third International …, 2018Springer
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 …
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.
Springer
Showing the best result for this search. See all results