Сравнение Boost.Odeint vs Scipy.integrate.odeint?

Question

Я недавно опрокинул библиотеку boost.odeint, и я удивлен количеством возможностей и настраиваемости. Однако, широко используя scipy.integrate.odeint (который по существу является оберткой для ODEPACK в fortran), мне интересно, как их производительность сравнивается. Я знаю, что boost.odeint также поставляется с распараллеливанием, что невозможно с помощью scipy (насколько я знаю), что значительно увеличило бы производительность, но я прошу об одном ядре. Однако, поскольку в этом случае мне придется обернуть boost.odeint (используя cython или boost.python) в python, возможно, кто-то из вас уже это сделал? Это было бы большим достижением, поскольку все возможности анализа намного более продвинуты в python.

Answer 1

Насколько я могу судить из сравнения списков доступных степперы для Boost.odeint и scipy.integrate.ode, единственный алгоритм реализуется как это Дорманда-Принс пятого порядка степпер, dopri5. Вы могли бы сравнить эффективность двух реализаций этого алгоритма в Python с помощью this Cython wrapper в Boost.odeint (он не выставляет все степперы , предоставленные Boost.odeint, но обнажает dopri5).

В зависимости от вашего определения «производительности тестирования» вы также можете сравнить различные алгоритмы, но это, очевидно, не то же самое, что сравнение двух реализаций того же алгоритма.