Paper detail

Non-stationary extremal eigenvalue approximations in iterative solutions of linear systems and estimators for relative error

Non-stationary approximations of the final value of a converging sequence are discussed, and we show that extremal eigenvalues can be reasonably estimated from the CG iterates without much computation at all. We introduce estimators of relative error for conjugate gradient (CG)methods that adopt past work on computationally efficient bounds of the absolute errors using quadrature formulas. The evaluation of the Gauss quadrature based estimates though, depends on a priori knowledge of extremal eigenvalues; and the upper bounds in particular that are useful as a stopping criterion fail in the absence of a reasonable underestimate of smallest eigenvalue. Estimators for relative errors in A-norm and their extension to errors in l2 norm are presented with numerical results. Estimating the relative error from the residue in an iterative solution is required for efficient solution of a large problem with even a moderately high condition. Specifically, in a problem of solving for vector x in Ax=b, the uncertainty between the strict upper bound in relative error [κ*||r(i)||/||b||] and its strict lower bound [||r(i)||/(κ*||b||)] is a factor of κ^2 (given residue r(i)= b-Ax(i) is the residual vector at ith iteration and κ the condition number of the square matrix A).