X-Git-Url: http://nsz.repo.hu/git/?p=musl;a=blobdiff_plain;f=src%2Fmath%2F__rem_pio2f.c;h=198ee91b16bc08cbfa1112905152f2a1a318f3b7;hp=965dc46ac53a01b0a5dc05b4a2d82efda6f23ad5;hb=9c6b1de0fb11d6e2927a19bbaf6345aedcc84297;hpb=b69f695acedd4ce2798ef9ea28d834ceccc789bd diff --git a/src/math/__rem_pio2f.c b/src/math/__rem_pio2f.c index 965dc46a..198ee91b 100644 --- a/src/math/__rem_pio2f.c +++ b/src/math/__rem_pio2f.c @@ -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);