[SciPy-dev] FAIL: test_yn_zeros, test_ynp_zeros, Negative-order Bessels

Pauli Virtanen pav@iki...
Mon Jan 26 03:23:37 CST 2009


Sun, 25 Jan 2009 20:03:26 +0100, Nils Wagner wrote:
> I do get three failures wrt. to special functions and python2.6

Too strict tolerances in the tests for the most part, I think. Will fix.
I presume this is a 64-bit Linux platform.

I'm in progress of addressing several bugs and improving the tests for 
the Bessel-related special functions in Scipy trunk. If someone wants to 
review the changesets I've committed recently, please go ahead.

> ======================================================================
> FAIL: test_yn_zeros (test_basic.TestBessel)
> ----------------------------------------------------------------------
[clip]
> Not equal to tolerance rtol=1e-19, atol=0

Too strict tolerance, I think. Some errors ~eps probably arise, depending 
on the compiler.

> ======================================================================
> FAIL: test_ynp_zeros (test_basic.TestBessel)
> ----------------------------------------------------------------------
[clip]
> AssertionError:
> Not equal to tolerance rtol=1e-07, atol=1e-15
> 
> (mismatch 100.0%)
>   x: array([  1.239e-10,  -8.119e-16,   3.608e-16, 5.898e-16,  
>   1.226e-15])
>   y: array(0)

Ditto. Except maybe for the first item. The question now is whether the 
problem is in the `yvp` routine or in `ynp_zeros`.

> ======================================================================
> FAIL: Negative-order Bessels
> ----------------------------------------------------------------------
[clip]
>      assert_tol_equal(iv(-0.5, 1   ), 1.231200214592967)
[clip]

Negative half-integers are a known failure of cephes/iv (also Scipy 
0.6.0), since hyperg has a pole at the corresponding points. But I think 
this test is commented out in the current SVN trunk, though.

I have a fix for this, but it's not committed yet. 

-- 
Pauli Virtanen



More information about the Scipy-dev mailing list