[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #47800] gammainc rounds down to zero if output
|
From: |
Nir Krakauer |
|
Subject: |
[Octave-bug-tracker] [bug #47800] gammainc rounds down to zero if output is below eps |
|
Date: |
Thu, 28 Apr 2016 18:11:04 +0000 |
|
User-agent: |
Mozilla/5.0 (Macintosh; Intel Mac OS X 10_11_3) AppleWebKit/601.4.4 (KHTML, like Gecko) Version/9.0.3 Safari/601.4.4 |
Follow-up Comment #1, bug #47800 (project octave):
The problem is really in
http://hg.savannah.gnu.org/hgweb/octave/file/@/libinterp/corefcn/gammainc.cc
There, gammainc(a, x, 'upper') is computed as 1 - gammainc(a, x), which
expectedly evaluates to 0 if gammainc(a, x) is within eps of 1.
If we want more accurate results, we need a separate function to compute the
complement.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?47800>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/