Inexact preconditioned conjugate gradient method with inner-outer iteration

Gene H. Golub, Qiang Ye

Research output: Contribution to journalArticlepeer-review

156 Scopus citations

Abstract

An important variation of preconditioned conjugate gradient algorithms is inexact preconditioner implemented with inner-outer iterations [G. H. Golub and M.L. Overton, Numerical Analysis, Lecture Notes in Math. 912, Springer, Berlin, New York, 1982], where the preconditioner is solved by an inner iteration to a prescribed precision. In this paper, we formulate an inexact preconditioned conjugate gradient algorithm for a symmetric positive definite system and analyze its convergence property. We establish a linear convergence result using a local relation of residual norms. We also analyze the algorithm using a global equation and show that the algorithm may have the superlinear convergence property when the inner iteration is solved to high accuracy. The analysis is in agreement with observed numerical behavior of the algorithm. In particular, it suggests a heuristic choice of the stopping threshold for the inner iteration. Numerical examples are given to show the effectiveness of this choice and to compare the convergence bound.

Original languageEnglish
Pages (from-to)1305-1320
Number of pages16
JournalUnknown Journal
Volume21
Issue number4
DOIs
StatePublished - Dec 1999

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Inexact preconditioned conjugate gradient method with inner-outer iteration'. Together they form a unique fingerprint.

Cite this