08-Jan-2023 07:15:04

gmres_test():
  MATLAB/Octave version 4.2.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| = 4.67745e-16
  Norm of error |x1-x2| = 4.40641e-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| = 2.90168e-07
  Norm of error |x1-x2| = 9.20567e-09

  Full iteration number = 2
  Norm of residual |Ax-b| = 1.81042e-07
  Norm of error |x1-x2| = 6.22492e-09

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

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

gmres_invol():
  Test gmres() on the invol matrix.
  Number of variables N = 20
  Estimated condition number = 4.82609e+20

  Full iteration number = 1
  Norm of residual |Ax-b| = 5135.42
  Norm of error |x1-x2| = 3.1935e+06

  Full iteration number = 2
  Norm of residual |Ax-b| = 5552.5
  Norm of error |x1-x2| = 3.1935e+06

  Full iteration number = 3
  Norm of residual |Ax-b| = 5552.5
  Norm of error |x1-x2| = 3.1935e+06

  Full iteration number = 4
  Norm of residual |Ax-b| = 5552.5
  Norm of error |x1-x2| = 3.1935e+06

gmres_test():
  Normal end of execution.

08-Jan-2023 07:15:04