floatscan: fix incorrect count of leading nonzero digits

[musl] / src / math / __cosl.c

diff --git a/src/math/__cosl.c b/src/math/__cosl.c
index 9ea51ec..9d32576 100644 (file)
--- a/src/math/__cosl.c
+++ b/src/math/__cosl.c
@@ -41,19 +41,8 @@
  * almost for free from the complications needed to search for the best
  * higher coefficients.
  */
-static const double one = 1.0;
-
-// FIXME
-/* Long double constants are slow on these arches, and broken on i386. */
-static const volatile double
-C1hi = 0.041666666666666664,            /*  0x15555555555555.0p-57 */
-C1lo = 2.2598839032744733e-18;          /*  0x14d80000000000.0p-111 */
-#define C1      ((long double)C1hi + C1lo)
-
-#if 0
 static const long double
 C1 =  0.0416666666666666666136L;        /*  0xaaaaaaaaaaaaaa9b.0p-68 */
-#endif
 
 static const double
 C2 = -0.0013888888888888874,            /* -0x16c16c16c16c10.0p-62 */
@@ -70,7 +59,7 @@ long double __cosl(long double x, long double y)
        z  = x*x;
        r  = z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*C7))))));
        hz = 0.5*z;
-       w  = one-hz;
-       return w + (((one-w)-hz) + (z*r-x*y));
+       w  = 1.0-hz;
+       return w + (((1.0-w)-hz) + (z*r-x*y));
 }
 #endif