[GiNaC-list] Re: Airy functions

Richard B. Kreckel kreckel at ginac.de
Sat Dec 29 21:10:01 CET 2007

Alexei Sheplyakov wrote:
>>   // approximation for smallest number so that 1+tol > 1
>>   numeric tol=ex_to<numeric>(3*pow(ex(10),ex(-(Digits-2)))); 
> Not really.

Let me just add that plain CLN does come with a better function for 
that. From the manual:
`cl_F float_epsilon (float_format_t f)'
   Returns the smallest floating point number e > 0 such that 1+e != 1.

Richard B. Kreckel

More information about the GiNaC-list mailing list