This is the mail archive of the
glibc-bugs@sourceware.org
mailing list for the glibc project.
[Bug math/3976] libm rounding modes do not work correctly for many archs
- From: "c_keil at yahoo dot de" <sourceware-bugzilla at sourceware dot org>
- To: glibc-bugs at sources dot redhat dot com
- Date: 27 Oct 2008 08:17:06 -0000
- Subject: [Bug math/3976] libm rounding modes do not work correctly for many archs
- References: <20070206131410.3976.madcoder@debian.org>
- Reply-to: sourceware-bugzilla at sourceware dot org
------- 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 are on the CC list for the bug, or are watching someone who is.