[SciPy-user] odeint rtol and atol default values
Travis Oliphant
oliphant.travis at ieee.org
Tue Jun 27 18:10:44 EDT 2006
Steve Schmerler 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's the sqrt(epsilon) where epsilon is the level of precision on C
double for IEEE double-precision standard.
-Travis
More information about the SciPy-User
mailing list