07-Jan-2022 20:33:22

gmres_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2.
  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 56)
  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 56)
  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 56)
  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.

07-Jan-2022 20:33:22