|
From: | Colin Macdonald |
Subject: | [Octave-bug-tracker] [bug #47800] gammainc(x, a, "upper") rounds down to zero if output is below eps |
Date: | Fri, 06 May 2016 05:10:11 +0000 |
User-agent: | Mozilla/5.0 (X11; Fedora; Linux x86_64; rv:45.0) Gecko/20100101 Firefox/45.0 |
Follow-up Comment #11, bug #47800 (project octave): In case it helps with development, I've just added gammainc to Symbolic. Here's a useful recipe for experimenting: >> xs = sym (20); x = double (xs); a = 1; >> A = gammainc (x, a, 'upper') A = 2.0612e-09 >> B = double (gammainc (xs, a, 'upper')) B = 2.0612e-09 >> relerr = (A-B)/abs(B) relerr = -2.0034e-08 I also put some 'xtest' BIST into '@sym/gammainc' that should pass once this bug is fixed. I'm worried about the state of special functions in Octave :( _______________________________________________________ Reply to this item at: <http://savannah.gnu.org/bugs/?47800> _______________________________________________ Message sent via/by Savannah http://savannah.gnu.org/
[Prev in Thread] | Current Thread | [Next in Thread] |