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.

As the scipy.signal.bessel docs note (permalink): The Bessel is inherently an analog filter. This function generates digital Bessel filters using the bilinear transform, which does not preserve the ...