[SciPy-user] scipy.test() floating point exception [old: svn commands]

[Steve Schmerler]

Hi

with the latest svn versions I get:

##########################################################################################
scipy.test(verbostiy=2)

[...]

check_fdtridfd (scipy.special.basic.test_basic.test_cephes) ... ok
check_fresnel (scipy.special.basic.test_basic.test_cephes) ... ok
check_gamma (scipy.special.basic.test_basic.test_cephes) ... ok
check_gammainc (scipy.special.basic.test_basic.test_cephes) ... ok
check_gammaincc (scipy.special.basic.test_basic.test_cephes) ... ok
check_gammainccinv (scipy.special.basic.test_basic.test_cephes) ... ok
check_gammaln (scipy.special.basic.test_basic.test_cephes) ... ok
check_gdtr (scipy.special.basic.test_basic.test_cephes) ... ok
check_gdtrc (scipy.special.basic.test_basic.test_cephes) ... ok
check_gdtria (scipy.special.basic.test_basic.test_cephes) ... ok
check_gdtrib (scipy.special.basic.test_basic.test_cephes) ... ok
check_gdtrix (scipy.special.basic.test_basic.test_cephes) ... ok
check_hankel1 
(scipy.special.basic.test_basic.test_cephes)Gleitkomma-Ausnahme
##########################################################################################

The test runs fine and crashes when the hankel1 function is called.

Calls to the other hankel function (defined in the _cephes module) also 
end up in float exceptions:

scipy.special._cephes.hankel1(1,1)
scipy.special._cephes.hankel1e(1,1)
scipy.special._cephes.hankel2(1,1)
scipy.special._cephes.hankel2e(1,1)

Outcommenting the tests of these functions in 
special/tests/test_basic.py makes place for more exceptions when calling 
other functions (i0e, i1e, ive, jve, k0, and probably many more), all 
from the _cephes module.

Any ideas?

cheers,
steve

-- 
Should array indices start at 0 or 1? My compromise of 0.5 was rejected 
without, I thought, proper consideration.
-- Stan Kelly-Bootle