|
From: | John W. Eaton |
Subject: | Re: improving special functions in GSoC 2017 |
Date: | Sun, 26 Jun 2016 02:52:52 -0400 |
User-agent: | Mozilla/5.0 (X11; Linux x86_64; rv:38.0) Gecko/20100101 Icedove/38.5.0 |
On 06/26/2016 02:08 AM, Colin Macdonald wrote:
On 25/06/16 02:54, Marco Caliari wrote:"Special functions are expected by users to just work".So I just found that our "besselj" gives NaNs for values like "1e10". >> besselj(1, 1e9) ans = -5.21042264155388e-06 >> besselj(1, 1e10) ans = NaN + NaNi https://savannah.gnu.org/bugs/index.php?48316 That's not some wild exotic function: its a common Bessel function! Why are we not just calling some free-licensed library that nails these things?
Amos was the best thing I knew of at the time. If there is something better now, or a standard library solution, then perhaps we should be using it. Propose a patch and a test suite.
jwe
[Prev in Thread] | Current Thread | [Next in Thread] |