[Bug math/3976] libm rounding modes do not work correctly for many archs
------- Additional Comments From c_keil at yahoo dot de 2008-10-27 08:17 -------
At least on my Core2 Duo it's also not working with still other values (32bit
version is ok):
gcc -m64 -lm -o a a.c; ./a 1 2
N: exp(1) = 2.7182818284590451
Z: exp(1) = 2.7182818284590451
D: exp(1) = 2.7182818284590451
U: exp(1) = 7.1387612927397726
N: exp(2) = 7.3890560989306504
Z: exp(2) = 4.0037745305985499
D: exp(2) = 4.0037745305985499
U: exp(2) = 7.3890560989306504
As far as I dug into it, the 32bit and 64bit versions use other code. 64bit
comes from the IBM Accurate Mathematical Library. For exp
sysdeps/ieee754/dbl-64/e_exp.c produces the wrong results.
--
What |Removed |Added
----------------------------------------------------------------------------
CC| |c_keil at yahoo dot de
http://sourceware.org/bugzilla/show_bug.cgi?id=3976
------- You are receiving this mail because: -------
You reported the bug, or are watching the reporter.
You are on the CC list for the bug, or are watching someone who is.
Reply to: