This is the mail archive of the
glibc-bugs@sourceware.org
mailing list for the glibc project.
[Bug math/10709] sin() returns incorrect rounding
- From: "zimmerma+gcc at loria dot fr" <sourceware-bugzilla at sourceware dot org>
- To: glibc-bugs at sources dot redhat dot com
- Date: 30 Sep 2009 13:07:33 -0000
- Subject: [Bug math/10709] sin() returns incorrect rounding
- References: <20090930125053.10709.zimmerma+gcc@loria.fr>
- Reply-to: sourceware-bugzilla at sourceware dot org
------- Additional Comments From zimmerma+gcc at loria dot fr 2009-09-30 13:07 -------
This bug is analyzed and a fix is proposed at
<http://www.loria.fr/~zimmerma/papers/bug10709.pdf>.
In short, it suffices to replace the line:
return (res==res+1.025*cor)? res : slow1(x);
in file glibc/sysdeps/ieee754/dbl-64/s_sin.c by:
return (res==res+1.096*cor)? res : slow1(x);
Paul Zimmermann
--
http://sourceware.org/bugzilla/show_bug.cgi?id=10709
------- You are receiving this mail because: -------
You are on the CC list for the bug, or are watching someone who is.