python scipy ode dopri5 'larger nmax needed'

Question

在使用 scipy 0.13.0、ode(f).set_integrator('dopri5') 时，我收到错误消息 -

需要更大的 nmax

我在 ode.py 中查找了 nmax 但我看不到该变量。我猜想集成调用的数量超过了允许的默认值。

如何增加 nmax 值？

Answer 1

nmax 是指求解器将采用的最大内部步数。默认值为 500。您可以使用 set_integrator 方法的 nsteps 参数更改它。例如

ode(f).set_integrator('dopri5', nsteps=1000)

(Fortran 代码调用此 NMAX，显然 Fortran 名称已复制到“dopri5”求解器的 python 代码中的错误消息中。在 ode类 API，所有求解器(“dopri5”、“vode”、“lsoda”等)一致调用此求解器参数 nsteps，因此 scipy 应该将 python 代码中使用的错误消息更改为 nsteps.)