[SciPy-user] odeint rtol and atol default values

Keith Goodman kwgoodman at gmail.com
Tue Jun 27 17:08:54 CDT 2006


On 6/27/06, Steve Schmerler <elcorto at gmx.net> wrote:
> David M. Cooke wrote:
>
> >>>> What are the default values for rtol and atol in
> >>>> scipy.integrate.odeint?
> >>>>
>
> >
> > Digging through __odepack.h finds that the defaults are 1.49012e-8.
> >
>
> Thanks. Maybe doc string should mention this.
>
> Not that I would care to much, but this seems a rather random value at
> the first glance. Is there a special reason for this val?

It looks like sqrt(eps) where eps is the machine precision.

BTW, can scipy return eps?



More information about the SciPy-user mailing list