我试图用python解决耦合的复杂ode系统西皮·奥德带zvode积分器。但是一旦我运行代码,这个错误消息就会出现。在
ZVODE-- At T(=R1) and step size H(=R2), the error test failed repeatedly or with abs(H) = HMIN. In above, R1 = 0.1018805870139D-15 R2 = 0.2392554739952D-22
我确实看过FORTRAN的源代码,但无法理解它的含义。在
如有任何帮助,我们将不胜感激。在
编辑:代码已包含在内。 我还试着打印出一些值,还为使用简单Euler方法的集成编写了单独的代码。从这些我有一种感觉,误差可能是由于数值超出范围,即大于10^308。(可能是由于某些参数的错误)。谁能确认一下吗?在
^{pr2}$
这意味着你的系统是僵硬的,步长控制器的启发式算法计算出它需要非常小的步长来保证所需的误差边界,但是步长变得非常小,因此所需的步数如此之大,以至于浮点噪声的积累变得更为主要,这意味着控制器将失去对误差累积的控制。似乎为了避免这种情况,控制器将步长限制在}值的一部分。在
2e-7
,比sqrt(mu)
,即{相关问题 更多 >
编程相关推荐