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

[Nir Krakauer]

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

Hi Marco-- I added two more methods to your m-file to deal with positive a and
x without overflow. Now all the tests pass, albeit at fractional tolerances of
~1000*eps. I also streamlined the processing of scaled* options.

Based on our results so far, I think we should first implement gammainc as an
m-file, since the available Fortran programs don't work for all cases and are
difficult to debug. Once we have a working implementation, we can look at
increasing efficiency by calling C++ or Fortran for some branches.

(file #37083)
    _______________________________________________________

Additional Item Attachment:

File name: newgammainc.m                  Size:5 KB


    _______________________________________________________

Reply to this item at:

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

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