scipy.special
次数 n と引数 x のベッセル関数が xjv(n,x)
でベクトル化されていることに気付きました。
In [14]: import scipy.special as sp
In [16]: sp.jv(1, range(3)) # n=1, [x=0,1,2]
Out[16]: array([ 0., 0.44005059, 0.57672481])
しかし、対応する球面ベッセル関数のベクトル化された形式はありませんsp.sph_jn
。
In [19]: sp.sph_jn(1,range(3))
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-19-ea59d2f45497> in <module>()
----> 1 sp.sph_jn(1,range(3)) #n=1, 3 value array
/home/glue/anaconda/envs/fibersim/lib/python2.7/site-packages/scipy/special/basic.pyc in sph_jn(n, z)
262 """
263 if not (isscalar(n) and isscalar(z)):
--> 264 raise ValueError("arguments must be scalars.")
265 if (n != floor(n)) or (n < 0):
266 raise ValueError("n must be a non-negative integer.")
ValueError: arguments must be scalars.
さらに、球ベッセル関数は、1 回のパスで N のすべての次数を計算します。したがって、n=5
引数 にベッセル関数が必要な場合x=10
、n=1,2,3,4,5 が返されます。実際には、jn とその導関数を 1 回のパスで返します。
In [21]: sp.sph_jn(5,10)
Out[21]:
(array([-0.05440211, 0.07846694, 0.07794219, -0.03949584, -0.10558929,
-0.05553451]),
array([-0.07846694, -0.0700955 , 0.05508428, 0.09374053, 0.0132988 ,
-0.07226858]))
API にこの非対称性が存在するのはなぜですか。ベクトル化された、または少なくともより迅速に (つまり cython で) 球状ベッセル関数を返すライブラリを知っている人はいますか?