function test18
%TEST18 test iterative refinement after backslash
%
% Example:
%   test18
% See also: testall

%   Copyright 2006, Timothy A. Davis.
%   http://www.cise.ufl.edu/research/sparse


index = UFget ;
[ignore f] = sort (max (index.nrows, index.ncols)) ;
f = f (1:100) ;

clf
% f = f(1)

for i = f
    Prob = UFget (i) ;
    disp (Prob) ;
    A = Prob.A ;
    [m n] = size (A) ;
    if (~isreal (A) || m ~= n)
	continue
    end

    b = rand (n,1) ;
    x = A\b ;
    r = b - A*x ;
    x = x + A\r ;

    fprintf ('\n%6.2e to %6.2e\n', norm (r), norm (b-A*x)) ;
end