On Feb 8, 2009, at 12:23 PM, Pauli Virtanen wrote:

However, the Boost library seems to have good implementations Bessel
(and some other) special functions:

http://svn.boost.org/svn/boost/trunk/boost/math/special_functions/detail/
http://www.boost.org/doc/libs/1_37_0/libs/math/doc/sf_and_dist/html/math_toolkit/special.html

Also the license seems Scipy-compatible:

http://www.boost.org/LICENSE_1_0.txt

So, I'd like to bring these over to Scipy, to replace some of the Cephes
routines.

The only problem is that being in Boost, they are written in C++, and I
guess we can't make Scipy to depend on it.

What about the GSL implementation of the Bessel function?  That's already in C and seems, in some sense, a more natural companion library to SciPy than boost.