English

Accurate error estimation in CG

Numerical Analysis 2021-01-12 v1 Numerical Analysis

Abstract

In practical computations, the (preconditioned) conjugate gradient (P)CG method is the iterative method of choice for solving systems of linear algebraic equations Ax=bAx=b with a real symmetric positive definite matrix AA. During the iterations it is important to monitor the quality of the approximate solution xkx_k so that the process could be stopped whenever xkx_k is accurate enough. One of the most relevant quantities for monitoring the quality of xkx_k is the squared AA-norm of the error vector xxkx-x_k. This quantity cannot be easily evaluated, however, it can be estimated. Many of the existing estimation techniques are inspired by the view of CG as a procedure for approximating a certain Riemann--Stieltjes integral. The most natural technique is based on the Gauss quadrature approximation and provides a lower bound on the quantity of interest. The bound can be cheaply evaluated using terms that have to be computed anyway in the forthcoming CG iterations. If the squared AA-norm of the error vector decreases rapidly, then the lower bound represents a tight estimate. In this paper we suggest a heuristic strategy aiming to answer the question of how many forthcoming CG iterations are needed to get an estimate with the prescribed accuracy. Numerical experiments demonstrate that the suggested strategy is efficient and robust.

Keywords

Cite

@article{arxiv.2101.03931,
  title  = {Accurate error estimation in CG},
  author = {Gérard Meurant and Jan Papež and Petr Tichý},
  journal= {arXiv preprint arXiv:2101.03931},
  year   = {2021}
}

Comments

21 pages, 8 figures

R2 v1 2026-06-23T21:59:39.805Z