[SciPy-user] odeint rtol and atol default values

Travis Oliphant oliphant.travis at ieee.org
Tue Jun 27 17:10:44 CDT 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