English

Error estimators and their analysis for CG, Bi-CG and GMRES

Numerical Analysis 2022-09-12 v3 Numerical Analysis

Abstract

The demands of accuracy in measurements and engineering models today, renders the condition number of problems larger. While a corresponding increase in the precision of floating point numbers ensured a stable computing, the uncertainty in convergence when using residue as a stopping criterion has increased. We present an analysis of the uncertainty in convergence when using relative residue as a stopping criterion for iterative solution of linear systems, and the resulting over/under computation for a given tolerance in error. This shows that error estimation is significant for an efficient or accurate solution even when the condition number of the matrix is not large. An O(1)\mathcal{O}(1) error estimator for iterations of the CG algorithm was proposed more than two decades ago. Recently, an O(k2)\mathcal{O}(k^2) error estimator was described for the GMRES algorithm which allows for non-symmetric linear systems as well, where kk is the iteration number. We suggest a minor modification in this GMRES error estimation for increased stability. In this work, we also propose an O(n)\mathcal{O}(n) error estimator for A-norm and l2l_{2} norm of the error vector in Bi-CG algorithm. The robust performance of these estimates as a stopping criterion results in increased savings and accuracy in computation, as condition number and size of problems increase.

Keywords

Cite

@article{arxiv.1705.08806,
  title  = {Error estimators and their analysis for CG, Bi-CG and GMRES},
  author = {Puneet Jain and Krishna Manglani and Murugesan Venkatapathi},
  journal= {arXiv preprint arXiv:1705.08806},
  year   = {2022}
}