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

Pauli Virtanen pav at iki.fi
Wed Feb 25 22:03:58 EST 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:

	http://codereview.appspot.com/20078

A git branch is here:

	http://github.com/pv/scipy-work/tree/ticket-503-special-iv-fix

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