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.