1 /* Print floating point number in hexadecimal notation according to ISO C99. 2 Copyright (C) 1997-2022 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 4 5 The GNU C Library is free software; you can redistribute it and/or 6 modify it under the terms of the GNU Lesser General Public 7 License as published by the Free Software Foundation; either 8 version 2.1 of the License, or (at your option) any later version. 9 10 The GNU C Library is distributed in the hope that it will be useful, 11 but WITHOUT ANY WARRANTY; without even the implied warranty of 12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 13 Lesser General Public License for more details. 14 15 You should have received a copy of the GNU Lesser General Public 16 License along with the GNU C Library; if not, see 17 <https://www.gnu.org/licenses/>. */ 18 19 #define PRINT_FPHEX_LONG_DOUBLE \ 20 do { \ 21 /* We have 105 bits of mantissa plus one implicit digit. Since \ 22 106 bits are representable without rest using hexadecimal \ 23 digits we use only the implicit digits for the number before \ 24 the decimal point. */ \ 25 unsigned long long int num0, num1; \ 26 unsigned long long hi, lo; \ 27 int ediff; \ 28 union ibm_extended_long_double u; \ 29 u.ld = fpnum.ldbl; \ 30 \ 31 assert (sizeof (long double) == 16); \ 32 \ 33 lo = ((long long)u.d[1].ieee.mantissa0 << 32) | u.d[1].ieee.mantissa1; \ 34 hi = ((long long)u.d[0].ieee.mantissa0 << 32) | u.d[0].ieee.mantissa1; \ 35 lo <<= 7; /* pre-shift lo to match ieee854. */ \ 36 /* If the lower double is not a denormal or zero then set the hidden \ 37 53rd bit. */ \ 38 if (u.d[1].ieee.exponent != 0) \ 39 lo |= (1ULL << (52 + 7)); \ 40 else \ 41 lo <<= 1; \ 42 /* The lower double is normalized separately from the upper. We \ 43 may need to adjust the lower manitissa to reflect this. */ \ 44 ediff = u.d[0].ieee.exponent - u.d[1].ieee.exponent - 53; \ 45 if (ediff > 63) \ 46 lo = 0; \ 47 else if (ediff > 0) \ 48 lo = lo >> ediff; \ 49 else if (ediff < 0) \ 50 lo = lo << -ediff; \ 51 if (u.d[0].ieee.negative != u.d[1].ieee.negative \ 52 && lo != 0) \ 53 { \ 54 lo = (1ULL << 60) - lo; \ 55 if (hi == 0L) \ 56 { \ 57 /* we have a borrow from the hidden bit, so shift left 1. */ \ 58 hi = 0xffffffffffffeLL | (lo >> 59); \ 59 lo = 0xfffffffffffffffLL & (lo << 1); \ 60 u.d[0].ieee.exponent--; \ 61 } \ 62 else \ 63 hi--; \ 64 } \ 65 num1 = (hi << 60) | lo; \ 66 num0 = hi >> 4; \ 67 \ 68 zero_mantissa = (num0|num1) == 0; \ 69 \ 70 if (sizeof (unsigned long int) > 6) \ 71 { \ 72 numstr = _itoa_word (num1, numbuf + sizeof numbuf, 16, \ 73 info->spec == 'A'); \ 74 wnumstr = _itowa_word (num1, \ 75 wnumbuf + sizeof (wnumbuf) / sizeof (wchar_t),\ 76 16, info->spec == 'A'); \ 77 } \ 78 else \ 79 { \ 80 numstr = _itoa (num1, numbuf + sizeof numbuf, 16, \ 81 info->spec == 'A'); \ 82 wnumstr = _itowa (num1, \ 83 wnumbuf + sizeof (wnumbuf) / sizeof (wchar_t), \ 84 16, info->spec == 'A'); \ 85 } \ 86 \ 87 while (numstr > numbuf + (sizeof numbuf - 64 / 4)) \ 88 { \ 89 *--numstr = '0'; \ 90 *--wnumstr = L'0'; \ 91 } \ 92 \ 93 if (sizeof (unsigned long int) > 6) \ 94 { \ 95 numstr = _itoa_word (num0, numstr, 16, info->spec == 'A'); \ 96 wnumstr = _itowa_word (num0, wnumstr, 16, info->spec == 'A'); \ 97 } \ 98 else \ 99 { \ 100 numstr = _itoa (num0, numstr, 16, info->spec == 'A'); \ 101 wnumstr = _itowa (num0, wnumstr, 16, info->spec == 'A'); \ 102 } \ 103 \ 104 /* Fill with zeroes. */ \ 105 while (numstr > numbuf + (sizeof numbuf - 112 / 4)) \ 106 { \ 107 *--numstr = '0'; \ 108 *--wnumstr = L'0'; \ 109 } \ 110 \ 111 leading = u.d[0].ieee.exponent == 0 ? '0' : '1'; \ 112 \ 113 exponent = u.d[0].ieee.exponent; \ 114 \ 115 if (exponent == 0) \ 116 { \ 117 if (zero_mantissa) \ 118 expnegative = 0; \ 119 else \ 120 { \ 121 /* This is a denormalized number. */ \ 122 expnegative = 1; \ 123 exponent = IEEE754_DOUBLE_BIAS - 1; \ 124 } \ 125 } \ 126 else if (exponent >= IEEE754_DOUBLE_BIAS) \ 127 { \ 128 expnegative = 0; \ 129 exponent -= IEEE754_DOUBLE_BIAS; \ 130 } \ 131 else \ 132 { \ 133 expnegative = 1; \ 134 exponent = -(exponent - IEEE754_DOUBLE_BIAS); \ 135 } \ 136 } while (0) 137 138 #include <stdio-common/printf_fphex.c> 139