A Variant of the ORTHOMIN(2) Method for Singular Linear Systems

For singular linear systems A x=b, ORTHOMIN(2) is known theoretically to attain the minimum residual min _xR nb–A x₂ under a certain condition. However, in the actual computation with finite precision arithmetic, the residual is often observed to be reduced further than the theoretically expected level. Therefore, we propose a variant of ORTHOMIN(2), which is mathematically equivalent to the original ORTHOMIN(2) method, but uses recurrence formulas that are different from those of ORTHOMIN(2); they contain alternative expressions for the auxiliary vector and the recurrence coefficients. Although our implementation has the same computational costs as ORTHOMIN(2), numerical experiments on singular systems show that our implementation is more accurate and less affected by rounding errors than ORTHOMIN(2).