微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

SciPy:集成期间solve_ivp的状态

如何解决SciPy:集成期间solve_ivp的状态

我正在使用scipy.integrate.solve_ivp在Python中运行长时间的ODE集成。在例程运行时,是否可以访问集成状态或检查例程位于哪个集成步骤?我的整合花费了比预期更长的时间,我想知道整合者是停留在某个步骤还是单个步骤花费的时间真的很长。

用于将来的任务;如果我将与solve_ivp的集成拆分为子间隔以在其间打印状态消息,这是否会干扰某些求解器的步长适应性?

感谢您的反馈!

解决方法

存在一个GitHub pull请求,用于向Solve_ivp添加一个详细选项,但这尚未完成。您可以自己修改scipy的resolve_ivp函数(应该很简单)来实现它,也可以只打印求解器给ODE函数指定的时间t。我就是做这个的。如果您的系统不是太小,那么您就不会因为打印而浪费很多时间。

按照建议的工作方式拆分集成,但是,如果每隔几个时间拆分一次,则由于求解器重新启动每个图块,您将浪费时间。这种影响与隐式算法有关,因为它们在每次启动时都会重新计算系统的雅可比行列式。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。