01-Jul-2023 12:03:55

cg_squared_test
  MATLAB/Octave version 5.2.0.
  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| = 9.29744e-12
  Norm of error |x1-x2| = 9.29757e-12

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| = 4.13133e-08
  Norm of error |x1-x2| = 2.22011e-09

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

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

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

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

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

  Full iteration number = 2
  Norm of residual |Ax-b| = 2.9669e+34
  Norm of error |x1-x2| = 1.83085e+24

  Full iteration number = 3
  Norm of residual |Ax-b| = 1.98931e+44
  Norm of error |x1-x2| = 2.94941e+32

  Full iteration number = 4
  Norm of residual |Ax-b| = 2.40647e+56
  Norm of error |x1-x2| = 1.33794e+46

cg_squared_test
  Normal end of execution.

01-Jul-2023 12:03:55