擅长:python、mysql、java
<p>我想你说得对。<code>GNLSE_RHS</code>显然是你的瓶颈。在</p>
<pre><code>52007 188.766 0.004 387.243 0.007 gnlse.py:234(GNLSE_RHS)
</code></pre>
<p>问题是你的<code>percall</code>非常短。我猜<code>GNLSE_RHS</code>是一个由<code>scipy</code>包装的fortran函数。如果是这样的话,就很难将其简化。在</p>
<p>我解决薛定谔方程(主要是虚时间传播)的方法一直是在选择最终用C实现之前,先用Python来理解算法</p>