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

[Lachlan Andrew]

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

The the handling of arrays is really nice -- not a single for loop!

Here is an updated version that copies most of the documentation from the
existing C++ implementation.

I hadn't realised previously that the current implementation is in C++.  That
presumably will make it much faster than the current one.  Would you propose
to replace the existing C++ function by this one, or would you want to call
the C++ function for tail=="lower"?  Or should this be translated this into
C++, now that you have the algorithms selected?

(file #37225)
    _______________________________________________________

Additional Item Attachment:

File name: ngammainc.m                    Size:9 KB


    _______________________________________________________

Reply to this item at:

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

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