On Error Estimation in the Conjugate Gradient Method: Normwise Backward Error
Abstrakt
Using an idea of Duff and Vömel [BIT, 42 (2002), pp. 300-322 ] we suggest a simple algorithm that incrementally estimates the 2-norm of Jacobi matrices that are available during the conjugate gradient (CG) computations. The estimate can be used, e.g., in stopping criteria based on the normwise backward error. Numerical experiments predict that the estimate approximates the 2-norm of A with a suffcient accuracy.