raw
smg_comms_c_wrappers    1 /* longlong.h -- definitions for mixed size 32/64 bit arithmetic.
smg_comms_c_wrappers 2 * Modified by No Such Labs. (C) 2015. See README.
smg_comms_c_wrappers 3 *
smg_comms_c_wrappers 4 * This file was originally part of Gnu Privacy Guard (GPG), ver. 1.4.10,
smg_comms_c_wrappers 5 * SHA256(gnupg-1.4.10.tar.gz):
smg_comms_c_wrappers 6 * 0bfd74660a2f6cedcf7d8256db4a63c996ffebbcdc2cf54397bfb72878c5a85a
smg_comms_c_wrappers 7 * (C) 1994-2005 Free Software Foundation, Inc.
smg_comms_c_wrappers 8 *
smg_comms_c_wrappers 9 * This program is free software: you can redistribute it and/or modify
smg_comms_c_wrappers 10 * it under the terms of the GNU General Public License as published by
smg_comms_c_wrappers 11 * the Free Software Foundation, either version 3 of the License, or
smg_comms_c_wrappers 12 * (at your option) any later version.
smg_comms_c_wrappers 13 *
smg_comms_c_wrappers 14 * This program is distributed in the hope that it will be useful,
smg_comms_c_wrappers 15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
smg_comms_c_wrappers 16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
smg_comms_c_wrappers 17 * GNU General Public License for more details.
smg_comms_c_wrappers 18 *
smg_comms_c_wrappers 19 * You should have received a copy of the GNU General Public License
smg_comms_c_wrappers 20 * along with this program. If not, see <http://www.gnu.org/licenses/>.
smg_comms_c_wrappers 21 */
smg_comms_c_wrappers 22
smg_comms_c_wrappers 23 /* You have to define the following before including this file:
smg_comms_c_wrappers 24
smg_comms_c_wrappers 25 UWtype -- An unsigned type, default type for operations (typically a "word")
smg_comms_c_wrappers 26 UHWtype -- An unsigned type, at least half the size of UWtype.
smg_comms_c_wrappers 27 UDWtype -- An unsigned type, at least twice as large a UWtype
smg_comms_c_wrappers 28 W_TYPE_SIZE -- size in bits of UWtype
smg_comms_c_wrappers 29
smg_comms_c_wrappers 30 SItype, USItype -- Signed and unsigned 32 bit types.
smg_comms_c_wrappers 31 DItype, UDItype -- Signed and unsigned 64 bit types.
smg_comms_c_wrappers 32
smg_comms_c_wrappers 33 On a 32 bit machine UWtype should typically be USItype;
smg_comms_c_wrappers 34 on a 64 bit machine, UWtype should typically be UDItype.
smg_comms_c_wrappers 35 */
smg_comms_c_wrappers 36
smg_comms_c_wrappers 37 #define __BITS4 (W_TYPE_SIZE / 4)
smg_comms_c_wrappers 38 #define __ll_B ((UWtype) 1 << (W_TYPE_SIZE / 2))
smg_comms_c_wrappers 39 #define __ll_lowpart(t) ((UWtype) (t) & (__ll_B - 1))
smg_comms_c_wrappers 40 #define __ll_highpart(t) ((UWtype) (t) >> (W_TYPE_SIZE / 2))
smg_comms_c_wrappers 41
smg_comms_c_wrappers 42 /* This is used to make sure no undesirable sharing between different libraries
smg_comms_c_wrappers 43 that use this file takes place. */
smg_comms_c_wrappers 44 #ifndef __MPN
smg_comms_c_wrappers 45 #define __MPN(x) __##x
smg_comms_c_wrappers 46 #endif
smg_comms_c_wrappers 47
smg_comms_c_wrappers 48 /***************************************
smg_comms_c_wrappers 49 *********** Generic Versions ********
smg_comms_c_wrappers 50 ***************************************/
smg_comms_c_wrappers 51 #if !defined (umul_ppmm) && defined (__umulsidi3)
smg_comms_c_wrappers 52 #define umul_ppmm(ph, pl, m0, m1) \
smg_comms_c_wrappers 53 { \
smg_comms_c_wrappers 54 UDWtype __ll = __umulsidi3 (m0, m1); \
smg_comms_c_wrappers 55 ph = (UWtype) (__ll >> W_TYPE_SIZE); \
smg_comms_c_wrappers 56 pl = (UWtype) __ll; \
smg_comms_c_wrappers 57 }
smg_comms_c_wrappers 58 #endif
smg_comms_c_wrappers 59
smg_comms_c_wrappers 60 #if !defined (__umulsidi3)
smg_comms_c_wrappers 61 #define __umulsidi3(u, v) \
smg_comms_c_wrappers 62 ({UWtype __hi, __lo; \
smg_comms_c_wrappers 63 umul_ppmm (__hi, __lo, u, v); \
smg_comms_c_wrappers 64 ((UDWtype) __hi << W_TYPE_SIZE) | __lo; })
smg_comms_c_wrappers 65 #endif
smg_comms_c_wrappers 66
smg_comms_c_wrappers 67 /* If this machine has no inline assembler, use C macros. */
smg_comms_c_wrappers 68
smg_comms_c_wrappers 69 #if !defined (add_ssaaaa)
smg_comms_c_wrappers 70 #define add_ssaaaa(sh, sl, ah, al, bh, bl) \
smg_comms_c_wrappers 71 do { \
smg_comms_c_wrappers 72 UWtype __x; \
smg_comms_c_wrappers 73 __x = (al) + (bl); \
smg_comms_c_wrappers 74 (sh) = (ah) + (bh) + (__x < (al)); \
smg_comms_c_wrappers 75 (sl) = __x; \
smg_comms_c_wrappers 76 } while (0)
smg_comms_c_wrappers 77 #endif
smg_comms_c_wrappers 78
smg_comms_c_wrappers 79 #if !defined (sub_ddmmss)
smg_comms_c_wrappers 80 #define sub_ddmmss(sh, sl, ah, al, bh, bl) \
smg_comms_c_wrappers 81 do { \
smg_comms_c_wrappers 82 UWtype __x; \
smg_comms_c_wrappers 83 __x = (al) - (bl); \
smg_comms_c_wrappers 84 (sh) = (ah) - (bh) - (__x > (al)); \
smg_comms_c_wrappers 85 (sl) = __x; \
smg_comms_c_wrappers 86 } while (0)
smg_comms_c_wrappers 87 #endif
smg_comms_c_wrappers 88
smg_comms_c_wrappers 89 #if !defined (umul_ppmm)
smg_comms_c_wrappers 90 #define umul_ppmm(w1, w0, u, v) \
smg_comms_c_wrappers 91 do { \
smg_comms_c_wrappers 92 UWtype __x0, __x1, __x2, __x3; \
smg_comms_c_wrappers 93 UHWtype __ul, __vl, __uh, __vh; \
smg_comms_c_wrappers 94 UWtype __u = (u), __v = (v); \
smg_comms_c_wrappers 95 \
smg_comms_c_wrappers 96 __ul = __ll_lowpart (__u); \
smg_comms_c_wrappers 97 __uh = __ll_highpart (__u); \
smg_comms_c_wrappers 98 __vl = __ll_lowpart (__v); \
smg_comms_c_wrappers 99 __vh = __ll_highpart (__v); \
smg_comms_c_wrappers 100 \
smg_comms_c_wrappers 101 __x0 = (UWtype) __ul * __vl; \
smg_comms_c_wrappers 102 __x1 = (UWtype) __ul * __vh; \
smg_comms_c_wrappers 103 __x2 = (UWtype) __uh * __vl; \
smg_comms_c_wrappers 104 __x3 = (UWtype) __uh * __vh; \
smg_comms_c_wrappers 105 \
smg_comms_c_wrappers 106 __x1 += __ll_highpart (__x0);/* this can't give carry */ \
smg_comms_c_wrappers 107 __x1 += __x2; /* but this indeed can */ \
smg_comms_c_wrappers 108 if (__x1 < __x2) /* did we get it? */ \
smg_comms_c_wrappers 109 __x3 += __ll_B; /* yes, add it in the proper pos. */ \
smg_comms_c_wrappers 110 \
smg_comms_c_wrappers 111 (w1) = __x3 + __ll_highpart (__x1); \
smg_comms_c_wrappers 112 (w0) = (__ll_lowpart (__x1) << W_TYPE_SIZE/2) + __ll_lowpart (__x0);\
smg_comms_c_wrappers 113 } while (0)
smg_comms_c_wrappers 114 #endif
smg_comms_c_wrappers 115
smg_comms_c_wrappers 116 #if !defined (umul_ppmm)
smg_comms_c_wrappers 117 #define smul_ppmm(w1, w0, u, v) \
smg_comms_c_wrappers 118 do { \
smg_comms_c_wrappers 119 UWtype __w1; \
smg_comms_c_wrappers 120 UWtype __m0 = (u), __m1 = (v); \
smg_comms_c_wrappers 121 umul_ppmm (__w1, w0, __m0, __m1); \
smg_comms_c_wrappers 122 (w1) = __w1 - (-(__m0 >> (W_TYPE_SIZE - 1)) & __m1) \
smg_comms_c_wrappers 123 - (-(__m1 >> (W_TYPE_SIZE - 1)) & __m0); \
smg_comms_c_wrappers 124 } while (0)
smg_comms_c_wrappers 125 #endif
smg_comms_c_wrappers 126
smg_comms_c_wrappers 127 /* Define this unconditionally, so it can be used for debugging. */
smg_comms_c_wrappers 128 #define __udiv_qrnnd_c(q, r, n1, n0, d) \
smg_comms_c_wrappers 129 do { \
smg_comms_c_wrappers 130 UWtype __d1, __d0, __q1, __q0, __r1, __r0, __m; \
smg_comms_c_wrappers 131 __d1 = __ll_highpart (d); \
smg_comms_c_wrappers 132 __d0 = __ll_lowpart (d); \
smg_comms_c_wrappers 133 \
smg_comms_c_wrappers 134 __r1 = (n1) % __d1; \
smg_comms_c_wrappers 135 __q1 = (n1) / __d1; \
smg_comms_c_wrappers 136 __m = (UWtype) __q1 * __d0; \
smg_comms_c_wrappers 137 __r1 = __r1 * __ll_B | __ll_highpart (n0); \
smg_comms_c_wrappers 138 if (__r1 < __m) \
smg_comms_c_wrappers 139 { \
smg_comms_c_wrappers 140 __q1--, __r1 += (d); \
smg_comms_c_wrappers 141 if (__r1 >= (d)) /* i.e. we didn't get carry when adding to __r1 */\
smg_comms_c_wrappers 142 if (__r1 < __m) \
smg_comms_c_wrappers 143 __q1--, __r1 += (d); \
smg_comms_c_wrappers 144 } \
smg_comms_c_wrappers 145 __r1 -= __m; \
smg_comms_c_wrappers 146 \
smg_comms_c_wrappers 147 __r0 = __r1 % __d1; \
smg_comms_c_wrappers 148 __q0 = __r1 / __d1; \
smg_comms_c_wrappers 149 __m = (UWtype) __q0 * __d0; \
smg_comms_c_wrappers 150 __r0 = __r0 * __ll_B | __ll_lowpart (n0); \
smg_comms_c_wrappers 151 if (__r0 < __m) \
smg_comms_c_wrappers 152 { \
smg_comms_c_wrappers 153 __q0--, __r0 += (d); \
smg_comms_c_wrappers 154 if (__r0 >= (d)) \
smg_comms_c_wrappers 155 if (__r0 < __m) \
smg_comms_c_wrappers 156 __q0--, __r0 += (d); \
smg_comms_c_wrappers 157 } \
smg_comms_c_wrappers 158 __r0 -= __m; \
smg_comms_c_wrappers 159 \
smg_comms_c_wrappers 160 (q) = (UWtype) __q1 * __ll_B | __q0; \
smg_comms_c_wrappers 161 (r) = __r0; \
smg_comms_c_wrappers 162 } while (0)
smg_comms_c_wrappers 163
smg_comms_c_wrappers 164 /* If the processor has no udiv_qrnnd but sdiv_qrnnd, go through
smg_comms_c_wrappers 165 __udiv_w_sdiv (defined in libgcc or elsewhere). */
smg_comms_c_wrappers 166 #if !defined (udiv_qrnnd) && defined (sdiv_qrnnd)
smg_comms_c_wrappers 167 #define udiv_qrnnd(q, r, nh, nl, d) \
smg_comms_c_wrappers 168 do { \
smg_comms_c_wrappers 169 UWtype __r; \
smg_comms_c_wrappers 170 (q) = __MPN(udiv_w_sdiv) (&__r, nh, nl, d); \
smg_comms_c_wrappers 171 (r) = __r; \
smg_comms_c_wrappers 172 } while (0)
smg_comms_c_wrappers 173 #endif
smg_comms_c_wrappers 174
smg_comms_c_wrappers 175 /* If udiv_qrnnd was not defined for this processor, use __udiv_qrnnd_c. */
smg_comms_c_wrappers 176 #if !defined (udiv_qrnnd)
smg_comms_c_wrappers 177 #define UDIV_NEEDS_NORMALIZATION 1
smg_comms_c_wrappers 178 #define udiv_qrnnd __udiv_qrnnd_c
smg_comms_c_wrappers 179 #endif
smg_comms_c_wrappers 180
smg_comms_c_wrappers 181 #if !defined (count_leading_zeros)
smg_comms_c_wrappers 182 extern
smg_comms_c_wrappers 183 #ifdef __STDC__
smg_comms_c_wrappers 184 const
smg_comms_c_wrappers 185 #endif
smg_comms_c_wrappers 186 unsigned char __clz_tab[];
smg_comms_c_wrappers 187 #define MPI_INTERNAL_NEED_CLZ_TAB 1
smg_comms_c_wrappers 188 #define count_leading_zeros(count, x) \
smg_comms_c_wrappers 189 do { \
smg_comms_c_wrappers 190 UWtype __xr = (x); \
smg_comms_c_wrappers 191 UWtype __a; \
smg_comms_c_wrappers 192 \
smg_comms_c_wrappers 193 if (W_TYPE_SIZE <= 32) \
smg_comms_c_wrappers 194 { \
smg_comms_c_wrappers 195 __a = __xr < ((UWtype) 1 << 2*__BITS4) \
smg_comms_c_wrappers 196 ? (__xr < ((UWtype) 1 << __BITS4) ? 0 : __BITS4) \
smg_comms_c_wrappers 197 : (__xr < ((UWtype) 1 << 3*__BITS4) ? 2*__BITS4 : 3*__BITS4);\
smg_comms_c_wrappers 198 } \
smg_comms_c_wrappers 199 else \
smg_comms_c_wrappers 200 { \
smg_comms_c_wrappers 201 for (__a = W_TYPE_SIZE - 8; __a > 0; __a -= 8) \
smg_comms_c_wrappers 202 if (((__xr >> __a) & 0xff) != 0) \
smg_comms_c_wrappers 203 break; \
smg_comms_c_wrappers 204 } \
smg_comms_c_wrappers 205 \
smg_comms_c_wrappers 206 (count) = W_TYPE_SIZE - (__clz_tab[__xr >> __a] + __a); \
smg_comms_c_wrappers 207 } while (0)
smg_comms_c_wrappers 208 /* This version gives a well-defined value for zero. */
smg_comms_c_wrappers 209 #define COUNT_LEADING_ZEROS_0 W_TYPE_SIZE
smg_comms_c_wrappers 210 #endif
smg_comms_c_wrappers 211
smg_comms_c_wrappers 212 #if !defined (count_trailing_zeros)
smg_comms_c_wrappers 213 /* Define count_trailing_zeros using count_leading_zeros. The latter might be
smg_comms_c_wrappers 214 defined in asm, but if it is not, the C version above is good enough. */
smg_comms_c_wrappers 215 #define count_trailing_zeros(count, x) \
smg_comms_c_wrappers 216 do { \
smg_comms_c_wrappers 217 UWtype __ctz_x = (x); \
smg_comms_c_wrappers 218 UWtype __ctz_c; \
smg_comms_c_wrappers 219 count_leading_zeros (__ctz_c, __ctz_x & -__ctz_x); \
smg_comms_c_wrappers 220 (count) = W_TYPE_SIZE - 1 - __ctz_c; \
smg_comms_c_wrappers 221 } while (0)
smg_comms_c_wrappers 222 #endif
smg_comms_c_wrappers 223
smg_comms_c_wrappers 224 #ifndef UDIV_NEEDS_NORMALIZATION
smg_comms_c_wrappers 225 #define UDIV_NEEDS_NORMALIZATION 0
smg_comms_c_wrappers 226 #endif