[SciPy-User] ODE integration : pre allocation

Fabrice Silva silva@lma.cnrs-mrs...
Wed Jul 4 06:07:15 CDT 2012


Hi folks,
I was wondering if the ode integrators in scipy.integrate.ode handle the
pre-allocation of the result of the rhs (i.e. y'=f(y,t) ).
A creation/allocation of the result array may lead to a strong
performance penalty, am I wrong ?

Is there a way to provide a reference to a pre-allocated array that
would only have to be filled ? For example, the callable signature would
be :
        f(t, y, y', *f_args)
where t, y and f_args are inputs, and y' is the output. When calling the
callable, y' would be a reference to the result array to be filled.

Any thought?
-- 
Fabrice Silva



More information about the SciPy-User mailing list