[SciPy-user] odeint rtol and atol default values
doug-scipy at sadahome.ca
Tue Jun 27 14:39:33 CDT 2006
I've been happily using the ode class with the vode integrator for a few
months now. I rewrote a model from Octave into Python/NumPy/SciPy.
Agreement between the Python/ode/vode code and the Octave one was good. The
model is substantially faster in SciPy than it was in Octave, but there are
a lot of factors that changed (processor, OS, etc.).
On 6/25/06, Steve Schmerler <elcorto at gmx.net> wrote:
> Nils Wagner wrote:
> > On Sun, 25 Jun 2006 19:17:16 +0200
> > Steve Schmerler <elcorto at gmx.net> wrote:
> >> Hi
> >> What are the default values for rtol and atol in
> >> scipy.integrate.odeint?
> > You can find them in scipy/Lib/integrate/ode.py
> > rtol=1e-6,atol=1e-12
> Sure? What you mention is in site-packages/scipy/integrate/ode.py where
> the class ode is defined. But I'm using integrate.odeint (in
> site-packages/scipy/integrate/odepack.py) and there is nothing mentioned
> about it. odeint is there called like this
> import _odepack
> output = _odepack.odeint(func, y0, t, args, Dfun, col_deriv, ml, mu,
> full_output, rtol, atol, tcrit, h0,
> ixpr, mxstep, mxhnil, mxordn, mxords)
> BTW, does anyone have experiences with the ode class (and it's vode (or
> cvode?) solver) rather than odeint's lsoda?
> Random number generation is the art of producing pure gibberish as
> quickly as possible.
> SciPy-user mailing list
> SciPy-user at scipy.net
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-user