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

Question

我最近偶然发现了 boost.odeint 库，我对可能性和可配置性的数量感到惊讶。但是，在广泛使用 scipy.integrate.odeint(它本质上是 Fortran 中 ODEPACK 的包装器)之后，我想知道它们的性能如何比较。我知道 boost.odeint 还带有并行化，这对于 scipy(据我所知)是不可能的，这会大大 boost 性能，但我要求的是单核案例。但是，由于在那种情况下我必须将 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)。

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