diff options
| author | Ulrich Drepper <drepper@redhat.com> | 1999-03-08 21:41:47 +0000 |
|---|---|---|
| committer | Ulrich Drepper <drepper@redhat.com> | 1999-03-08 21:41:47 +0000 |
| commit | 5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715 (patch) | |
| tree | f33e10763cba76eb5dd8c5c86db7f7b07b9cc1cc /stdlib/strtod.c | |
| parent | 29329ef54739d84393b59058d83ba23e40ca88e9 (diff) | |
| download | glibc-5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715.zip glibc-5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715.tar.gz glibc-5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715.tar.bz2 | |
Update.
1999-03-09 Andreas Schwab <schwab@issan.cs.uni-dortmund.de>
* stdlib/strtod.c (round_and_return): Fix rounding of denormals.
(DENORM_EXP): New define. Use it as the exponent for a denormal.
* stdlib/strtold.c: Moved to ...
* sysdeps/generic/strtold.c: ... here.
* sysdeps/m68k/strtold.c: New file.
1999-03-09 Andreas Schwab <schwab@issan.cs.uni-dortmund.de>
* stdio-common/printf_fphex.c (__printf_fphex): Print no decimal
point if the precision is zero, even if the mantissa is non-zero.
Diffstat (limited to 'stdlib/strtod.c')
| -rw-r--r-- | stdlib/strtod.c | 12 |
1 files changed, 9 insertions, 3 deletions
diff --git a/stdlib/strtod.c b/stdlib/strtod.c index 5b15859..0de96a5 100644 --- a/stdlib/strtod.c +++ b/stdlib/strtod.c @@ -226,7 +226,7 @@ round_and_return (mp_limb_t *retval, int exponent, int negative, int i; round_limb = retval[RETURN_LIMB_SIZE - 1]; - round_bit = BITS_PER_MP_LIMB - 1; + round_bit = (MANT_DIG - 1) % BITS_PER_MP_LIMB; for (i = 0; i < RETURN_LIMB_SIZE; ++i) more_bits |= retval[i] != 0; MPN_ZERO (retval, RETURN_LIMB_SIZE); @@ -254,7 +254,13 @@ round_and_return (mp_limb_t *retval, int exponent, int negative, round_bit = shift - 1; (void) __mpn_rshift (retval, retval, RETURN_LIMB_SIZE, shift); } - exponent = MIN_EXP - 2; + /* This is a hook for the m68k long double format, where the + exponent bias is the same for normalized and denormalized + numbers. */ +#ifndef DENORM_EXP +# define DENORM_EXP (MIN_EXP - 2) +#endif + exponent = DENORM_EXP; } if ((round_limb & (((mp_limb_t) 1) << round_bit)) != 0 @@ -273,7 +279,7 @@ round_and_return (mp_limb_t *retval, int exponent, int negative, retval[RETURN_LIMB_SIZE - 1] |= ((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB); } - else if (exponent == MIN_EXP - 2 + else if (exponent == DENORM_EXP && (retval[RETURN_LIMB_SIZE - 1] & (((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB))) != 0) |