cexp error

Stephen Montgomery-Smith stephen at missouri.edu
Sat Sep 15 00:26:57 UTC 2012


On 09/14/2012 01:05 AM, Bruce Evans wrote:

> My tests also determined the exact minimum for all multiples up to the
> thresholds for "large" multiples in e_rem_pio2.c, except for ld128,
> to verify that there are enough bits in the special approximations for
> Pi/2 there, except for ld128.  The maximum multiples handled there are
> 2**28*Pi/2, except for ld128 they are 2**45*Pi/2.  2**45 is too many
> to check exhaustively.

But presumably one could prove a result to this effect using continued 
fractions, right?



More information about the freebsd-numerics mailing list