0

我有 2*m+3 个STIFF常微分方程要求解。我尝试使用 MATLAB ode15s for m=1 解决它，它工作正常。但我想使用 Sundials 包 CVODE 来求解方程。在尝试这样做时，我使用了“向后微分公式”和牛顿迭代。我不提供雅可比并允许以数字方式计算它。但它不起作用并显示错误：

[CVODE WARNING] CVode
 Internal t = 0 and h = 0 are such that t + h = t on the next step. The solver will continue anyway.

[CVODE ERROR] CVode
 At t = 0 and h = 0, the correction convergence test failed repeatedly or with :h: = hmin.

SUNDIALS_ERROR: CVode<> failed with flag -4
  1. 我相信 CVODE 使用与 ode15s 相同的反向微分。那为什么它不起作用？
  2. 我应该尝试在 CVODE 中使用带有预处理的 Krylov 求解器吗？

期待任何帮助。谢谢你。

4

0 回答 0