This is the mail archive of the cygwin mailing list for the Cygwin project.

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

Re: 1.5.24: sin() bug

Dmitry Karasik wrote:

> #include <math.h>
> #include <stdio.h>
> int main( int argc, char ** argv)
> {
>         double g = (double) 3.1415926535897900074;
>         printf("sin(%.10g)=%.10g\n", g, sin(g));
> }
> output is :
> sin(3.141592654)=3.231089149e-15
> whereas all other sin() implementation I could find ( freebsd, linux, msvc)
> report this:
> sin(3.141592654)=3.231085104e-015
> the difference is in 7th digit, and is significant for double precision.

This appears to be due to an approximation in __kernel_sin in which if
|x| < 2**-27 the approximation sin(x) ~ x is used.


This code seems to be completely untouched since the original import
into CVS in 2000, so I'm guessing it's a very old implementation.  I'm
not sure how we'd go about fixing this other than to look for a more
sophisticated implementation and contribute it to newlib.

In any case this is more of a newlib issue than a Cygwin issue as we are
just a consumer of this code, so followups should go to the newlib list.


Unsubscribe info:
Problem reports:

Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]