14-May-2025 08:21:35

bicg_test():
  MATLAB/Octave version 6.4.0.
  Test bicg().
  To avoid conflicts with the built-in bicg() function,
  we refer to our code as "bicg_mine()"

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

  Norm of residual |Ax-b| = 2.38788e-09
  Norm of error |x1-x2| = 2.38788e-09

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

  Full iteration number = 1
  Norm of residual |Ax-b| = 0.000154786
  Norm of error |x1-x2| = 7.48222e-07

  Full iteration number = 2
  Norm of residual |Ax-b| = 2.99865e-06
  Norm of error |x1-x2| = 1.66455e-08

bicg_invol():
  Test bicg() on the invol matrix.
  Number of variables N = 20
  Estimated condition number = 4.8261e+20

  Full iteration number = 1
  Norm of residual |Ax-b| = 2.88953e+13
  Norm of error |x1-x2| = 2.44696e+14

  Full iteration number = 2
  Norm of residual |Ax-b| = 1.53197e+15
  Norm of error |x1-x2| = 6.0167e+14

bicg_test():
  Normal end of execution.

14-May-2025 08:21:35