[SciPy-Dev] Consensus on the future of integrate.ode

Geoff Oxberry goxberry at mit.edu
Thu Sep 5 17:50:10 EDT 2013


Please, please don't implement anything in the SciPy ODE solvers with Boost.
odeint unless it's a separate interface that requires extreme generality.

Boost.odeint uses Boost.uBLAS for its linear algebra (see
http://headmyshoulder.github.io/odeint-v2/doc/boost_numeric_odeint/concepts/implicit_system.html),
and Boost.uBLAS is slow because it tries to be general (everything is
templated), and it has debugging features added by default.

A lot of users (myself included) don't need that generality, and would
prefer the extra performance that you'd get from a code that can use a
faster dense linear algebra library (like an optimized LAPACK, Eigen,
Elemental, maybe MTL4).

Geoff

On Thu, Sep 5, 2013 at 5:25 AM, Arnd Baecker <arnd.baecker at web.de> wrote:

> Could odeint be useful in this context?
>
>    http://headmyshoulder.github.io/odeint-v2/
>
> "Odeint is a modern C++ library for numerically solving Ordinary
> Differential Equations. It is developed in a generic way using Template
> Metaprogramming which leads to extraordinary high flexibility at top
> performance."
>
> Best, Arnd
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-dev
>



-- 
Geoffrey Oxberry, Ph.D., E.I.T.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20130905/bac26f00/attachment.html>


More information about the SciPy-Dev mailing list