kc3-lang/gnulib/lib/float.c
Branch
-
Show log
Commit
-
Hash : 1ce8849f
Author :
Date : 2025-07-28T17:50:30
float-h: change IBM long double to match GCC 15 This is in response to the GCC developers; see <https://gcc.gnu.org/bugzilla/show_bug.cgi?id=120993>. * lib/float.c (gl_LDBL_MAX): Remove the implementation for PowerPC, as it should no longer be needed. * lib/float.in.h (LDBL_MIN_EXP, LDBL_MIN_10_EXP, LDBL_MIN) (LDBL_MAX, LDBL_EPSILON, LDBL_NORM_MAX): On PowerPC with IBM long double, simplify by using the GCC 15 values unconditionally. These are the correct values according to the GCC developers and there seems little point to disagreeing with current GCC about obsolescent arithmetic that is so problematic in practice. * tests/test-float-h.c (test_long_double): Relax tests of LDBL_MAX when !LDBL_IS_IEC_60559, as the tests would now fail on PowerPC and they were not portable in that case anyway.
-
lib/float.c
-
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117
/* Auxiliary definitions for <float.h>. Copyright (C) 2011-2025 Free Software Foundation, Inc. Written by Bruno Haible <bruno@clisp.org>, 2011. This file 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. This file 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 this program. If not, see <https://www.gnu.org/licenses/>. */ #include <config.h> /* Specification. */ #include <float.h> #if GNULIB_defined_long_double_union # ifdef __i386__ const union gl_long_double_union gl_LDBL_MAX = { { 0xFFFFFFFF, 0xFFFFFFFF, 32766 } }; # endif # if defined __i386__ && (defined __FreeBSD__ || defined __DragonFly__) /* We can't even simply evaluate the formula (LDBL_MIN / 9223372036854775808.0L) at run time, because it would require BEGIN_LONG_DOUBLE_ROUNDING / END_LONG_DOUBLE_ROUNDING invocations. It simpler to just write down the representation of LDBL_TRUE_MIN, based on <https://en.wikipedia.org/wiki/Extended_precision#x86_extended_precision_format>. */ const union gl_long_double_union gl_LDBL_TRUE_MIN = { { 0x00000001, 0x00000000, 0 } }; # endif #endif #if GNULIB_defined_FLT_SNAN /* Define like memory_positive_SNaNf(), see signed-snan.h and snan.h, or like setpayloadsigf() with an arbitrary payload. */ gl_FLT_SNAN_t gl_FLT_SNAN = # if FLT_MANT_DIG == 24 # if defined __hppa || (defined __mips__ && !MIPS_NAN2008_FLOAT) || defined __sh__ /* sign bit: 0, 8 exponent bits: all 1, next bit: 1, payload: 0b10...0 */ { { 0x7FE00000U } } # else /* sign bit: 0, 8 exponent bits: all 1, next bit: 0, payload: 0b10...0 */ { { 0x7FA00000U } } # endif # endif ; #endif #if GNULIB_defined_DBL_SNAN /* Define like memory_positive_SNaNd(), see signed-snan.h and snan.h, or like setpayloadsig() with an arbitrary payload. */ gl_DBL_SNAN_t gl_DBL_SNAN = # if DBL_MANT_DIG == 53 # if defined __hppa || (defined __mips__ && !MIPS_NAN2008_FLOAT) || defined __sh__ /* sign bit: 0, 11 exponent bits: all 1, next bit: 1, payload: 0b10...0 */ { { 0x7FFC000000000000ULL } } # else /* sign bit: 0, 11 exponent bits: all 1, next bit: 0, payload: 0b10...0 */ { { 0x7FF4000000000000ULL } } # endif # endif ; #endif #if GNULIB_defined_LDBL_SNAN # ifdef WORDS_BIGENDIAN # define TWO(hi,lo) { hi, lo } # else # define TWO(hi,lo) { lo, hi } # endif /* Define like memory_positive_SNaNl(), see signed-snan.h and snan.h, or like setpayloadsigl() with an arbitrary payload. */ gl_LDBL_SNAN_t gl_LDBL_SNAN = # if LDBL_MANT_DIG == 53 /* on arm, hppa, mips, sh, but also MSVC */ # if defined __hppa || (defined __mips__ && !MIPS_NAN2008_FLOAT) || defined __sh__ /* sign bit: 0, 11 exponent bits: all 1, next bit: 1, payload: 0b10...0 */ { { 0x7FFC000000000000ULL } } # else /* sign bit: 0, 11 exponent bits: all 1, next bit: 0, payload: 0b10...0 */ { { 0x7FF4000000000000ULL } } # endif # elif LDBL_MANT_DIG == 64 /* on i386, x86_64, ia64, m68k */ # if defined __m68k__ /* sign bit: 0, 15 exponent bits: all 1, 16 gap bits: all 0, always=1 bit: 1, next bit: 0, payload: 0b10...0 */ { { 0x7FFF0000ULL, 0xA0000000ULL, 0x00000000ULL } } # else /* sign bit: 0, 15 exponent bits: all 1, always=1 bit: 1, next bit: 0, payload: 0b10...0 (see <https://en.wikipedia.org/wiki/Extended_precision#x86_extended_precision_format>) */ { TWO (0x00007FFFULL, 0xA000000000000000ULL) } # endif # elif LDBL_MANT_DIG == 106 /* on powerpc, powerpc64, powerpc64le */ /* most-significant double: sign bit: 0, 11 exponent bits: all 1, next bit: 0, payload: 0b10...0, least-significant double: 0.0 */ { { 0x7FF4000000000000ULL, 0ULL } } # elif LDBL_MANT_DIG == 113 /* on alpha, arm64, loongarch64, mips64, riscv64, s390x, sparc64 */ # if (defined __mips__ && !MIPS_NAN2008_FLOAT) /* sign bit: 0, 15 exponent bits: all 1, next bit: 1, payload: 0b10...0 */ { TWO (0x7FFFC00000000000ULL, 0ULL) } # else /* sign bit: 0, 15 exponent bits: all 1, next bit: 0, payload: 0b10...0 */ { TWO (0x7FFF400000000000ULL, 0ULL) } # endif # endif ; #endif /* This declaration is solely to ensure that after preprocessing this file is never empty. */ typedef int dummy;