Abstract
In this note, we examine a modified conjugate gradient procedure for solving \(A\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{x} = \underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\) in which the approximation space is based upon the Krylov space (\(A\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{x} = \underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\)) associated with \(\sqrt A\) and \(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\). We show that, given initial vectors \(\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\) and \(\sqrt A \underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\) (possibly computed at some expense), the best fit solution in \(K^k \sqrt A ,\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{b}\) can be computed using a finite-term recurrence requiring only one multiplication by A per iteration. The initial convergence rate appears, as expected, to be twice as fast as that of the standard conjugate gradient method, but stability problems cause the convergence to be degraded.
Similar content being viewed by others
REFERENCES
Birkhoff, G., and Lynch, R. E. (1984) Numerical Solution of Elliptic Problems, Siam, Philadelphia.
Golub, G. H., and Van Loan, C. H. (1989). Matrix Computations, Johns Hopkins University Press, Baltimore, Maryland.
Lanczos, C. (1952) Solution of Systems of linear equations by minimized iterations, J. Res. Nat. Bureau of Standards, 49, 33-53.
O'Leary, D. P. (1980). The block conjugate gradient algorithm and related methods, Linear Algebra and Its Appl. 29, 293-322.
Simon, H. D. (1984). The Lanczos method with partial reorthogonalization, Math. Comp. 42, 115-142.
van der Vorst, H. A. (1987). An iterative solution method for solving f(A)x = b, using Krylov subspace information obtained for the symmetric positive definite matrix A, J. Comp. Appl. Math. 18, 249-263.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Gottlieb, S., Fischer, P.F. Modified Conjugate Gradient Method for the Solution of Ax=b. Journal of Scientific Computing 13, 173–183 (1998). https://doi.org/10.1023/A:1023222110984
Issue Date:
DOI: https://doi.org/10.1023/A:1023222110984