cvs commit: src/lib/msun/src e_expf.c

[Bruce Evans]

On Wed, 30 Nov 2005, Bruce Evans wrote:

>  Modified files:
>    lib/msun/src         e_expf.c
>  Log:
>  ...  The end result is
>  that expf() is now perfectly rounded (to nearest) except in 52561 cases
>  instead of except in 67027 cases, and the maximum error is 0.5013 ulps
>  instead of 0.5023 ulps.

Oops, this is only for i386 where there is considerable extra precision.
On A64, the maximum error actually mysteriously increased a little (to
0.9101 ulps), and there are a lot more imperfectly rounded cases
(17928157 now).  The algorithm is only supposed to give < 1 ulps.

Bruce