numbaでscipy.specialの関数をつかう

numbaでscipy.specialの関数を使うには、現状自分で拡張を書く必要がある。

support for scipy.special functions : feature request · Issue #3086 · numba/numba · GitHub

↓のノートブックでは球面調和関数用のAPIを作っている。 球面調和関数の返り値は複素数で、関数の返り値がcomplexの場合は実装がnumbaだけで閉じなくて難しいらしい。 ただ今の場合、陪ルジャンドル多項式(これは実)さえ求まればそこから球面調和関数も出せるので、上のissueで示されている方法に従えばよい。