07-Jan-2022 17:59:28

cg_squared_test
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2.
  Test cg_squared().

cg_squared_helmert
  Test cg_squared() on the Helmert matrix.
  Number of variables N = 10
  Estimated condition number = 1

  Norm of residual ||Ax-b|| = 1.15599e-11
  Norm of error ||x1-x2|| = 1.156e-11

cg_squared_lesp
  Test cg_squared() on the Lesp matrix.
  Number of variables N = 100
  Estimated condition number = 67.4852

  Full iteration number = 1
  Norm of residual ||Ax-b|| = 3.96347e-08
  Norm of error ||x1-x2|| = 1.8333e-09

  Full iteration number = 2
  Norm of residual ||Ax-b|| = 3.96347e-08
  Norm of error ||x1-x2|| = 1.8333e-09

  Full iteration number = 3
  Norm of residual ||Ax-b|| = 3.96347e-08
  Norm of error ||x1-x2|| = 1.8333e-09

  Full iteration number = 4
  Norm of residual ||Ax-b|| = 3.96347e-08
  Norm of error ||x1-x2|| = 1.8333e-09

cg_squared_invol
  Test cg_squared() on the invol matrix.
  Number of variables N = 20
  Estimated condition number = 8.21156e+20

  Full iteration number = 1
  Norm of residual ||Ax-b|| = 8.52547e+25
  Norm of error ||x1-x2|| = 1.56454e+16

  Full iteration number = 2
  Norm of residual ||Ax-b|| = 1.26774e+39
  Norm of error ||x1-x2|| = 3.76153e+28

  Full iteration number = 3
  Norm of residual ||Ax-b|| = 5.62734e+48
  Norm of error ||x1-x2|| = 1.78318e+37

  Full iteration number = 4
  Norm of residual ||Ax-b|| = 2.21496e+61
  Norm of error ||x1-x2|| = 4.26449e+49

cg_squared_test
  Normal end of execution.

07-Jan-2022 17:59:28