Round-Off Error Analysis of a New Class of Conjugate Gradient Algorithms.
Abstract
We perform the rounding error analysis of the conjugate gradient algorithms for the solution of a large system of linear equations Ax = b where A is an hermitian and positive definite matrix. We propose a new class of conjugate gradient algorithms and prove that in the spectral norm the relative error of the computed sequence (x sub k) (in floating point arithmetic) depends at worst on zeta eta to the 3/2 power where zeta is the relative computer precision and eta is the condition number of A. We show that the residual vectors r sub k - Ax sub k-b are at worst of order eta (A) abs. val. x sub k. We point out that with iterative refinement these algorithms are numerically stable. If zeta eta-squared is at most of order unity, then they are also well-behaved. (Author)
Document Details
- Document Type
- Technical Report
- Publication Date
- Dec 01, 1978
- Accession Number
- ADA064293
Entities
People
- H. Wozniakowski
Organizations
- Carnegie Mellon University