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

[Lachlan Andrew]

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

Hmmm...  Before this can be applied, it will definitely need to handle
vectors, even if only by loops.

There are many special cases in this code.  I'd imagine that at least some of
them can be vectorised efficiently.  The iterative components are probably the
most expensive, and they seem to be vectorisable, although iteration would
need to continue until all had converged, and so the actual FLOP count could
be much higher than currently.

I'll add it to my to-do list, but will need reminding...

    _______________________________________________________

Reply to this item at:

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

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