summaryrefslogtreecommitdiff
path: root/stdlib/strtod.c
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>1999-03-08 21:41:47 +0000
committerUlrich Drepper <drepper@redhat.com>1999-03-08 21:41:47 +0000
commit5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715 (patch)
treef33e10763cba76eb5dd8c5c86db7f7b07b9cc1cc /stdlib/strtod.c
parent29329ef54739d84393b59058d83ba23e40ca88e9 (diff)
downloadglibc-5be4316ea2f4bb8a59a7861b5dd3c5ebfea6a715.tar.gz
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.c12
1 files changed, 9 insertions, 3 deletions
diff --git a/stdlib/strtod.c b/stdlib/strtod.c
index 5b15859066..0de96a571b 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)