比较 Boost.Odeint 与 Scipy.integrate.odeint？

Question

我最近偶然发现了 boost.odeint 库，我对它的可能性和可配置性的数量感到惊讶。但是，在广泛使用 scipy.integrate.odeint （本质上是 fortran 中 ODEPACK 的包装器）之后，我想知道它们的性能比较如何。我知道 boost.odeint 也带有并行化，这对于 scipy （据我所知）来说是不可能的，这会大大提高性能，但我要求的是单核案例。 但是，由于在这种情况下我必须将 boost.odeint（使用 cython 或 boost.python）包装到 python 中，也许你们中的某些人已经这样做了？这将是一个伟大的成就，因为所有的分析可能性都在 python 中更高级。

Answer 1

据我所知，通过比较可用步进器列表 Boost.odeint 和 scipy.integrate.ode，两者实现的唯一算法是 Dormand-Prince 五阶步进电机，dopri5。你可以比较 在 Python 中使用该算法的两种实现的效率 this Cython wrapper 到 Boost.odeint（它不会暴露所有 Boost.odeint 提供的步进器，但确实暴露了 dopri5)。

根据您对“测试性能”的定义，您还可以比较 不同的算法，但这显然与比较不同 同一算法的两种实现。