比較Boost.Odeint vs 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太子五階步進，dopri5。您可以使用 this Cython wrapper Boost.odeint（它並不公開所有由Boost.odeint提供的 步進器，但確實會公開dopri5），在Python中比較此算法的兩種實現的效率。

根據你對「測試性能」的定義，你也可以比較不同的算法，但是這顯然不同於相同算法的兩種實現。