擅长:python、mysql、java
<p>如果有人还感兴趣的话,我找到了一个比特德·普利克(Ted Pudlik)快17倍的解决方案。我使用的事实是,n阶球面贝塞尔函数本质上是n阶标准贝塞尔函数的1/sqrt(x)乘以已经矢量化的n+1/2阶标准贝塞尔函数:</p>
<pre><code>import numpy as np
from scipy import special
sphj_bessel = lambda n, z: special.jv(n+1/2,z)*np.sqrt(np.pi/2)/(np.sqrt(z))
</code></pre>
<p>我得到了以下时间安排:</p>
^{pr2}$