mirror of
https://sourceware.org/git/glibc.git
synced 2024-12-04 19:00:09 +00:00
1b6adf888d
http://sourceware.org/ml/libc-alpha/2013-08/msg00081.html This is the first of a series of patches to ban ieee854_long_double and the ieee854_long_double macros when using IBM long double. union ieee854_long_double just isn't correct for IBM long double, especially when little-endian, and pretending it is OK has allowed a number of bugs to remain undetected in sysdeps/ieee754/ldbl-128ibm/. This changes the few places in generic code that use it. * stdio-common/printf_size.c (__printf_size): Don't use union ieee854_long_double in fpnum union. * stdio-common/printf_fphex.c (__printf_fphex): Likewise. Use signbit macro to retrieve sign from long double. * stdio-common/printf_fp.c (___printf_fp): Use signbit macro to retrieve sign from long double. * sysdeps/ieee754/ldbl-128ibm/printf_fphex.c: Adjust for fpnum change. * sysdeps/ieee754/ldbl-128/printf_fphex.c: Likewise. * sysdeps/ieee754/ldbl-96/printf_fphex.c: Likewise. * sysdeps/x86_64/fpu/printf_fphex.c: Likewise. * math/test-misc.c (main): Don't use union ieee854_long_double. ports/ * sysdeps/ia64/fpu/printf_fphex.c: Adjust for fpnum change.
94 lines
3.6 KiB
C
94 lines
3.6 KiB
C
/* Print floating point number in hexadecimal notation according to ISO C99.
|
|
Copyright (C) 1997-2013 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
|
|
modify it under the terms of the GNU Lesser General Public
|
|
License as published by the Free Software Foundation; either
|
|
version 2.1 of the License, or (at your option) any later version.
|
|
|
|
The GNU C Library is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Lesser General Public License for more details.
|
|
|
|
You should have received a copy of the GNU Lesser General Public
|
|
License along with the GNU C Library; if not, see
|
|
<http://www.gnu.org/licenses/>. */
|
|
|
|
#ifndef LONG_DOUBLE_DENORM_BIAS
|
|
# define LONG_DOUBLE_DENORM_BIAS (IEEE854_LONG_DOUBLE_BIAS - 1)
|
|
#endif
|
|
|
|
#define PRINT_FPHEX_LONG_DOUBLE \
|
|
do { \
|
|
/* The "strange" 80 bit format on ix86 and m68k has an explicit \
|
|
leading digit in the 64 bit mantissa. */ \
|
|
unsigned long long int num; \
|
|
union ieee854_long_double u; \
|
|
u.d = fpnum.ldbl; \
|
|
\
|
|
num = (((unsigned long long int) u.ieee.mantissa0) << 32 \
|
|
| u.ieee.mantissa1); \
|
|
\
|
|
zero_mantissa = num == 0; \
|
|
\
|
|
if (sizeof (unsigned long int) > 6) \
|
|
{ \
|
|
numstr = _itoa_word (num, numbuf + sizeof numbuf, 16, \
|
|
info->spec == 'A'); \
|
|
wnumstr = _itowa_word (num, \
|
|
wnumbuf + sizeof (wnumbuf) / sizeof (wchar_t),\
|
|
16, info->spec == 'A'); \
|
|
} \
|
|
else \
|
|
{ \
|
|
numstr = _itoa (num, numbuf + sizeof numbuf, 16, info->spec == 'A');\
|
|
wnumstr = _itowa (num, \
|
|
wnumbuf + sizeof (wnumbuf) / sizeof (wchar_t), \
|
|
16, info->spec == 'A'); \
|
|
} \
|
|
\
|
|
/* Fill with zeroes. */ \
|
|
while (numstr > numbuf + (sizeof numbuf - 64 / 4)) \
|
|
{ \
|
|
*--numstr = '0'; \
|
|
*--wnumstr = L'0'; \
|
|
} \
|
|
\
|
|
/* We use a full nibble for the leading digit. */ \
|
|
leading = *numstr++; \
|
|
wnumstr++; \
|
|
\
|
|
/* We have 3 bits from the mantissa in the leading nibble. \
|
|
Therefore we are here using `IEEE854_LONG_DOUBLE_BIAS + 3'. */ \
|
|
exponent = u.ieee.exponent; \
|
|
\
|
|
if (exponent == 0) \
|
|
{ \
|
|
if (zero_mantissa) \
|
|
expnegative = 0; \
|
|
else \
|
|
{ \
|
|
/* This is a denormalized number. */ \
|
|
expnegative = 1; \
|
|
/* This is a hook for the m68k long double format, where the \
|
|
exponent bias is the same for normalized and denormalized \
|
|
numbers. */ \
|
|
exponent = LONG_DOUBLE_DENORM_BIAS + 3; \
|
|
} \
|
|
} \
|
|
else if (exponent >= IEEE854_LONG_DOUBLE_BIAS + 3) \
|
|
{ \
|
|
expnegative = 0; \
|
|
exponent -= IEEE854_LONG_DOUBLE_BIAS + 3; \
|
|
} \
|
|
else \
|
|
{ \
|
|
expnegative = 1; \
|
|
exponent = -(exponent - (IEEE854_LONG_DOUBLE_BIAS + 3)); \
|
|
} \
|
|
} while (0)
|
|
|
|
#include <stdio-common/printf_fphex.c>
|