2003-01-20 Segher Boessenkool <segher@koffie.nl>
authorroland <roland>
Mon, 20 Jan 2003 05:25:30 +0000 (05:25 +0000)
committerroland <roland>
Mon, 20 Jan 2003 05:25:30 +0000 (05:25 +0000)
* sysdeps/ieee754/dbl-64/sincos32.c (__mpcos): Really
compute the cosine, not the sine, even if x > 0.8.

sysdeps/ieee754/dbl-64/sincos32.c

index 989982e..a4f896a 100644 (file)
@@ -214,7 +214,7 @@ double __mpcos(double x, double dx) {
   __add(&a,&b,&c,p);
   if (x>0.8)
   { __sub(&hp,&c,&b,p);
-    __c32(&b,&a,&c,p);
+    __c32(&b,&c,&a,p);
   }
   else __c32(&c,&a,&b,p);     /* a = cos(x+dx)     */
   __mp_dbl(&a,&y,p);