Is it possible to expose the low-level special functions from scipy.special.cython_special as C symbols in a shared library (e.g., libscipy_special.so), similar to how libraries like GSL (libgsl.so) ...
Spherical Bessel functions break when called for complex arguments with zero imaginary part when checking the type of the argument. Changing any (z<0) to any (real (z+0j)<0) fixes the problem.