[SciPy-dev] RFR 503, 849: more robust implementation of real Bessel I_v

Pauli Virtanen pav@iki...
Wed Feb 25 21:03:58 CST 2009

(For trying out the code review tool...)

Scipy bugs #503 and #849 are due to a non-robust implementation of Bessel 
I function in Cephes. The following changes address this, by using an 
implementation from the Boost library, converted to C:


A git branch is here:


I'm thinking that this code would be ready to be committed in. I'm not 
aware of bugs in it.

Pauli Virtanen

More information about the Scipy-dev mailing list