This is the mail archive of the
newlib@sourceware.org
mailing list for the newlib project.
inaccuracy in cos()
- From: Martin Koeppe <mkoeppe at gmx dot de>
- To: newlib at sources dot redhat dot com
- Date: Fri, 19 Oct 2007 13:38:57 +0200 (CEST)
- Subject: inaccuracy in cos()
Hi,
just found some inaccuracy in mathfp/cos(): cos(0.0) returns
0.999999999994129362690387097246 instead of 1.0. For 64-bit double
precision (53 bit mantissa, 18 decimal digits) this could be
considered as wrong.
Also I missed the fast assembler versions in machine/i386 of at least
sin(), cos(), and sqrt(). Assembler instructions exist for these
since at least i387 AFAIK.
Martin