24-Dec-2022 09:16:53

bicg_test():
  MATLAB/Octave version 4.2.2.
  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| = 7.3458e-11
  Norm of error |x1-x2| = 7.3458e-11

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| = 7.5412e-06
  Norm of error |x1-x2| = 4.31818e-08

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

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

  Full iteration number = 1
  Norm of residual |Ax-b| = 3.17592e+13
  Norm of error |x1-x2| = 4.42822e+15

  Full iteration number = 2
  Norm of residual |Ax-b| = 2.61161e+21
  Norm of error |x1-x2| = 3.9769e+16

bicg_test():
  Normal end of execution.

24-Dec-2022 09:16:54