G-2016-70
Estimates of the 2-norm forward error for SYMMLQ and CG
, , and BibTeX reference
For positive definite linear systems (or semidefinite consistent systems), we use Gauss-Radau quadrature to obtain a cheaply computable upper bound on the 2-norm error of SYMMLQ iterates. The close relationship between SYMMLQ and CG iterates lets us construct an upper bound on the 2-norm error for CG. For indefinite systems, the upper bound becomes an estimate of the 2-norm SYMMLQ error. Numerical experiments demonstrate that the bounds and estimates are remarkably tight.
Published September 2016 , 15 pages
Research Axis
Research application
Publication
Jan 2019
, , and
SIAM Journal on Matrix Analysis and Applications, 40(1), 235–253, 2019
BibTeX reference