fpu/e_expl.c for i686
Stephen L Moshier
moshier@mediaone.net
Mon May 7 12:57:00 GMT 2001
> expl (1000) is .19e435 and can only be represented with long doubles.
> How much would the error be for expl (100)?
Instead of exp(x) the program computes exp(x(1 + dx)) where dx
is the relative input error introduced, in this case by the range
reduction method.
exp(x(1 + dx)) = exp(x) (1 + x dx + ...)
So if x is 100 and dx is 1/2 ulp, the relative error would be
about 50 ulps (long double ulps).
More information about the Libc-alpha
mailing list