diff options
| author | Andreas Schwab <schwab@linux-m68k.org> | 2012-04-03 18:38:46 +0200 |
|---|---|---|
| committer | Andreas Schwab <schwab@linux-m68k.org> | 2012-04-28 22:21:27 +0200 |
| commit | 7e0d315da842379f077fc61dca72460fb0e9263e (patch) | |
| tree | 252f9f59a4b17e807b33979c51ecdff86a07720c /sysdeps/ieee754 | |
| parent | 6d5c57fabd1a97b110b00e59fd52f50ce5bbb1a3 (diff) | |
| download | glibc-7e0d315da842379f077fc61dca72460fb0e9263e.tar.gz glibc-7e0d315da842379f077fc61dca72460fb0e9263e.tar.xz glibc-7e0d315da842379f077fc61dca72460fb0e9263e.zip | |
Fix formatting of denormal IBM long double numbers
Diffstat (limited to 'sysdeps/ieee754')
| -rw-r--r-- | sysdeps/ieee754/ldbl-128ibm/ldbl2mpn.c | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/ldbl2mpn.c b/sysdeps/ieee754/ldbl-128ibm/ldbl2mpn.c index 3162bbd9a5..83162c5b32 100644 --- a/sysdeps/ieee754/ldbl-128ibm/ldbl2mpn.c +++ b/sysdeps/ieee754/ldbl-128ibm/ldbl2mpn.c @@ -1,5 +1,4 @@ -/* Copyright (C) 1995,1996,1997,1998,1999,2002,2003,2006 - Free Software Foundation, Inc. +/* Copyright (C) 1995-2012 Free Software Foundation, Inc. This file is part of the GNU C Library. The GNU C Library is free software; you can redistribute it and/or @@ -104,7 +103,10 @@ __mpn_extract_long_double (mp_ptr res_ptr, mp_size_t size, else { /* It is a denormal number, meaning it has no implicit leading - one bit, and its exponent is in fact the format minimum. */ + one bit, and its exponent is in fact the format minimum. We + use DBL_MIN_EXP instead of LDBL_MIN_EXP below because the + latter describes the properties of both parts together, but + the exponent is computed from the high part only. */ int cnt; #if N == 2 @@ -115,7 +117,7 @@ __mpn_extract_long_double (mp_ptr res_ptr, mp_size_t size, res_ptr[N - 1] = res_ptr[N - 1] << cnt | (res_ptr[0] >> (BITS_PER_MP_LIMB - cnt)); res_ptr[0] <<= cnt; - *expt = LDBL_MIN_EXP - 1 - cnt; + *expt = DBL_MIN_EXP - 1 - cnt; } else { @@ -130,7 +132,7 @@ __mpn_extract_long_double (mp_ptr res_ptr, mp_size_t size, res_ptr[N - 1] = res_ptr[0] >> (NUM_LEADING_ZEROS - cnt); res_ptr[0] <<= BITS_PER_MP_LIMB - (NUM_LEADING_ZEROS - cnt); } - *expt = LDBL_MIN_EXP - 1 + *expt = DBL_MIN_EXP - 1 - (BITS_PER_MP_LIMB - NUM_LEADING_ZEROS) - cnt; } #else @@ -161,7 +163,7 @@ __mpn_extract_long_double (mp_ptr res_ptr, mp_size_t size, for (; k >= 0; k--) res_ptr[k] = 0; - *expt = LDBL_MIN_EXP - 1 - l * BITS_PER_MP_LIMB - cnt; + *expt = DBL_MIN_EXP - 1 - l * BITS_PER_MP_LIMB - cnt; #endif } } |