1*820c1a8dSHiroo HAYASHI /* intprops.h -- properties of integer types 2*820c1a8dSHiroo HAYASHI 3*820c1a8dSHiroo HAYASHI Copyright (C) 2001-2021 Free Software Foundation, Inc. 4*820c1a8dSHiroo HAYASHI 5*820c1a8dSHiroo HAYASHI This program is free software: you can redistribute it and/or modify it 6*820c1a8dSHiroo HAYASHI under the terms of the GNU Lesser General Public License as published 7*820c1a8dSHiroo HAYASHI by the Free Software Foundation; either version 2.1 of the License, or 8*820c1a8dSHiroo HAYASHI (at your option) any later version. 9*820c1a8dSHiroo HAYASHI 10*820c1a8dSHiroo HAYASHI This program is distributed in the hope that it will be useful, 11*820c1a8dSHiroo HAYASHI but WITHOUT ANY WARRANTY; without even the implied warranty of 12*820c1a8dSHiroo HAYASHI MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 13*820c1a8dSHiroo HAYASHI GNU Lesser General Public License for more details. 14*820c1a8dSHiroo HAYASHI 15*820c1a8dSHiroo HAYASHI You should have received a copy of the GNU Lesser General Public License 16*820c1a8dSHiroo HAYASHI along with this program. If not, see <https://www.gnu.org/licenses/>. */ 17*820c1a8dSHiroo HAYASHI 18*820c1a8dSHiroo HAYASHI /* Written by Paul Eggert. */ 19*820c1a8dSHiroo HAYASHI 20*820c1a8dSHiroo HAYASHI #ifndef _GL_INTPROPS_H 21*820c1a8dSHiroo HAYASHI #define _GL_INTPROPS_H 22*820c1a8dSHiroo HAYASHI 23*820c1a8dSHiroo HAYASHI #include <limits.h> 24*820c1a8dSHiroo HAYASHI 25*820c1a8dSHiroo HAYASHI /* Return a value with the common real type of E and V and the value of V. 26*820c1a8dSHiroo HAYASHI Do not evaluate E. */ 27*820c1a8dSHiroo HAYASHI #define _GL_INT_CONVERT(e, v) ((1 ? 0 : (e)) + (v)) 28*820c1a8dSHiroo HAYASHI 29*820c1a8dSHiroo HAYASHI /* Act like _GL_INT_CONVERT (E, -V) but work around a bug in IRIX 6.5 cc; see 30*820c1a8dSHiroo HAYASHI <https://lists.gnu.org/r/bug-gnulib/2011-05/msg00406.html>. */ 31*820c1a8dSHiroo HAYASHI #define _GL_INT_NEGATE_CONVERT(e, v) ((1 ? 0 : (e)) - (v)) 32*820c1a8dSHiroo HAYASHI 33*820c1a8dSHiroo HAYASHI /* The extra casts in the following macros work around compiler bugs, 34*820c1a8dSHiroo HAYASHI e.g., in Cray C 5.0.3.0. */ 35*820c1a8dSHiroo HAYASHI 36*820c1a8dSHiroo HAYASHI /* True if the arithmetic type T is an integer type. bool counts as 37*820c1a8dSHiroo HAYASHI an integer. */ 38*820c1a8dSHiroo HAYASHI #define TYPE_IS_INTEGER(t) ((t) 1.5 == 1) 39*820c1a8dSHiroo HAYASHI 40*820c1a8dSHiroo HAYASHI /* True if the real type T is signed. */ 41*820c1a8dSHiroo HAYASHI #define TYPE_SIGNED(t) (! ((t) 0 < (t) -1)) 42*820c1a8dSHiroo HAYASHI 43*820c1a8dSHiroo HAYASHI /* Return 1 if the real expression E, after promotion, has a 44*820c1a8dSHiroo HAYASHI signed or floating type. Do not evaluate E. */ 45*820c1a8dSHiroo HAYASHI #define EXPR_SIGNED(e) (_GL_INT_NEGATE_CONVERT (e, 1) < 0) 46*820c1a8dSHiroo HAYASHI 47*820c1a8dSHiroo HAYASHI 48*820c1a8dSHiroo HAYASHI /* Minimum and maximum values for integer types and expressions. */ 49*820c1a8dSHiroo HAYASHI 50*820c1a8dSHiroo HAYASHI /* The width in bits of the integer type or expression T. 51*820c1a8dSHiroo HAYASHI Do not evaluate T. T must not be a bit-field expression. 52*820c1a8dSHiroo HAYASHI Padding bits are not supported; this is checked at compile-time below. */ 53*820c1a8dSHiroo HAYASHI #define TYPE_WIDTH(t) (sizeof (t) * CHAR_BIT) 54*820c1a8dSHiroo HAYASHI 55*820c1a8dSHiroo HAYASHI /* The maximum and minimum values for the integer type T. */ 56*820c1a8dSHiroo HAYASHI #define TYPE_MINIMUM(t) ((t) ~ TYPE_MAXIMUM (t)) 57*820c1a8dSHiroo HAYASHI #define TYPE_MAXIMUM(t) \ 58*820c1a8dSHiroo HAYASHI ((t) (! TYPE_SIGNED (t) \ 59*820c1a8dSHiroo HAYASHI ? (t) -1 \ 60*820c1a8dSHiroo HAYASHI : ((((t) 1 << (TYPE_WIDTH (t) - 2)) - 1) * 2 + 1))) 61*820c1a8dSHiroo HAYASHI 62*820c1a8dSHiroo HAYASHI /* The maximum and minimum values for the type of the expression E, 63*820c1a8dSHiroo HAYASHI after integer promotion. E is not evaluated. */ 64*820c1a8dSHiroo HAYASHI #define _GL_INT_MINIMUM(e) \ 65*820c1a8dSHiroo HAYASHI (EXPR_SIGNED (e) \ 66*820c1a8dSHiroo HAYASHI ? ~ _GL_SIGNED_INT_MAXIMUM (e) \ 67*820c1a8dSHiroo HAYASHI : _GL_INT_CONVERT (e, 0)) 68*820c1a8dSHiroo HAYASHI #define _GL_INT_MAXIMUM(e) \ 69*820c1a8dSHiroo HAYASHI (EXPR_SIGNED (e) \ 70*820c1a8dSHiroo HAYASHI ? _GL_SIGNED_INT_MAXIMUM (e) \ 71*820c1a8dSHiroo HAYASHI : _GL_INT_NEGATE_CONVERT (e, 1)) 72*820c1a8dSHiroo HAYASHI #define _GL_SIGNED_INT_MAXIMUM(e) \ 73*820c1a8dSHiroo HAYASHI (((_GL_INT_CONVERT (e, 1) << (TYPE_WIDTH (+ (e)) - 2)) - 1) * 2 + 1) 74*820c1a8dSHiroo HAYASHI 75*820c1a8dSHiroo HAYASHI /* Work around OpenVMS incompatibility with C99. */ 76*820c1a8dSHiroo HAYASHI #if !defined LLONG_MAX && defined __INT64_MAX 77*820c1a8dSHiroo HAYASHI # define LLONG_MAX __INT64_MAX 78*820c1a8dSHiroo HAYASHI # define LLONG_MIN __INT64_MIN 79*820c1a8dSHiroo HAYASHI #endif 80*820c1a8dSHiroo HAYASHI 81*820c1a8dSHiroo HAYASHI /* This include file assumes that signed types are two's complement without 82*820c1a8dSHiroo HAYASHI padding bits; the above macros have undefined behavior otherwise. 83*820c1a8dSHiroo HAYASHI If this is a problem for you, please let us know how to fix it for your host. 84*820c1a8dSHiroo HAYASHI This assumption is tested by the intprops-tests module. */ 85*820c1a8dSHiroo HAYASHI 86*820c1a8dSHiroo HAYASHI /* Does the __typeof__ keyword work? This could be done by 87*820c1a8dSHiroo HAYASHI 'configure', but for now it's easier to do it by hand. */ 88*820c1a8dSHiroo HAYASHI #if (2 <= __GNUC__ \ 89*820c1a8dSHiroo HAYASHI || (4 <= __clang_major__) \ 90*820c1a8dSHiroo HAYASHI || (1210 <= __IBMC__ && defined __IBM__TYPEOF__) \ 91*820c1a8dSHiroo HAYASHI || (0x5110 <= __SUNPRO_C && !__STDC__)) 92*820c1a8dSHiroo HAYASHI # define _GL_HAVE___TYPEOF__ 1 93*820c1a8dSHiroo HAYASHI #else 94*820c1a8dSHiroo HAYASHI # define _GL_HAVE___TYPEOF__ 0 95*820c1a8dSHiroo HAYASHI #endif 96*820c1a8dSHiroo HAYASHI 97*820c1a8dSHiroo HAYASHI /* Return 1 if the integer type or expression T might be signed. Return 0 98*820c1a8dSHiroo HAYASHI if it is definitely unsigned. T must not be a bit-field expression. 99*820c1a8dSHiroo HAYASHI This macro does not evaluate its argument, and expands to an 100*820c1a8dSHiroo HAYASHI integer constant expression. */ 101*820c1a8dSHiroo HAYASHI #if _GL_HAVE___TYPEOF__ 102*820c1a8dSHiroo HAYASHI # define _GL_SIGNED_TYPE_OR_EXPR(t) TYPE_SIGNED (__typeof__ (t)) 103*820c1a8dSHiroo HAYASHI #else 104*820c1a8dSHiroo HAYASHI # define _GL_SIGNED_TYPE_OR_EXPR(t) 1 105*820c1a8dSHiroo HAYASHI #endif 106*820c1a8dSHiroo HAYASHI 107*820c1a8dSHiroo HAYASHI /* Bound on length of the string representing an unsigned integer 108*820c1a8dSHiroo HAYASHI value representable in B bits. log10 (2.0) < 146/485. The 109*820c1a8dSHiroo HAYASHI smallest value of B where this bound is not tight is 2621. */ 110*820c1a8dSHiroo HAYASHI #define INT_BITS_STRLEN_BOUND(b) (((b) * 146 + 484) / 485) 111*820c1a8dSHiroo HAYASHI 112*820c1a8dSHiroo HAYASHI /* Bound on length of the string representing an integer type or expression T. 113*820c1a8dSHiroo HAYASHI T must not be a bit-field expression. 114*820c1a8dSHiroo HAYASHI 115*820c1a8dSHiroo HAYASHI Subtract 1 for the sign bit if T is signed, and then add 1 more for 116*820c1a8dSHiroo HAYASHI a minus sign if needed. 117*820c1a8dSHiroo HAYASHI 118*820c1a8dSHiroo HAYASHI Because _GL_SIGNED_TYPE_OR_EXPR sometimes returns 1 when its argument is 119*820c1a8dSHiroo HAYASHI unsigned, this macro may overestimate the true bound by one byte when 120*820c1a8dSHiroo HAYASHI applied to unsigned types of size 2, 4, 16, ... bytes. */ 121*820c1a8dSHiroo HAYASHI #define INT_STRLEN_BOUND(t) \ 122*820c1a8dSHiroo HAYASHI (INT_BITS_STRLEN_BOUND (TYPE_WIDTH (t) - _GL_SIGNED_TYPE_OR_EXPR (t)) \ 123*820c1a8dSHiroo HAYASHI + _GL_SIGNED_TYPE_OR_EXPR (t)) 124*820c1a8dSHiroo HAYASHI 125*820c1a8dSHiroo HAYASHI /* Bound on buffer size needed to represent an integer type or expression T, 126*820c1a8dSHiroo HAYASHI including the terminating null. T must not be a bit-field expression. */ 127*820c1a8dSHiroo HAYASHI #define INT_BUFSIZE_BOUND(t) (INT_STRLEN_BOUND (t) + 1) 128*820c1a8dSHiroo HAYASHI 129*820c1a8dSHiroo HAYASHI 130*820c1a8dSHiroo HAYASHI /* Range overflow checks. 131*820c1a8dSHiroo HAYASHI 132*820c1a8dSHiroo HAYASHI The INT_<op>_RANGE_OVERFLOW macros return 1 if the corresponding C 133*820c1a8dSHiroo HAYASHI operators might not yield numerically correct answers due to 134*820c1a8dSHiroo HAYASHI arithmetic overflow. They do not rely on undefined or 135*820c1a8dSHiroo HAYASHI implementation-defined behavior. Their implementations are simple 136*820c1a8dSHiroo HAYASHI and straightforward, but they are harder to use and may be less 137*820c1a8dSHiroo HAYASHI efficient than the INT_<op>_WRAPV, INT_<op>_OK, and 138*820c1a8dSHiroo HAYASHI INT_<op>_OVERFLOW macros described below. 139*820c1a8dSHiroo HAYASHI 140*820c1a8dSHiroo HAYASHI Example usage: 141*820c1a8dSHiroo HAYASHI 142*820c1a8dSHiroo HAYASHI long int i = ...; 143*820c1a8dSHiroo HAYASHI long int j = ...; 144*820c1a8dSHiroo HAYASHI if (INT_MULTIPLY_RANGE_OVERFLOW (i, j, LONG_MIN, LONG_MAX)) 145*820c1a8dSHiroo HAYASHI printf ("multiply would overflow"); 146*820c1a8dSHiroo HAYASHI else 147*820c1a8dSHiroo HAYASHI printf ("product is %ld", i * j); 148*820c1a8dSHiroo HAYASHI 149*820c1a8dSHiroo HAYASHI Restrictions on *_RANGE_OVERFLOW macros: 150*820c1a8dSHiroo HAYASHI 151*820c1a8dSHiroo HAYASHI These macros do not check for all possible numerical problems or 152*820c1a8dSHiroo HAYASHI undefined or unspecified behavior: they do not check for division 153*820c1a8dSHiroo HAYASHI by zero, for bad shift counts, or for shifting negative numbers. 154*820c1a8dSHiroo HAYASHI 155*820c1a8dSHiroo HAYASHI These macros may evaluate their arguments zero or multiple times, 156*820c1a8dSHiroo HAYASHI so the arguments should not have side effects. The arithmetic 157*820c1a8dSHiroo HAYASHI arguments (including the MIN and MAX arguments) must be of the same 158*820c1a8dSHiroo HAYASHI integer type after the usual arithmetic conversions, and the type 159*820c1a8dSHiroo HAYASHI must have minimum value MIN and maximum MAX. Unsigned types should 160*820c1a8dSHiroo HAYASHI use a zero MIN of the proper type. 161*820c1a8dSHiroo HAYASHI 162*820c1a8dSHiroo HAYASHI Because all arguments are subject to integer promotions, these 163*820c1a8dSHiroo HAYASHI macros typically do not work on types narrower than 'int'. 164*820c1a8dSHiroo HAYASHI 165*820c1a8dSHiroo HAYASHI These macros are tuned for constant MIN and MAX. For commutative 166*820c1a8dSHiroo HAYASHI operations such as A + B, they are also tuned for constant B. */ 167*820c1a8dSHiroo HAYASHI 168*820c1a8dSHiroo HAYASHI /* Return 1 if A + B would overflow in [MIN,MAX] arithmetic. 169*820c1a8dSHiroo HAYASHI See above for restrictions. */ 170*820c1a8dSHiroo HAYASHI #define INT_ADD_RANGE_OVERFLOW(a, b, min, max) \ 171*820c1a8dSHiroo HAYASHI ((b) < 0 \ 172*820c1a8dSHiroo HAYASHI ? (a) < (min) - (b) \ 173*820c1a8dSHiroo HAYASHI : (max) - (b) < (a)) 174*820c1a8dSHiroo HAYASHI 175*820c1a8dSHiroo HAYASHI /* Return 1 if A - B would overflow in [MIN,MAX] arithmetic. 176*820c1a8dSHiroo HAYASHI See above for restrictions. */ 177*820c1a8dSHiroo HAYASHI #define INT_SUBTRACT_RANGE_OVERFLOW(a, b, min, max) \ 178*820c1a8dSHiroo HAYASHI ((b) < 0 \ 179*820c1a8dSHiroo HAYASHI ? (max) + (b) < (a) \ 180*820c1a8dSHiroo HAYASHI : (a) < (min) + (b)) 181*820c1a8dSHiroo HAYASHI 182*820c1a8dSHiroo HAYASHI /* Return 1 if - A would overflow in [MIN,MAX] arithmetic. 183*820c1a8dSHiroo HAYASHI See above for restrictions. */ 184*820c1a8dSHiroo HAYASHI #define INT_NEGATE_RANGE_OVERFLOW(a, min, max) \ 185*820c1a8dSHiroo HAYASHI ((min) < 0 \ 186*820c1a8dSHiroo HAYASHI ? (a) < - (max) \ 187*820c1a8dSHiroo HAYASHI : 0 < (a)) 188*820c1a8dSHiroo HAYASHI 189*820c1a8dSHiroo HAYASHI /* Return 1 if A * B would overflow in [MIN,MAX] arithmetic. 190*820c1a8dSHiroo HAYASHI See above for restrictions. Avoid && and || as they tickle 191*820c1a8dSHiroo HAYASHI bugs in Sun C 5.11 2010/08/13 and other compilers; see 192*820c1a8dSHiroo HAYASHI <https://lists.gnu.org/r/bug-gnulib/2011-05/msg00401.html>. */ 193*820c1a8dSHiroo HAYASHI #define INT_MULTIPLY_RANGE_OVERFLOW(a, b, min, max) \ 194*820c1a8dSHiroo HAYASHI ((b) < 0 \ 195*820c1a8dSHiroo HAYASHI ? ((a) < 0 \ 196*820c1a8dSHiroo HAYASHI ? (a) < (max) / (b) \ 197*820c1a8dSHiroo HAYASHI : (b) == -1 \ 198*820c1a8dSHiroo HAYASHI ? 0 \ 199*820c1a8dSHiroo HAYASHI : (min) / (b) < (a)) \ 200*820c1a8dSHiroo HAYASHI : (b) == 0 \ 201*820c1a8dSHiroo HAYASHI ? 0 \ 202*820c1a8dSHiroo HAYASHI : ((a) < 0 \ 203*820c1a8dSHiroo HAYASHI ? (a) < (min) / (b) \ 204*820c1a8dSHiroo HAYASHI : (max) / (b) < (a))) 205*820c1a8dSHiroo HAYASHI 206*820c1a8dSHiroo HAYASHI /* Return 1 if A / B would overflow in [MIN,MAX] arithmetic. 207*820c1a8dSHiroo HAYASHI See above for restrictions. Do not check for division by zero. */ 208*820c1a8dSHiroo HAYASHI #define INT_DIVIDE_RANGE_OVERFLOW(a, b, min, max) \ 209*820c1a8dSHiroo HAYASHI ((min) < 0 && (b) == -1 && (a) < - (max)) 210*820c1a8dSHiroo HAYASHI 211*820c1a8dSHiroo HAYASHI /* Return 1 if A % B would overflow in [MIN,MAX] arithmetic. 212*820c1a8dSHiroo HAYASHI See above for restrictions. Do not check for division by zero. 213*820c1a8dSHiroo HAYASHI Mathematically, % should never overflow, but on x86-like hosts 214*820c1a8dSHiroo HAYASHI INT_MIN % -1 traps, and the C standard permits this, so treat this 215*820c1a8dSHiroo HAYASHI as an overflow too. */ 216*820c1a8dSHiroo HAYASHI #define INT_REMAINDER_RANGE_OVERFLOW(a, b, min, max) \ 217*820c1a8dSHiroo HAYASHI INT_DIVIDE_RANGE_OVERFLOW (a, b, min, max) 218*820c1a8dSHiroo HAYASHI 219*820c1a8dSHiroo HAYASHI /* Return 1 if A << B would overflow in [MIN,MAX] arithmetic. 220*820c1a8dSHiroo HAYASHI See above for restrictions. Here, MIN and MAX are for A only, and B need 221*820c1a8dSHiroo HAYASHI not be of the same type as the other arguments. The C standard says that 222*820c1a8dSHiroo HAYASHI behavior is undefined for shifts unless 0 <= B < wordwidth, and that when 223*820c1a8dSHiroo HAYASHI A is negative then A << B has undefined behavior and A >> B has 224*820c1a8dSHiroo HAYASHI implementation-defined behavior, but do not check these other 225*820c1a8dSHiroo HAYASHI restrictions. */ 226*820c1a8dSHiroo HAYASHI #define INT_LEFT_SHIFT_RANGE_OVERFLOW(a, b, min, max) \ 227*820c1a8dSHiroo HAYASHI ((a) < 0 \ 228*820c1a8dSHiroo HAYASHI ? (a) < (min) >> (b) \ 229*820c1a8dSHiroo HAYASHI : (max) >> (b) < (a)) 230*820c1a8dSHiroo HAYASHI 231*820c1a8dSHiroo HAYASHI /* True if __builtin_add_overflow (A, B, P) and __builtin_sub_overflow 232*820c1a8dSHiroo HAYASHI (A, B, P) work when P is non-null. */ 233*820c1a8dSHiroo HAYASHI /* __builtin_{add,sub}_overflow exists but is not reliable in GCC 5.x and 6.x, 234*820c1a8dSHiroo HAYASHI see <https://gcc.gnu.org/bugzilla/show_bug.cgi?id=98269>. */ 235*820c1a8dSHiroo HAYASHI #if 7 <= __GNUC__ && !defined __ICC 236*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_ADD_OVERFLOW 1 237*820c1a8dSHiroo HAYASHI #elif defined __has_builtin 238*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_ADD_OVERFLOW __has_builtin (__builtin_add_overflow) 239*820c1a8dSHiroo HAYASHI #else 240*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_ADD_OVERFLOW 0 241*820c1a8dSHiroo HAYASHI #endif 242*820c1a8dSHiroo HAYASHI 243*820c1a8dSHiroo HAYASHI /* True if __builtin_mul_overflow (A, B, P) works when P is non-null. */ 244*820c1a8dSHiroo HAYASHI #ifdef __clang__ 245*820c1a8dSHiroo HAYASHI /* Work around Clang bug <https://bugs.llvm.org/show_bug.cgi?id=16404>. */ 246*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_MUL_OVERFLOW 0 247*820c1a8dSHiroo HAYASHI #else 248*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_MUL_OVERFLOW _GL_HAS_BUILTIN_ADD_OVERFLOW 249*820c1a8dSHiroo HAYASHI #endif 250*820c1a8dSHiroo HAYASHI 251*820c1a8dSHiroo HAYASHI /* True if __builtin_add_overflow_p (A, B, C) works, and similarly for 252*820c1a8dSHiroo HAYASHI __builtin_sub_overflow_p and __builtin_mul_overflow_p. */ 253*820c1a8dSHiroo HAYASHI #if defined __clang__ || defined __ICC 254*820c1a8dSHiroo HAYASHI /* Clang 11 lacks __builtin_mul_overflow_p, and even if it did it 255*820c1a8dSHiroo HAYASHI would presumably run afoul of Clang bug 16404. ICC 2021.1's 256*820c1a8dSHiroo HAYASHI __builtin_add_overflow_p etc. are not treated as integral constant 257*820c1a8dSHiroo HAYASHI expressions even when all arguments are. */ 258*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_OVERFLOW_P 0 259*820c1a8dSHiroo HAYASHI #elif defined __has_builtin 260*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_OVERFLOW_P __has_builtin (__builtin_mul_overflow_p) 261*820c1a8dSHiroo HAYASHI #else 262*820c1a8dSHiroo HAYASHI # define _GL_HAS_BUILTIN_OVERFLOW_P (7 <= __GNUC__) 263*820c1a8dSHiroo HAYASHI #endif 264*820c1a8dSHiroo HAYASHI 265*820c1a8dSHiroo HAYASHI /* The _GL*_OVERFLOW macros have the same restrictions as the 266*820c1a8dSHiroo HAYASHI *_RANGE_OVERFLOW macros, except that they do not assume that operands 267*820c1a8dSHiroo HAYASHI (e.g., A and B) have the same type as MIN and MAX. Instead, they assume 268*820c1a8dSHiroo HAYASHI that the result (e.g., A + B) has that type. */ 269*820c1a8dSHiroo HAYASHI #if _GL_HAS_BUILTIN_OVERFLOW_P 270*820c1a8dSHiroo HAYASHI # define _GL_ADD_OVERFLOW(a, b, min, max) \ 271*820c1a8dSHiroo HAYASHI __builtin_add_overflow_p (a, b, (__typeof__ ((a) + (b))) 0) 272*820c1a8dSHiroo HAYASHI # define _GL_SUBTRACT_OVERFLOW(a, b, min, max) \ 273*820c1a8dSHiroo HAYASHI __builtin_sub_overflow_p (a, b, (__typeof__ ((a) - (b))) 0) 274*820c1a8dSHiroo HAYASHI # define _GL_MULTIPLY_OVERFLOW(a, b, min, max) \ 275*820c1a8dSHiroo HAYASHI __builtin_mul_overflow_p (a, b, (__typeof__ ((a) * (b))) 0) 276*820c1a8dSHiroo HAYASHI #else 277*820c1a8dSHiroo HAYASHI # define _GL_ADD_OVERFLOW(a, b, min, max) \ 278*820c1a8dSHiroo HAYASHI ((min) < 0 ? INT_ADD_RANGE_OVERFLOW (a, b, min, max) \ 279*820c1a8dSHiroo HAYASHI : (a) < 0 ? (b) <= (a) + (b) \ 280*820c1a8dSHiroo HAYASHI : (b) < 0 ? (a) <= (a) + (b) \ 281*820c1a8dSHiroo HAYASHI : (a) + (b) < (b)) 282*820c1a8dSHiroo HAYASHI # define _GL_SUBTRACT_OVERFLOW(a, b, min, max) \ 283*820c1a8dSHiroo HAYASHI ((min) < 0 ? INT_SUBTRACT_RANGE_OVERFLOW (a, b, min, max) \ 284*820c1a8dSHiroo HAYASHI : (a) < 0 ? 1 \ 285*820c1a8dSHiroo HAYASHI : (b) < 0 ? (a) - (b) <= (a) \ 286*820c1a8dSHiroo HAYASHI : (a) < (b)) 287*820c1a8dSHiroo HAYASHI # define _GL_MULTIPLY_OVERFLOW(a, b, min, max) \ 288*820c1a8dSHiroo HAYASHI (((min) == 0 && (((a) < 0 && 0 < (b)) || ((b) < 0 && 0 < (a)))) \ 289*820c1a8dSHiroo HAYASHI || INT_MULTIPLY_RANGE_OVERFLOW (a, b, min, max)) 290*820c1a8dSHiroo HAYASHI #endif 291*820c1a8dSHiroo HAYASHI #define _GL_DIVIDE_OVERFLOW(a, b, min, max) \ 292*820c1a8dSHiroo HAYASHI ((min) < 0 ? (b) == _GL_INT_NEGATE_CONVERT (min, 1) && (a) < - (max) \ 293*820c1a8dSHiroo HAYASHI : (a) < 0 ? (b) <= (a) + (b) - 1 \ 294*820c1a8dSHiroo HAYASHI : (b) < 0 && (a) + (b) <= (a)) 295*820c1a8dSHiroo HAYASHI #define _GL_REMAINDER_OVERFLOW(a, b, min, max) \ 296*820c1a8dSHiroo HAYASHI ((min) < 0 ? (b) == _GL_INT_NEGATE_CONVERT (min, 1) && (a) < - (max) \ 297*820c1a8dSHiroo HAYASHI : (a) < 0 ? (a) % (b) != ((max) - (b) + 1) % (b) \ 298*820c1a8dSHiroo HAYASHI : (b) < 0 && ! _GL_UNSIGNED_NEG_MULTIPLE (a, b, max)) 299*820c1a8dSHiroo HAYASHI 300*820c1a8dSHiroo HAYASHI /* Return a nonzero value if A is a mathematical multiple of B, where 301*820c1a8dSHiroo HAYASHI A is unsigned, B is negative, and MAX is the maximum value of A's 302*820c1a8dSHiroo HAYASHI type. A's type must be the same as (A % B)'s type. Normally (A % 303*820c1a8dSHiroo HAYASHI -B == 0) suffices, but things get tricky if -B would overflow. */ 304*820c1a8dSHiroo HAYASHI #define _GL_UNSIGNED_NEG_MULTIPLE(a, b, max) \ 305*820c1a8dSHiroo HAYASHI (((b) < -_GL_SIGNED_INT_MAXIMUM (b) \ 306*820c1a8dSHiroo HAYASHI ? (_GL_SIGNED_INT_MAXIMUM (b) == (max) \ 307*820c1a8dSHiroo HAYASHI ? (a) \ 308*820c1a8dSHiroo HAYASHI : (a) % (_GL_INT_CONVERT (a, _GL_SIGNED_INT_MAXIMUM (b)) + 1)) \ 309*820c1a8dSHiroo HAYASHI : (a) % - (b)) \ 310*820c1a8dSHiroo HAYASHI == 0) 311*820c1a8dSHiroo HAYASHI 312*820c1a8dSHiroo HAYASHI /* Check for integer overflow, and report low order bits of answer. 313*820c1a8dSHiroo HAYASHI 314*820c1a8dSHiroo HAYASHI The INT_<op>_OVERFLOW macros return 1 if the corresponding C operators 315*820c1a8dSHiroo HAYASHI might not yield numerically correct answers due to arithmetic overflow. 316*820c1a8dSHiroo HAYASHI The INT_<op>_WRAPV macros compute the low-order bits of the sum, 317*820c1a8dSHiroo HAYASHI difference, and product of two C integers, and return 1 if these 318*820c1a8dSHiroo HAYASHI low-order bits are not numerically correct. 319*820c1a8dSHiroo HAYASHI These macros work correctly on all known practical hosts, and do not rely 320*820c1a8dSHiroo HAYASHI on undefined behavior due to signed arithmetic overflow. 321*820c1a8dSHiroo HAYASHI 322*820c1a8dSHiroo HAYASHI Example usage, assuming A and B are long int: 323*820c1a8dSHiroo HAYASHI 324*820c1a8dSHiroo HAYASHI if (INT_MULTIPLY_OVERFLOW (a, b)) 325*820c1a8dSHiroo HAYASHI printf ("result would overflow\n"); 326*820c1a8dSHiroo HAYASHI else 327*820c1a8dSHiroo HAYASHI printf ("result is %ld (no overflow)\n", a * b); 328*820c1a8dSHiroo HAYASHI 329*820c1a8dSHiroo HAYASHI Example usage with WRAPV flavor: 330*820c1a8dSHiroo HAYASHI 331*820c1a8dSHiroo HAYASHI long int result; 332*820c1a8dSHiroo HAYASHI bool overflow = INT_MULTIPLY_WRAPV (a, b, &result); 333*820c1a8dSHiroo HAYASHI printf ("result is %ld (%s)\n", result, 334*820c1a8dSHiroo HAYASHI overflow ? "after overflow" : "no overflow"); 335*820c1a8dSHiroo HAYASHI 336*820c1a8dSHiroo HAYASHI Restrictions on these macros: 337*820c1a8dSHiroo HAYASHI 338*820c1a8dSHiroo HAYASHI These macros do not check for all possible numerical problems or 339*820c1a8dSHiroo HAYASHI undefined or unspecified behavior: they do not check for division 340*820c1a8dSHiroo HAYASHI by zero, for bad shift counts, or for shifting negative numbers. 341*820c1a8dSHiroo HAYASHI 342*820c1a8dSHiroo HAYASHI These macros may evaluate their arguments zero or multiple times, so the 343*820c1a8dSHiroo HAYASHI arguments should not have side effects. 344*820c1a8dSHiroo HAYASHI 345*820c1a8dSHiroo HAYASHI The WRAPV macros are not constant expressions. They support only 346*820c1a8dSHiroo HAYASHI +, binary -, and *. 347*820c1a8dSHiroo HAYASHI 348*820c1a8dSHiroo HAYASHI Because the WRAPV macros convert the result, they report overflow 349*820c1a8dSHiroo HAYASHI in different circumstances than the OVERFLOW macros do. For 350*820c1a8dSHiroo HAYASHI example, in the typical case with 16-bit 'short' and 32-bit 'int', 351*820c1a8dSHiroo HAYASHI if A, B and R are all of type 'short' then INT_ADD_OVERFLOW (A, B) 352*820c1a8dSHiroo HAYASHI returns false because the addition cannot overflow after A and B 353*820c1a8dSHiroo HAYASHI are converted to 'int', whereas INT_ADD_WRAPV (A, B, &R) returns 354*820c1a8dSHiroo HAYASHI true or false depending on whether the sum fits into 'short'. 355*820c1a8dSHiroo HAYASHI 356*820c1a8dSHiroo HAYASHI These macros are tuned for their last input argument being a constant. 357*820c1a8dSHiroo HAYASHI 358*820c1a8dSHiroo HAYASHI Return 1 if the integer expressions A * B, A - B, -A, A * B, A / B, 359*820c1a8dSHiroo HAYASHI A % B, and A << B would overflow, respectively. */ 360*820c1a8dSHiroo HAYASHI 361*820c1a8dSHiroo HAYASHI #define INT_ADD_OVERFLOW(a, b) \ 362*820c1a8dSHiroo HAYASHI _GL_BINARY_OP_OVERFLOW (a, b, _GL_ADD_OVERFLOW) 363*820c1a8dSHiroo HAYASHI #define INT_SUBTRACT_OVERFLOW(a, b) \ 364*820c1a8dSHiroo HAYASHI _GL_BINARY_OP_OVERFLOW (a, b, _GL_SUBTRACT_OVERFLOW) 365*820c1a8dSHiroo HAYASHI #if _GL_HAS_BUILTIN_OVERFLOW_P 366*820c1a8dSHiroo HAYASHI # define INT_NEGATE_OVERFLOW(a) INT_SUBTRACT_OVERFLOW (0, a) 367*820c1a8dSHiroo HAYASHI #else 368*820c1a8dSHiroo HAYASHI # define INT_NEGATE_OVERFLOW(a) \ 369*820c1a8dSHiroo HAYASHI INT_NEGATE_RANGE_OVERFLOW (a, _GL_INT_MINIMUM (a), _GL_INT_MAXIMUM (a)) 370*820c1a8dSHiroo HAYASHI #endif 371*820c1a8dSHiroo HAYASHI #define INT_MULTIPLY_OVERFLOW(a, b) \ 372*820c1a8dSHiroo HAYASHI _GL_BINARY_OP_OVERFLOW (a, b, _GL_MULTIPLY_OVERFLOW) 373*820c1a8dSHiroo HAYASHI #define INT_DIVIDE_OVERFLOW(a, b) \ 374*820c1a8dSHiroo HAYASHI _GL_BINARY_OP_OVERFLOW (a, b, _GL_DIVIDE_OVERFLOW) 375*820c1a8dSHiroo HAYASHI #define INT_REMAINDER_OVERFLOW(a, b) \ 376*820c1a8dSHiroo HAYASHI _GL_BINARY_OP_OVERFLOW (a, b, _GL_REMAINDER_OVERFLOW) 377*820c1a8dSHiroo HAYASHI #define INT_LEFT_SHIFT_OVERFLOW(a, b) \ 378*820c1a8dSHiroo HAYASHI INT_LEFT_SHIFT_RANGE_OVERFLOW (a, b, \ 379*820c1a8dSHiroo HAYASHI _GL_INT_MINIMUM (a), _GL_INT_MAXIMUM (a)) 380*820c1a8dSHiroo HAYASHI 381*820c1a8dSHiroo HAYASHI /* Return 1 if the expression A <op> B would overflow, 382*820c1a8dSHiroo HAYASHI where OP_RESULT_OVERFLOW (A, B, MIN, MAX) does the actual test, 383*820c1a8dSHiroo HAYASHI assuming MIN and MAX are the minimum and maximum for the result type. 384*820c1a8dSHiroo HAYASHI Arguments should be free of side effects. */ 385*820c1a8dSHiroo HAYASHI #define _GL_BINARY_OP_OVERFLOW(a, b, op_result_overflow) \ 386*820c1a8dSHiroo HAYASHI op_result_overflow (a, b, \ 387*820c1a8dSHiroo HAYASHI _GL_INT_MINIMUM (_GL_INT_CONVERT (a, b)), \ 388*820c1a8dSHiroo HAYASHI _GL_INT_MAXIMUM (_GL_INT_CONVERT (a, b))) 389*820c1a8dSHiroo HAYASHI 390*820c1a8dSHiroo HAYASHI /* Store the low-order bits of A + B, A - B, A * B, respectively, into *R. 391*820c1a8dSHiroo HAYASHI Return 1 if the result overflows. See above for restrictions. */ 392*820c1a8dSHiroo HAYASHI #if _GL_HAS_BUILTIN_ADD_OVERFLOW 393*820c1a8dSHiroo HAYASHI # define INT_ADD_WRAPV(a, b, r) __builtin_add_overflow (a, b, r) 394*820c1a8dSHiroo HAYASHI # define INT_SUBTRACT_WRAPV(a, b, r) __builtin_sub_overflow (a, b, r) 395*820c1a8dSHiroo HAYASHI #else 396*820c1a8dSHiroo HAYASHI # define INT_ADD_WRAPV(a, b, r) \ 397*820c1a8dSHiroo HAYASHI _GL_INT_OP_WRAPV (a, b, r, +, _GL_INT_ADD_RANGE_OVERFLOW) 398*820c1a8dSHiroo HAYASHI # define INT_SUBTRACT_WRAPV(a, b, r) \ 399*820c1a8dSHiroo HAYASHI _GL_INT_OP_WRAPV (a, b, r, -, _GL_INT_SUBTRACT_RANGE_OVERFLOW) 400*820c1a8dSHiroo HAYASHI #endif 401*820c1a8dSHiroo HAYASHI #if _GL_HAS_BUILTIN_MUL_OVERFLOW 402*820c1a8dSHiroo HAYASHI # if ((9 < __GNUC__ + (3 <= __GNUC_MINOR__) \ 403*820c1a8dSHiroo HAYASHI || (__GNUC__ == 8 && 4 <= __GNUC_MINOR__)) \ 404*820c1a8dSHiroo HAYASHI && !defined __ICC) 405*820c1a8dSHiroo HAYASHI # define INT_MULTIPLY_WRAPV(a, b, r) __builtin_mul_overflow (a, b, r) 406*820c1a8dSHiroo HAYASHI # else 407*820c1a8dSHiroo HAYASHI /* Work around GCC bug 91450. */ 408*820c1a8dSHiroo HAYASHI # define INT_MULTIPLY_WRAPV(a, b, r) \ 409*820c1a8dSHiroo HAYASHI ((!_GL_SIGNED_TYPE_OR_EXPR (*(r)) && EXPR_SIGNED (a) && EXPR_SIGNED (b) \ 410*820c1a8dSHiroo HAYASHI && _GL_INT_MULTIPLY_RANGE_OVERFLOW (a, b, 0, (__typeof__ (*(r))) -1)) \ 411*820c1a8dSHiroo HAYASHI ? ((void) __builtin_mul_overflow (a, b, r), 1) \ 412*820c1a8dSHiroo HAYASHI : __builtin_mul_overflow (a, b, r)) 413*820c1a8dSHiroo HAYASHI # endif 414*820c1a8dSHiroo HAYASHI #else 415*820c1a8dSHiroo HAYASHI # define INT_MULTIPLY_WRAPV(a, b, r) \ 416*820c1a8dSHiroo HAYASHI _GL_INT_OP_WRAPV (a, b, r, *, _GL_INT_MULTIPLY_RANGE_OVERFLOW) 417*820c1a8dSHiroo HAYASHI #endif 418*820c1a8dSHiroo HAYASHI 419*820c1a8dSHiroo HAYASHI /* Nonzero if this compiler has GCC bug 68193 or Clang bug 25390. See: 420*820c1a8dSHiroo HAYASHI https://gcc.gnu.org/bugzilla/show_bug.cgi?id=68193 421*820c1a8dSHiroo HAYASHI https://llvm.org/bugs/show_bug.cgi?id=25390 422*820c1a8dSHiroo HAYASHI For now, assume all versions of GCC-like compilers generate bogus 423*820c1a8dSHiroo HAYASHI warnings for _Generic. This matters only for compilers that 424*820c1a8dSHiroo HAYASHI lack relevant builtins. */ 425*820c1a8dSHiroo HAYASHI #if __GNUC__ || defined __clang__ 426*820c1a8dSHiroo HAYASHI # define _GL__GENERIC_BOGUS 1 427*820c1a8dSHiroo HAYASHI #else 428*820c1a8dSHiroo HAYASHI # define _GL__GENERIC_BOGUS 0 429*820c1a8dSHiroo HAYASHI #endif 430*820c1a8dSHiroo HAYASHI 431*820c1a8dSHiroo HAYASHI /* Store the low-order bits of A <op> B into *R, where OP specifies 432*820c1a8dSHiroo HAYASHI the operation and OVERFLOW the overflow predicate. Return 1 if the 433*820c1a8dSHiroo HAYASHI result overflows. See above for restrictions. */ 434*820c1a8dSHiroo HAYASHI #if 201112 <= __STDC_VERSION__ && !_GL__GENERIC_BOGUS 435*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV(a, b, r, op, overflow) \ 436*820c1a8dSHiroo HAYASHI (_Generic \ 437*820c1a8dSHiroo HAYASHI (*(r), \ 438*820c1a8dSHiroo HAYASHI signed char: \ 439*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 440*820c1a8dSHiroo HAYASHI signed char, SCHAR_MIN, SCHAR_MAX), \ 441*820c1a8dSHiroo HAYASHI unsigned char: \ 442*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 443*820c1a8dSHiroo HAYASHI unsigned char, 0, UCHAR_MAX), \ 444*820c1a8dSHiroo HAYASHI short int: \ 445*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 446*820c1a8dSHiroo HAYASHI short int, SHRT_MIN, SHRT_MAX), \ 447*820c1a8dSHiroo HAYASHI unsigned short int: \ 448*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 449*820c1a8dSHiroo HAYASHI unsigned short int, 0, USHRT_MAX), \ 450*820c1a8dSHiroo HAYASHI int: \ 451*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 452*820c1a8dSHiroo HAYASHI int, INT_MIN, INT_MAX), \ 453*820c1a8dSHiroo HAYASHI unsigned int: \ 454*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 455*820c1a8dSHiroo HAYASHI unsigned int, 0, UINT_MAX), \ 456*820c1a8dSHiroo HAYASHI long int: \ 457*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 458*820c1a8dSHiroo HAYASHI long int, LONG_MIN, LONG_MAX), \ 459*820c1a8dSHiroo HAYASHI unsigned long int: \ 460*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 461*820c1a8dSHiroo HAYASHI unsigned long int, 0, ULONG_MAX), \ 462*820c1a8dSHiroo HAYASHI long long int: \ 463*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long long int, \ 464*820c1a8dSHiroo HAYASHI long long int, LLONG_MIN, LLONG_MAX), \ 465*820c1a8dSHiroo HAYASHI unsigned long long int: \ 466*820c1a8dSHiroo HAYASHI _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long long int, \ 467*820c1a8dSHiroo HAYASHI unsigned long long int, 0, ULLONG_MAX))) 468*820c1a8dSHiroo HAYASHI #else 469*820c1a8dSHiroo HAYASHI /* Store the low-order bits of A <op> B into *R, where OP specifies 470*820c1a8dSHiroo HAYASHI the operation and OVERFLOW the overflow predicate. If *R is 471*820c1a8dSHiroo HAYASHI signed, its type is ST with bounds SMIN..SMAX; otherwise its type 472*820c1a8dSHiroo HAYASHI is UT with bounds U..UMAX. ST and UT are narrower than int. 473*820c1a8dSHiroo HAYASHI Return 1 if the result overflows. See above for restrictions. */ 474*820c1a8dSHiroo HAYASHI # if _GL_HAVE___TYPEOF__ 475*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV_SMALLISH(a,b,r,op,overflow,st,smin,smax,ut,umax) \ 476*820c1a8dSHiroo HAYASHI (TYPE_SIGNED (__typeof__ (*(r))) \ 477*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, st, smin, smax) \ 478*820c1a8dSHiroo HAYASHI : _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, ut, 0, umax)) 479*820c1a8dSHiroo HAYASHI # else 480*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV_SMALLISH(a,b,r,op,overflow,st,smin,smax,ut,umax) \ 481*820c1a8dSHiroo HAYASHI (overflow (a, b, smin, smax) \ 482*820c1a8dSHiroo HAYASHI ? (overflow (a, b, 0, umax) \ 483*820c1a8dSHiroo HAYASHI ? (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a,b,op,unsigned,st), 1) \ 484*820c1a8dSHiroo HAYASHI : (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a,b,op,unsigned,st)) < 0) \ 485*820c1a8dSHiroo HAYASHI : (overflow (a, b, 0, umax) \ 486*820c1a8dSHiroo HAYASHI ? (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a,b,op,unsigned,st)) >= 0 \ 487*820c1a8dSHiroo HAYASHI : (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a,b,op,unsigned,st), 0))) 488*820c1a8dSHiroo HAYASHI # endif 489*820c1a8dSHiroo HAYASHI 490*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV(a, b, r, op, overflow) \ 491*820c1a8dSHiroo HAYASHI (sizeof *(r) == sizeof (signed char) \ 492*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_WRAPV_SMALLISH (a, b, r, op, overflow, \ 493*820c1a8dSHiroo HAYASHI signed char, SCHAR_MIN, SCHAR_MAX, \ 494*820c1a8dSHiroo HAYASHI unsigned char, UCHAR_MAX) \ 495*820c1a8dSHiroo HAYASHI : sizeof *(r) == sizeof (short int) \ 496*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_WRAPV_SMALLISH (a, b, r, op, overflow, \ 497*820c1a8dSHiroo HAYASHI short int, SHRT_MIN, SHRT_MAX, \ 498*820c1a8dSHiroo HAYASHI unsigned short int, USHRT_MAX) \ 499*820c1a8dSHiroo HAYASHI : sizeof *(r) == sizeof (int) \ 500*820c1a8dSHiroo HAYASHI ? (EXPR_SIGNED (*(r)) \ 501*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 502*820c1a8dSHiroo HAYASHI int, INT_MIN, INT_MAX) \ 503*820c1a8dSHiroo HAYASHI : _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned int, \ 504*820c1a8dSHiroo HAYASHI unsigned int, 0, UINT_MAX)) \ 505*820c1a8dSHiroo HAYASHI : _GL_INT_OP_WRAPV_LONGISH(a, b, r, op, overflow)) 506*820c1a8dSHiroo HAYASHI # ifdef LLONG_MAX 507*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV_LONGISH(a, b, r, op, overflow) \ 508*820c1a8dSHiroo HAYASHI (sizeof *(r) == sizeof (long int) \ 509*820c1a8dSHiroo HAYASHI ? (EXPR_SIGNED (*(r)) \ 510*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 511*820c1a8dSHiroo HAYASHI long int, LONG_MIN, LONG_MAX) \ 512*820c1a8dSHiroo HAYASHI : _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 513*820c1a8dSHiroo HAYASHI unsigned long int, 0, ULONG_MAX)) \ 514*820c1a8dSHiroo HAYASHI : (EXPR_SIGNED (*(r)) \ 515*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long long int, \ 516*820c1a8dSHiroo HAYASHI long long int, LLONG_MIN, LLONG_MAX) \ 517*820c1a8dSHiroo HAYASHI : _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long long int, \ 518*820c1a8dSHiroo HAYASHI unsigned long long int, 0, ULLONG_MAX))) 519*820c1a8dSHiroo HAYASHI # else 520*820c1a8dSHiroo HAYASHI # define _GL_INT_OP_WRAPV_LONGISH(a, b, r, op, overflow) \ 521*820c1a8dSHiroo HAYASHI (EXPR_SIGNED (*(r)) \ 522*820c1a8dSHiroo HAYASHI ? _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 523*820c1a8dSHiroo HAYASHI long int, LONG_MIN, LONG_MAX) \ 524*820c1a8dSHiroo HAYASHI : _GL_INT_OP_CALC (a, b, r, op, overflow, unsigned long int, \ 525*820c1a8dSHiroo HAYASHI unsigned long int, 0, ULONG_MAX)) 526*820c1a8dSHiroo HAYASHI # endif 527*820c1a8dSHiroo HAYASHI #endif 528*820c1a8dSHiroo HAYASHI 529*820c1a8dSHiroo HAYASHI /* Store the low-order bits of A <op> B into *R, where the operation 530*820c1a8dSHiroo HAYASHI is given by OP. Use the unsigned type UT for calculation to avoid 531*820c1a8dSHiroo HAYASHI overflow problems. *R's type is T, with extrema TMIN and TMAX. 532*820c1a8dSHiroo HAYASHI T must be a signed integer type. Return 1 if the result overflows. */ 533*820c1a8dSHiroo HAYASHI #define _GL_INT_OP_CALC(a, b, r, op, overflow, ut, t, tmin, tmax) \ 534*820c1a8dSHiroo HAYASHI (overflow (a, b, tmin, tmax) \ 535*820c1a8dSHiroo HAYASHI ? (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a, b, op, ut, t), 1) \ 536*820c1a8dSHiroo HAYASHI : (*(r) = _GL_INT_OP_WRAPV_VIA_UNSIGNED (a, b, op, ut, t), 0)) 537*820c1a8dSHiroo HAYASHI 538*820c1a8dSHiroo HAYASHI /* Return the low-order bits of A <op> B, where the operation is given 539*820c1a8dSHiroo HAYASHI by OP. Use the unsigned type UT for calculation to avoid undefined 540*820c1a8dSHiroo HAYASHI behavior on signed integer overflow, and convert the result to type T. 541*820c1a8dSHiroo HAYASHI UT is at least as wide as T and is no narrower than unsigned int, 542*820c1a8dSHiroo HAYASHI T is two's complement, and there is no padding or trap representations. 543*820c1a8dSHiroo HAYASHI Assume that converting UT to T yields the low-order bits, as is 544*820c1a8dSHiroo HAYASHI done in all known two's-complement C compilers. E.g., see: 545*820c1a8dSHiroo HAYASHI https://gcc.gnu.org/onlinedocs/gcc/Integers-implementation.html 546*820c1a8dSHiroo HAYASHI 547*820c1a8dSHiroo HAYASHI According to the C standard, converting UT to T yields an 548*820c1a8dSHiroo HAYASHI implementation-defined result or signal for values outside T's 549*820c1a8dSHiroo HAYASHI range. However, code that works around this theoretical problem 550*820c1a8dSHiroo HAYASHI runs afoul of a compiler bug in Oracle Studio 12.3 x86. See: 551*820c1a8dSHiroo HAYASHI https://lists.gnu.org/r/bug-gnulib/2017-04/msg00049.html 552*820c1a8dSHiroo HAYASHI As the compiler bug is real, don't try to work around the 553*820c1a8dSHiroo HAYASHI theoretical problem. */ 554*820c1a8dSHiroo HAYASHI 555*820c1a8dSHiroo HAYASHI #define _GL_INT_OP_WRAPV_VIA_UNSIGNED(a, b, op, ut, t) \ 556*820c1a8dSHiroo HAYASHI ((t) ((ut) (a) op (ut) (b))) 557*820c1a8dSHiroo HAYASHI 558*820c1a8dSHiroo HAYASHI /* Return true if the numeric values A + B, A - B, A * B fall outside 559*820c1a8dSHiroo HAYASHI the range TMIN..TMAX. Arguments should be integer expressions 560*820c1a8dSHiroo HAYASHI without side effects. TMIN should be signed and nonpositive. 561*820c1a8dSHiroo HAYASHI TMAX should be positive, and should be signed unless TMIN is zero. */ 562*820c1a8dSHiroo HAYASHI #define _GL_INT_ADD_RANGE_OVERFLOW(a, b, tmin, tmax) \ 563*820c1a8dSHiroo HAYASHI ((b) < 0 \ 564*820c1a8dSHiroo HAYASHI ? (((tmin) \ 565*820c1a8dSHiroo HAYASHI ? ((EXPR_SIGNED (_GL_INT_CONVERT (a, (tmin) - (b))) || (b) < (tmin)) \ 566*820c1a8dSHiroo HAYASHI && (a) < (tmin) - (b)) \ 567*820c1a8dSHiroo HAYASHI : (a) <= -1 - (b)) \ 568*820c1a8dSHiroo HAYASHI || ((EXPR_SIGNED (a) ? 0 <= (a) : (tmax) < (a)) && (tmax) < (a) + (b))) \ 569*820c1a8dSHiroo HAYASHI : (a) < 0 \ 570*820c1a8dSHiroo HAYASHI ? (((tmin) \ 571*820c1a8dSHiroo HAYASHI ? ((EXPR_SIGNED (_GL_INT_CONVERT (b, (tmin) - (a))) || (a) < (tmin)) \ 572*820c1a8dSHiroo HAYASHI && (b) < (tmin) - (a)) \ 573*820c1a8dSHiroo HAYASHI : (b) <= -1 - (a)) \ 574*820c1a8dSHiroo HAYASHI || ((EXPR_SIGNED (_GL_INT_CONVERT (a, b)) || (tmax) < (b)) \ 575*820c1a8dSHiroo HAYASHI && (tmax) < (a) + (b))) \ 576*820c1a8dSHiroo HAYASHI : (tmax) < (b) || (tmax) - (b) < (a)) 577*820c1a8dSHiroo HAYASHI #define _GL_INT_SUBTRACT_RANGE_OVERFLOW(a, b, tmin, tmax) \ 578*820c1a8dSHiroo HAYASHI (((a) < 0) == ((b) < 0) \ 579*820c1a8dSHiroo HAYASHI ? ((a) < (b) \ 580*820c1a8dSHiroo HAYASHI ? !(tmin) || -1 - (tmin) < (b) - (a) - 1 \ 581*820c1a8dSHiroo HAYASHI : (tmax) < (a) - (b)) \ 582*820c1a8dSHiroo HAYASHI : (a) < 0 \ 583*820c1a8dSHiroo HAYASHI ? ((!EXPR_SIGNED (_GL_INT_CONVERT ((a) - (tmin), b)) && (a) - (tmin) < 0) \ 584*820c1a8dSHiroo HAYASHI || (a) - (tmin) < (b)) \ 585*820c1a8dSHiroo HAYASHI : ((! (EXPR_SIGNED (_GL_INT_CONVERT (tmax, b)) \ 586*820c1a8dSHiroo HAYASHI && EXPR_SIGNED (_GL_INT_CONVERT ((tmax) + (b), a))) \ 587*820c1a8dSHiroo HAYASHI && (tmax) <= -1 - (b)) \ 588*820c1a8dSHiroo HAYASHI || (tmax) + (b) < (a))) 589*820c1a8dSHiroo HAYASHI #define _GL_INT_MULTIPLY_RANGE_OVERFLOW(a, b, tmin, tmax) \ 590*820c1a8dSHiroo HAYASHI ((b) < 0 \ 591*820c1a8dSHiroo HAYASHI ? ((a) < 0 \ 592*820c1a8dSHiroo HAYASHI ? (EXPR_SIGNED (_GL_INT_CONVERT (tmax, b)) \ 593*820c1a8dSHiroo HAYASHI ? (a) < (tmax) / (b) \ 594*820c1a8dSHiroo HAYASHI : ((INT_NEGATE_OVERFLOW (b) \ 595*820c1a8dSHiroo HAYASHI ? _GL_INT_CONVERT (b, tmax) >> (TYPE_WIDTH (+ (b)) - 1) \ 596*820c1a8dSHiroo HAYASHI : (tmax) / -(b)) \ 597*820c1a8dSHiroo HAYASHI <= -1 - (a))) \ 598*820c1a8dSHiroo HAYASHI : INT_NEGATE_OVERFLOW (_GL_INT_CONVERT (b, tmin)) && (b) == -1 \ 599*820c1a8dSHiroo HAYASHI ? (EXPR_SIGNED (a) \ 600*820c1a8dSHiroo HAYASHI ? 0 < (a) + (tmin) \ 601*820c1a8dSHiroo HAYASHI : 0 < (a) && -1 - (tmin) < (a) - 1) \ 602*820c1a8dSHiroo HAYASHI : (tmin) / (b) < (a)) \ 603*820c1a8dSHiroo HAYASHI : (b) == 0 \ 604*820c1a8dSHiroo HAYASHI ? 0 \ 605*820c1a8dSHiroo HAYASHI : ((a) < 0 \ 606*820c1a8dSHiroo HAYASHI ? (INT_NEGATE_OVERFLOW (_GL_INT_CONVERT (a, tmin)) && (a) == -1 \ 607*820c1a8dSHiroo HAYASHI ? (EXPR_SIGNED (b) ? 0 < (b) + (tmin) : -1 - (tmin) < (b) - 1) \ 608*820c1a8dSHiroo HAYASHI : (tmin) / (a) < (b)) \ 609*820c1a8dSHiroo HAYASHI : (tmax) / (b) < (a))) 610*820c1a8dSHiroo HAYASHI 611*820c1a8dSHiroo HAYASHI /* The following macros compute A + B, A - B, and A * B, respectively. 612*820c1a8dSHiroo HAYASHI If no overflow occurs, they set *R to the result and return 1; 613*820c1a8dSHiroo HAYASHI otherwise, they return 0 and may modify *R. 614*820c1a8dSHiroo HAYASHI 615*820c1a8dSHiroo HAYASHI Example usage: 616*820c1a8dSHiroo HAYASHI 617*820c1a8dSHiroo HAYASHI long int result; 618*820c1a8dSHiroo HAYASHI if (INT_ADD_OK (a, b, &result)) 619*820c1a8dSHiroo HAYASHI printf ("result is %ld\n", result); 620*820c1a8dSHiroo HAYASHI else 621*820c1a8dSHiroo HAYASHI printf ("overflow\n"); 622*820c1a8dSHiroo HAYASHI 623*820c1a8dSHiroo HAYASHI A, B, and *R should be integers; they need not be the same type, 624*820c1a8dSHiroo HAYASHI and they need not be all signed or all unsigned. 625*820c1a8dSHiroo HAYASHI 626*820c1a8dSHiroo HAYASHI These macros work correctly on all known practical hosts, and do not rely 627*820c1a8dSHiroo HAYASHI on undefined behavior due to signed arithmetic overflow. 628*820c1a8dSHiroo HAYASHI 629*820c1a8dSHiroo HAYASHI These macros are not constant expressions. 630*820c1a8dSHiroo HAYASHI 631*820c1a8dSHiroo HAYASHI These macros may evaluate their arguments zero or multiple times, so the 632*820c1a8dSHiroo HAYASHI arguments should not have side effects. 633*820c1a8dSHiroo HAYASHI 634*820c1a8dSHiroo HAYASHI These macros are tuned for B being a constant. */ 635*820c1a8dSHiroo HAYASHI 636*820c1a8dSHiroo HAYASHI #define INT_ADD_OK(a, b, r) ! INT_ADD_WRAPV (a, b, r) 637*820c1a8dSHiroo HAYASHI #define INT_SUBTRACT_OK(a, b, r) ! INT_SUBTRACT_WRAPV (a, b, r) 638*820c1a8dSHiroo HAYASHI #define INT_MULTIPLY_OK(a, b, r) ! INT_MULTIPLY_WRAPV (a, b, r) 639*820c1a8dSHiroo HAYASHI 640*820c1a8dSHiroo HAYASHI #endif /* _GL_INTPROPS_H */ 641