1/* Software floating-point emulation. 2 Definitions for IEEE Quad Precision. 3 Copyright (C) 1997-2014 Free Software Foundation, Inc. 4 This file is part of the GNU C Library. 5 Contributed by Richard Henderson (rth@cygnus.com), 6 Jakub Jelinek (jj@ultra.linux.cz), 7 David S. Miller (davem@redhat.com) and 8 Peter Maydell (pmaydell@chiark.greenend.org.uk). 9 10 The GNU C Library is free software; you can redistribute it and/or 11 modify it under the terms of the GNU Lesser General Public 12 License as published by the Free Software Foundation; either 13 version 2.1 of the License, or (at your option) any later version. 14 15 In addition to the permissions in the GNU Lesser General Public 16 License, the Free Software Foundation gives you unlimited 17 permission to link the compiled version of this file into 18 combinations with other programs, and to distribute those 19 combinations without any restriction coming from the use of this 20 file. (The Lesser General Public License restrictions do apply in 21 other respects; for example, they cover modification of the file, 22 and distribution when not linked into a combine executable.) 23 24 The GNU C Library is distributed in the hope that it will be useful, 25 but WITHOUT ANY WARRANTY; without even the implied warranty of 26 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 27 Lesser General Public License for more details. 28 29 You should have received a copy of the GNU Lesser General Public 30 License along with the GNU C Library; if not, see 31 <http://www.gnu.org/licenses/>. */ 32 33#if _FP_W_TYPE_SIZE < 32 34# error "Here's a nickel, kid. Go buy yourself a real computer." 35#endif 36 37#if _FP_W_TYPE_SIZE < 64 38# define _FP_FRACTBITS_Q (4*_FP_W_TYPE_SIZE) 39# define _FP_FRACTBITS_DW_Q (8*_FP_W_TYPE_SIZE) 40#else 41# define _FP_FRACTBITS_Q (2*_FP_W_TYPE_SIZE) 42# define _FP_FRACTBITS_DW_Q (4*_FP_W_TYPE_SIZE) 43#endif 44 45#define _FP_FRACBITS_Q 113 46#define _FP_FRACXBITS_Q (_FP_FRACTBITS_Q - _FP_FRACBITS_Q) 47#define _FP_WFRACBITS_Q (_FP_WORKBITS + _FP_FRACBITS_Q) 48#define _FP_WFRACXBITS_Q (_FP_FRACTBITS_Q - _FP_WFRACBITS_Q) 49#define _FP_EXPBITS_Q 15 50#define _FP_EXPBIAS_Q 16383 51#define _FP_EXPMAX_Q 32767 52 53#define _FP_QNANBIT_Q \ 54 ((_FP_W_TYPE) 1 << (_FP_FRACBITS_Q-2) % _FP_W_TYPE_SIZE) 55#define _FP_QNANBIT_SH_Q \ 56 ((_FP_W_TYPE) 1 << (_FP_FRACBITS_Q-2+_FP_WORKBITS) % _FP_W_TYPE_SIZE) 57#define _FP_IMPLBIT_Q \ 58 ((_FP_W_TYPE) 1 << (_FP_FRACBITS_Q-1) % _FP_W_TYPE_SIZE) 59#define _FP_IMPLBIT_SH_Q \ 60 ((_FP_W_TYPE) 1 << (_FP_FRACBITS_Q-1+_FP_WORKBITS) % _FP_W_TYPE_SIZE) 61#define _FP_OVERFLOW_Q \ 62 ((_FP_W_TYPE) 1 << (_FP_WFRACBITS_Q % _FP_W_TYPE_SIZE)) 63 64#define _FP_WFRACBITS_DW_Q (2 * _FP_WFRACBITS_Q) 65#define _FP_WFRACXBITS_DW_Q (_FP_FRACTBITS_DW_Q - _FP_WFRACBITS_DW_Q) 66#define _FP_HIGHBIT_DW_Q \ 67 ((_FP_W_TYPE) 1 << (_FP_WFRACBITS_DW_Q - 1) % _FP_W_TYPE_SIZE) 68 69typedef float TFtype __attribute__ ((mode (TF))); 70 71#if _FP_W_TYPE_SIZE < 64 72 73union _FP_UNION_Q 74{ 75 TFtype flt; 76 struct _FP_STRUCT_LAYOUT 77 { 78# if __BYTE_ORDER == __BIG_ENDIAN 79 unsigned sign : 1; 80 unsigned exp : _FP_EXPBITS_Q; 81 unsigned long frac3 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0)-(_FP_W_TYPE_SIZE * 3); 82 unsigned long frac2 : _FP_W_TYPE_SIZE; 83 unsigned long frac1 : _FP_W_TYPE_SIZE; 84 unsigned long frac0 : _FP_W_TYPE_SIZE; 85# else 86 unsigned long frac0 : _FP_W_TYPE_SIZE; 87 unsigned long frac1 : _FP_W_TYPE_SIZE; 88 unsigned long frac2 : _FP_W_TYPE_SIZE; 89 unsigned long frac3 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0)-(_FP_W_TYPE_SIZE * 3); 90 unsigned exp : _FP_EXPBITS_Q; 91 unsigned sign : 1; 92# endif /* not bigendian */ 93 } bits __attribute__ ((packed)); 94}; 95 96 97# define FP_DECL_Q(X) _FP_DECL (4, X) 98# define FP_UNPACK_RAW_Q(X, val) _FP_UNPACK_RAW_4 (Q, X, (val)) 99# define FP_UNPACK_RAW_QP(X, val) _FP_UNPACK_RAW_4_P (Q, X, (val)) 100# define FP_PACK_RAW_Q(val, X) _FP_PACK_RAW_4 (Q, (val), X) 101# define FP_PACK_RAW_QP(val, X) \ 102 do \ 103 { \ 104 if (!FP_INHIBIT_RESULTS) \ 105 _FP_PACK_RAW_4_P (Q, (val), X); \ 106 } \ 107 while (0) 108 109# define FP_UNPACK_Q(X, val) \ 110 do \ 111 { \ 112 _FP_UNPACK_RAW_4 (Q, X, (val)); \ 113 _FP_UNPACK_CANONICAL (Q, 4, X); \ 114 } \ 115 while (0) 116 117# define FP_UNPACK_QP(X, val) \ 118 do \ 119 { \ 120 _FP_UNPACK_RAW_4_P (Q, X, (val)); \ 121 _FP_UNPACK_CANONICAL (Q, 4, X); \ 122 } \ 123 while (0) 124 125# define FP_UNPACK_SEMIRAW_Q(X, val) \ 126 do \ 127 { \ 128 _FP_UNPACK_RAW_4 (Q, X, (val)); \ 129 _FP_UNPACK_SEMIRAW (Q, 4, X); \ 130 } \ 131 while (0) 132 133# define FP_UNPACK_SEMIRAW_QP(X, val) \ 134 do \ 135 { \ 136 _FP_UNPACK_RAW_4_P (Q, X, (val)); \ 137 _FP_UNPACK_SEMIRAW (Q, 4, X); \ 138 } \ 139 while (0) 140 141# define FP_PACK_Q(val, X) \ 142 do \ 143 { \ 144 _FP_PACK_CANONICAL (Q, 4, X); \ 145 _FP_PACK_RAW_4 (Q, (val), X); \ 146 } \ 147 while (0) 148 149# define FP_PACK_QP(val, X) \ 150 do \ 151 { \ 152 _FP_PACK_CANONICAL (Q, 4, X); \ 153 if (!FP_INHIBIT_RESULTS) \ 154 _FP_PACK_RAW_4_P (Q, (val), X); \ 155 } \ 156 while (0) 157 158# define FP_PACK_SEMIRAW_Q(val, X) \ 159 do \ 160 { \ 161 _FP_PACK_SEMIRAW (Q, 4, X); \ 162 _FP_PACK_RAW_4 (Q, (val), X); \ 163 } \ 164 while (0) 165 166# define FP_PACK_SEMIRAW_QP(val, X) \ 167 do \ 168 { \ 169 _FP_PACK_SEMIRAW (Q, 4, X); \ 170 if (!FP_INHIBIT_RESULTS) \ 171 _FP_PACK_RAW_4_P (Q, (val), X); \ 172 } \ 173 while (0) 174 175# define FP_ISSIGNAN_Q(X) _FP_ISSIGNAN (Q, 4, X) 176# define FP_NEG_Q(R, X) _FP_NEG (Q, 4, R, X) 177# define FP_ADD_Q(R, X, Y) _FP_ADD (Q, 4, R, X, Y) 178# define FP_SUB_Q(R, X, Y) _FP_SUB (Q, 4, R, X, Y) 179# define FP_MUL_Q(R, X, Y) _FP_MUL (Q, 4, R, X, Y) 180# define FP_DIV_Q(R, X, Y) _FP_DIV (Q, 4, R, X, Y) 181# define FP_SQRT_Q(R, X) _FP_SQRT (Q, 4, R, X) 182# define _FP_SQRT_MEAT_Q(R, S, T, X, Q) _FP_SQRT_MEAT_4 (R, S, T, X, (Q)) 183# define FP_FMA_Q(R, X, Y, Z) _FP_FMA (Q, 4, 8, R, X, Y, Z) 184 185# define FP_CMP_Q(r, X, Y, un, ex) _FP_CMP (Q, 4, (r), X, Y, (un), (ex)) 186# define FP_CMP_EQ_Q(r, X, Y, ex) _FP_CMP_EQ (Q, 4, (r), X, Y, (ex)) 187# define FP_CMP_UNORD_Q(r, X, Y, ex) _FP_CMP_UNORD (Q, 4, (r), X, Y, (ex)) 188 189# define FP_TO_INT_Q(r, X, rsz, rsg) _FP_TO_INT (Q, 4, (r), X, (rsz), (rsg)) 190# define FP_FROM_INT_Q(X, r, rs, rt) _FP_FROM_INT (Q, 4, X, (r), (rs), rt) 191 192# define _FP_FRAC_HIGH_Q(X) _FP_FRAC_HIGH_4 (X) 193# define _FP_FRAC_HIGH_RAW_Q(X) _FP_FRAC_HIGH_4 (X) 194 195# define _FP_FRAC_HIGH_DW_Q(X) _FP_FRAC_HIGH_8 (X) 196 197#else /* not _FP_W_TYPE_SIZE < 64 */ 198union _FP_UNION_Q 199{ 200 TFtype flt /* __attribute__ ((mode (TF))) */ ; 201 struct _FP_STRUCT_LAYOUT 202 { 203 _FP_W_TYPE a, b; 204 } longs; 205 struct _FP_STRUCT_LAYOUT 206 { 207# if __BYTE_ORDER == __BIG_ENDIAN 208 unsigned sign : 1; 209 unsigned exp : _FP_EXPBITS_Q; 210 _FP_W_TYPE frac1 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0) - _FP_W_TYPE_SIZE; 211 _FP_W_TYPE frac0 : _FP_W_TYPE_SIZE; 212# else 213 _FP_W_TYPE frac0 : _FP_W_TYPE_SIZE; 214 _FP_W_TYPE frac1 : _FP_FRACBITS_Q - (_FP_IMPLBIT_Q != 0) - _FP_W_TYPE_SIZE; 215 unsigned exp : _FP_EXPBITS_Q; 216 unsigned sign : 1; 217# endif 218 } bits; 219}; 220 221# define FP_DECL_Q(X) _FP_DECL (2, X) 222# define FP_UNPACK_RAW_Q(X, val) _FP_UNPACK_RAW_2 (Q, X, (val)) 223# define FP_UNPACK_RAW_QP(X, val) _FP_UNPACK_RAW_2_P (Q, X, (val)) 224# define FP_PACK_RAW_Q(val, X) _FP_PACK_RAW_2 (Q, (val), X) 225# define FP_PACK_RAW_QP(val, X) \ 226 do \ 227 { \ 228 if (!FP_INHIBIT_RESULTS) \ 229 _FP_PACK_RAW_2_P (Q, (val), X); \ 230 } \ 231 while (0) 232 233# define FP_UNPACK_Q(X, val) \ 234 do \ 235 { \ 236 _FP_UNPACK_RAW_2 (Q, X, (val)); \ 237 _FP_UNPACK_CANONICAL (Q, 2, X); \ 238 } \ 239 while (0) 240 241# define FP_UNPACK_QP(X, val) \ 242 do \ 243 { \ 244 _FP_UNPACK_RAW_2_P (Q, X, (val)); \ 245 _FP_UNPACK_CANONICAL (Q, 2, X); \ 246 } \ 247 while (0) 248 249# define FP_UNPACK_SEMIRAW_Q(X, val) \ 250 do \ 251 { \ 252 _FP_UNPACK_RAW_2 (Q, X, (val)); \ 253 _FP_UNPACK_SEMIRAW (Q, 2, X); \ 254 } \ 255 while (0) 256 257# define FP_UNPACK_SEMIRAW_QP(X, val) \ 258 do \ 259 { \ 260 _FP_UNPACK_RAW_2_P (Q, X, (val)); \ 261 _FP_UNPACK_SEMIRAW (Q, 2, X); \ 262 } \ 263 while (0) 264 265# define FP_PACK_Q(val, X) \ 266 do \ 267 { \ 268 _FP_PACK_CANONICAL (Q, 2, X); \ 269 _FP_PACK_RAW_2 (Q, (val), X); \ 270 } \ 271 while (0) 272 273# define FP_PACK_QP(val, X) \ 274 do \ 275 { \ 276 _FP_PACK_CANONICAL (Q, 2, X); \ 277 if (!FP_INHIBIT_RESULTS) \ 278 _FP_PACK_RAW_2_P (Q, (val), X); \ 279 } \ 280 while (0) 281 282# define FP_PACK_SEMIRAW_Q(val, X) \ 283 do \ 284 { \ 285 _FP_PACK_SEMIRAW (Q, 2, X); \ 286 _FP_PACK_RAW_2 (Q, (val), X); \ 287 } \ 288 while (0) 289 290# define FP_PACK_SEMIRAW_QP(val, X) \ 291 do \ 292 { \ 293 _FP_PACK_SEMIRAW (Q, 2, X); \ 294 if (!FP_INHIBIT_RESULTS) \ 295 _FP_PACK_RAW_2_P (Q, (val), X); \ 296 } \ 297 while (0) 298 299# define FP_ISSIGNAN_Q(X) _FP_ISSIGNAN (Q, 2, X) 300# define FP_NEG_Q(R, X) _FP_NEG (Q, 2, R, X) 301# define FP_ADD_Q(R, X, Y) _FP_ADD (Q, 2, R, X, Y) 302# define FP_SUB_Q(R, X, Y) _FP_SUB (Q, 2, R, X, Y) 303# define FP_MUL_Q(R, X, Y) _FP_MUL (Q, 2, R, X, Y) 304# define FP_DIV_Q(R, X, Y) _FP_DIV (Q, 2, R, X, Y) 305# define FP_SQRT_Q(R, X) _FP_SQRT (Q, 2, R, X) 306# define _FP_SQRT_MEAT_Q(R, S, T, X, Q) _FP_SQRT_MEAT_2 (R, S, T, X, (Q)) 307# define FP_FMA_Q(R, X, Y, Z) _FP_FMA (Q, 2, 4, R, X, Y, Z) 308 309# define FP_CMP_Q(r, X, Y, un, ex) _FP_CMP (Q, 2, (r), X, Y, (un), (ex)) 310# define FP_CMP_EQ_Q(r, X, Y, ex) _FP_CMP_EQ (Q, 2, (r), X, Y, (ex)) 311# define FP_CMP_UNORD_Q(r, X, Y, ex) _FP_CMP_UNORD (Q, 2, (r), X, Y, (ex)) 312 313# define FP_TO_INT_Q(r, X, rsz, rsg) _FP_TO_INT (Q, 2, (r), X, (rsz), (rsg)) 314# define FP_FROM_INT_Q(X, r, rs, rt) _FP_FROM_INT (Q, 2, X, (r), (rs), rt) 315 316# define _FP_FRAC_HIGH_Q(X) _FP_FRAC_HIGH_2 (X) 317# define _FP_FRAC_HIGH_RAW_Q(X) _FP_FRAC_HIGH_2 (X) 318 319# define _FP_FRAC_HIGH_DW_Q(X) _FP_FRAC_HIGH_4 (X) 320 321#endif /* not _FP_W_TYPE_SIZE < 64 */ 322