如何解决solve_ivp 集成的边界
我有一个带有物理解释的 ODE（用于相空间密度），其中值应始终为非负。不幸的是，solve_ivp 使用了足够大的时间步长，因此这些值变为负数。强制较小的时间步长会导致计算时间的大量增加。所以我正在寻找其他方法来强制执行非否定性：
- 我查看了 solve_bvp，但它似乎只能强制执行一些线性 S*y = 0 而不是 y > 0
- 有没有办法改变solve_ivp集成中的哪些步骤被接受/放弃？
- 或者，当值接近于零时，有没有办法调整步长？
感谢您的帮助，
尼克
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。
