[SciPy-dev] FAIL: test_yn_zeros, test_ynp_zeros, Negative-order Bessels
Pauli Virtanen
pav at iki.fi
Mon Jan 26 04:28:47 EST 2009
Mon, 26 Jan 2009 09:03:19 +0100, Nils Wagner wrote:
[clip]
> ======================================================================
> FAIL: test_iv_cephes_vs_amos (test_basic.TestBessel)
> ----------------------------------------------------------------------
[clip]
> Not equal to tolerance rtol=1e-08, atol=1e-305 (-120, -11)
> (mismatch 100.0%)
> x: array(1.3384173609003782e-110)
> y: array((1.3384173859242368e-110+0j))
The way cephes/iv computes the value of Iv is not too accurate for large
orders or large x. I have some code to compute the result in this range
from proper asymptotic expansions, but also it isn't committed yet.
--
Pauli Virtanen
More information about the SciPy-Dev
mailing list