[Octave-bug-tracker] [bug #53258] test pcg.m fail

[Carlo de Falco]

Follow-up Comment #17, bug #53258 (project octave):

Thinking better, even if the problem never occurred with Cristiano's code, I
agree with Rik that we can not expect error
to be below 1e-10 if we solve with a tolerance of 1e-6 on the residual.

So, if we really want to check error, either we should set tolerance to 0 to
make sure all N+1 iterations are run and then check the error but ignore the
convergence flag, or we should more carefully set the pcg tolerance parameter
so that it makes sense in comparison with the error bound used in the
assertion.

But probably Rik's idea of just checking the residual makes more sense here
...

Still I believe the new version of pcg is working better so it should be
checked in ASAP.



    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53258>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/