我的问题是关于当前的scipy颂歌求解器。从scipy doc page,它们的用法是:矢量化SciPy颂歌求解器
# A problem to integrate and the corresponding jacobian:
from scipy.integrate import ode
y0, t0 = [1.0j, 2.0], 0
def f(t, y, arg1):
return [1j*arg1*y[0] + y[1], -arg1*y[1]**2]
def jac(t, y, arg1):
return [[1j*arg1, 1], [0, -arg1*2*y[1]]]
# The integration:
r = ode(f, jac).set_integrator('zvode', method='bdf', with_jacobian=True)
r.set_initial_value(y0, t0).set_f_params(2.0).set_jac_params(2.0)
t1 = 10
dt = 1
while r.successful() and r.t < t1:
r.integrate(r.t+dt)
print("%g %g" % (r.t, r.y))
我的问题是:它使用了很多蟒蛇循环(while循环)基本上减慢程序下来。我可以尝试编写C代码并使用ctypes使其更快,但我无法在scipy中访问像dopri5这样的漂亮算法(除非我自己实现它)。
是否有任何矢量化的编码方式,使其更快?
谢谢!
非常感谢!这对我来说是完全意义上的。 –