diff options
Diffstat (limited to 'newlib/libm/common/sf_lrint.c')
-rw-r--r-- | newlib/libm/common/sf_lrint.c | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/newlib/libm/common/sf_lrint.c b/newlib/libm/common/sf_lrint.c index a58508f02..427bbb655 100644 --- a/newlib/libm/common/sf_lrint.c +++ b/newlib/libm/common/sf_lrint.c @@ -28,6 +28,10 @@ static const float #else static float #endif +/* Adding a float, x, to 2^23 will cause the result to be rounded based on + the fractional part of x, according to the implementation's current rounding + mode. 2^23 is the smallest float that can be represented using all 23 significant + digits. */ TWO23[2]={ 8.3886080000e+06, /* 0x4b000000 */ -8.3886080000e+06, /* 0xcb000000 */ @@ -54,7 +58,7 @@ TWO23[2]={ /* Extract exponent field. */ j0 = ((i0 & 0x7f800000) >> 23) - 127; - if (j0 < (sizeof (long int) * 8) - 1) + if (j0 < (int)(sizeof (long int) * 8) - 1) { if (j0 < -1) return 0; @@ -65,6 +69,10 @@ TWO23[2]={ w = TWO23[sx] + x; t = w - TWO23[sx]; GET_FLOAT_WORD (i0, t); + /* Detect the all-zeros representation of plus and + minus zero, which fails the calculation below. */ + if ((i0 & ~(1 << 31)) == 0) + return 0; j0 = ((i0 >> 23) & 0xff) - 0x7f; i0 &= 0x7fffff; i0 |= 0x800000; |