2003-01-20 Segher Boessenkool <segher@koffie.nl>

	* sysdeps/ieee754/dbl-64/sincos32.c (__mpcos): Really
	compute the cosine, not the sine, even if x > 0.8.

diff --git a/sysdeps/ieee754/dbl-64/sincos32.c b/sysdeps/ieee754/dbl-64/sincos32.c
index 989982e..a4f896a 100644 (file)
--- a/sysdeps/ieee754/dbl-64/sincos32.c
+++ b/sysdeps/ieee754/dbl-64/sincos32.c
@@ -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);