math: fix comment in __rem_pio2f.c
authorSzabolcs Nagy <nsz@port70.net>
Wed, 12 Dec 2012 00:28:22 +0000 (01:28 +0100)
committerSzabolcs Nagy <nsz@port70.net>
Wed, 12 Dec 2012 00:28:22 +0000 (01:28 +0100)
src/math/__rem_pio2f.c

index 965dc46ac53a01b0a5dc05b4a2d82efda6f23ad5..198ee91b16bc08cbfa1112905152f2a1a318f3b7 100644 (file)
@@ -24,7 +24,7 @@
 
 /*
  * invpio2:  53 bits of 2/pi
- * pio2_1:   first  33 bit of pi/2
+ * pio2_1:   first 25 bits of pi/2
  * pio2_1t:  pi/2 - pio2_1
  */
 static const double
@@ -41,7 +41,7 @@ int __rem_pio2f(float x, double *y)
 
        GET_FLOAT_WORD(hx, x);
        ix = hx & 0x7fffffff;
-       /* 33+53 bit pi is good enough for medium size */
+       /* 25+53 bit pi is good enough for medium size */
        if (ix < 0x4dc90fdb) {  /* |x| ~< 2^28*(pi/2), medium size */
                /* Use a specialized rint() to get fn.  Assume round-to-nearest. */
                STRICT_ASSIGN(double, fn, x*invpio2 + 0x1.8p52);