[SciPy-user] odeint rtol and atol default values

Keith Goodman kwgoodman at gmail.com
Tue Jun 27 17:14:57 CDT 2006


On 6/27/06, Keith Goodman <kwgoodman at gmail.com> wrote:
> 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?

P.S. http://www.library.cornell.edu/nr/bookcpdf/c5-7.pdf



More information about the SciPy-user mailing list