c++ - Boost OdeInt 中是否有 lsode 模拟？

Question

里面有boost 1.53 和 OdeInt。它有一系列数值 ode 求解器。但是我想知道其中是否有lsode模拟 - 对于您已保证分析解决方案的情况？

Answer 1

我不完全理解你的问题。LSODE 求解 dx/dt = f(x,t) [1] 形式的任何通用 ODE，就像 odeint 中的方法一样。然而，LSODE 检查问题的刚度并选择正确的方案（非刚度的预测校正器和刚度问题的反向微分）。odeint 不提供这种自动选择，而是依赖于用户为问题选择正确的方案。据我所知，这两种实现都没有使用可能的现有分析解决方案。

[1] https://computation.llnl.gov/casc/odepack/odepack_home.html

Answer 2

大多数 ODE 求解器确实暗示了对原始 LSODE（Fortran 实现）的某种限制。目前我所知道的最接近的替代方案是来自 GNU 的 Octave。