• STFC home
  The Science and Technology Facilities Council
• HSL home
  About the HSL Mathematical Software Library
• HSL 2013
  Our current release
• HSL archive
  Older packages available at no cost
• Academic access
  Free access for teaching and research
• The Numerical Analysis Group
  About the research group behind HSL
• Contact us
  Our contact details

MA60 Iterative refinement and error estimation

This subroutine performs iterative refinement and provides information on the error in the resulting solution of a set of $n$ sparse linear equations $Ax=b$, either as an upper bound on the actual error (forward error) or as a bound on the perturbation to the matrix elements which would make the solution obtained the exact solution of a perturbed problem with the same sparsity pattern as the original matrix (backward error). All estimates are provided in the $\infty$-norm, that is $\left|\right|x|{|}_{\infty }=\underset{i}{max}|x_i|$.

There is an option for avoiding iterative refinement, in which case bounds on the errors associated with the given solution are provided.