cvs commit: src/lib/msun/src e_expf.c
bde at zeta.org.au
Wed Nov 30 12:48:24 GMT 2005
On Wed, 30 Nov 2005, Bruce Evans wrote:
> Modified files:
> lib/msun/src e_expf.c
> ... 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.
More information about the cvs-src