Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

13-May-2025 14:27:14

gmres_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7.
  Test gmres().

gmres_helmert():
  Test gmres() on the Helmert matrix.
  Number of variables N = 10
  Estimated condition number = 1

  Norm of residual |Ax-b| = 7.4476e-16
  Norm of error |x1-x2| = 8.45521e-16

gmres_lesp():
  Test gmres() on the Lesp matrix.
  Number of variables N = 100
  Estimated condition number = 67.4852

  Full iteration number = 1
  Norm of residual |Ax-b| = 3.53595e-07
  Norm of error |x1-x2| = 1.52433e-08

  Full iteration number = 2
  Norm of residual |Ax-b| = 2.19289e-07
  Norm of error |x1-x2| = 1.01868e-08

  Full iteration number = 3
  Norm of residual |Ax-b| = 1.5477e-07
  Norm of error |x1-x2| = 7.47486e-09

  Full iteration number = 4
  Norm of residual |Ax-b| = 1.2063e-07
  Norm of error |x1-x2| = 6.04709e-09

gmres_invol():
  Test gmres() on the invol matrix.
  Number of variables N = 20
  Estimated condition number = 8.21156e+20
[Warning: Matrix is close to singular or badly scaled. Results may be
inaccurate. RCOND =  1.586431e-16.] 
[> In gmres (line 83)
In gmres_invol (line 55)
In gmres_test (line 40)
In run (line 91)
] 

  Full iteration number = 1
  Norm of residual |Ax-b| = 1.05627e+06
  Norm of error |x1-x2| = 4.77326e+08
[Warning: Matrix is close to singular or badly scaled. Results may be
inaccurate. RCOND =  1.115746e-16.] 
[> In gmres (line 83)
In gmres_invol (line 55)
In gmres_test (line 40)
In run (line 91)
] 

  Full iteration number = 2
  Norm of residual |Ax-b| = 1.18384e+06
  Norm of error |x1-x2| = 5.30637e+08

  Full iteration number = 3
  Norm of residual |Ax-b| = 540157
  Norm of error |x1-x2| = 5.06188e+08
[Warning: Matrix is close to singular or badly scaled. Results may be
inaccurate. RCOND =  2.248226e-17.] 
[> In gmres (line 83)
In gmres_invol (line 55)
In gmres_test (line 40)
In run (line 91)
] 

  Full iteration number = 4
  Norm of residual |Ax-b| = 3.34684e+06
  Norm of error |x1-x2| = 1.31099e+09

gmres_test():
  Normal end of execution.

13-May-2025 14:27:14