[Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps

[Marco Caliari]

Follow-up Comment #37, bug #47800 (project octave):

Excellent, Nir. I see that you do not remove converged values in
__gammainc_lentz as you did for gammainc_s (jj subindex). Isn't it dangerous?
For instance, with your previous file #37329, if I do


gammainc(rand(4),rand(4))


quite often I have no convergence. The problem seems to me that at a certain
iteration, some values have converged, but at the next iteration their error
is a bit above the tolerance and other converged, and so on. However, I cannot
reproduce this behavior with your new file #37342.

    _______________________________________________________

Reply to this item at:

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

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