1 /* Software floating-point emulation. 2 Copyright (C) 1997-2016 Free Software Foundation, Inc. 3 This file is part of the GNU C Library. 4 Contributed by Richard Henderson (rth@cygnus.com), 5 Jakub Jelinek (jj@ultra.linux.cz), 6 David S. Miller (davem@redhat.com) and 7 Peter Maydell (pmaydell@chiark.greenend.org.uk). 8 9 The GNU C Library is free software; you can redistribute it and/or 10 modify it under the terms of the GNU Lesser General Public 11 License as published by the Free Software Foundation; either 12 version 2.1 of the License, or (at your option) any later version. 13 14 In addition to the permissions in the GNU Lesser General Public 15 License, the Free Software Foundation gives you unlimited 16 permission to link the compiled version of this file into 17 combinations with other programs, and to distribute those 18 combinations without any restriction coming from the use of this 19 file. (The Lesser General Public License restrictions do apply in 20 other respects; for example, they cover modification of the file, 21 and distribution when not linked into a combine executable.) 22 23 The GNU C Library is distributed in the hope that it will be useful, 24 but WITHOUT ANY WARRANTY; without even the implied warranty of 25 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 26 Lesser General Public License for more details. 27 28 You should have received a copy of the GNU Lesser General Public 29 License along with the GNU C Library; if not, see 30 <http://www.gnu.org/licenses/>. */ 31 32 #ifndef SOFT_FP_H 33 #define SOFT_FP_H 1 34 35 #ifdef _LIBC 36 # include <sfp-machine.h> 37 #elif defined __KERNEL__ 38 /* The Linux kernel uses asm/ names for architecture-specific 39 files. */ 40 # include <asm/sfp-machine.h> 41 #else 42 # include "sfp-machine.h" 43 #endif 44 45 /* Allow sfp-machine to have its own byte order definitions. */ 46 #ifndef __BYTE_ORDER 47 # ifdef _LIBC 48 # include <endian.h> 49 # else 50 # error "endianness not defined by sfp-machine.h" 51 # endif 52 #endif 53 54 /* For unreachable default cases in switch statements over bitwise OR 55 of FP_CLS_* values. */ 56 #if (defined __GNUC__ \ 57 && (__GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 5))) 58 # define _FP_UNREACHABLE __builtin_unreachable () 59 #else 60 # define _FP_UNREACHABLE abort () 61 #endif 62 63 #if ((defined __GNUC__ \ 64 && (__GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6))) \ 65 || (defined __STDC_VERSION__ && __STDC_VERSION__ >= 201112L)) 66 # define _FP_STATIC_ASSERT(expr, msg) \ 67 _Static_assert ((expr), msg) 68 #else 69 # define _FP_STATIC_ASSERT(expr, msg) \ 70 extern int (*__Static_assert_function (void)) \ 71 [!!sizeof (struct { int __error_if_negative: (expr) ? 2 : -1; })] 72 #endif 73 74 /* In the Linux kernel, some architectures have a single function that 75 uses different kinds of unpacking and packing depending on the 76 instruction being emulated, meaning it is not readily visible to 77 the compiler that variables from _FP_DECL and _FP_FRAC_DECL_* 78 macros are only used in cases where they were initialized. */ 79 #ifdef __KERNEL__ 80 # define _FP_ZERO_INIT = 0 81 #else 82 # define _FP_ZERO_INIT 83 #endif 84 85 #define _FP_WORKBITS 3 86 #define _FP_WORK_LSB ((_FP_W_TYPE) 1 << 3) 87 #define _FP_WORK_ROUND ((_FP_W_TYPE) 1 << 2) 88 #define _FP_WORK_GUARD ((_FP_W_TYPE) 1 << 1) 89 #define _FP_WORK_STICKY ((_FP_W_TYPE) 1 << 0) 90 91 #ifndef FP_RND_NEAREST 92 # define FP_RND_NEAREST 0 93 # define FP_RND_ZERO 1 94 # define FP_RND_PINF 2 95 # define FP_RND_MINF 3 96 #endif 97 #ifndef FP_ROUNDMODE 98 # define FP_ROUNDMODE FP_RND_NEAREST 99 #endif 100 101 /* By default don't care about exceptions. */ 102 #ifndef FP_EX_INVALID 103 # define FP_EX_INVALID 0 104 #endif 105 #ifndef FP_EX_OVERFLOW 106 # define FP_EX_OVERFLOW 0 107 #endif 108 #ifndef FP_EX_UNDERFLOW 109 # define FP_EX_UNDERFLOW 0 110 #endif 111 #ifndef FP_EX_DIVZERO 112 # define FP_EX_DIVZERO 0 113 #endif 114 #ifndef FP_EX_INEXACT 115 # define FP_EX_INEXACT 0 116 #endif 117 #ifndef FP_EX_DENORM 118 # define FP_EX_DENORM 0 119 #endif 120 121 /* Sub-exceptions of "invalid". */ 122 /* Signaling NaN operand. */ 123 #ifndef FP_EX_INVALID_SNAN 124 # define FP_EX_INVALID_SNAN 0 125 #endif 126 /* Inf * 0. */ 127 #ifndef FP_EX_INVALID_IMZ 128 # define FP_EX_INVALID_IMZ 0 129 #endif 130 /* fma (Inf, 0, c). */ 131 #ifndef FP_EX_INVALID_IMZ_FMA 132 # define FP_EX_INVALID_IMZ_FMA 0 133 #endif 134 /* Inf - Inf. */ 135 #ifndef FP_EX_INVALID_ISI 136 # define FP_EX_INVALID_ISI 0 137 #endif 138 /* 0 / 0. */ 139 #ifndef FP_EX_INVALID_ZDZ 140 # define FP_EX_INVALID_ZDZ 0 141 #endif 142 /* Inf / Inf. */ 143 #ifndef FP_EX_INVALID_IDI 144 # define FP_EX_INVALID_IDI 0 145 #endif 146 /* sqrt (negative). */ 147 #ifndef FP_EX_INVALID_SQRT 148 # define FP_EX_INVALID_SQRT 0 149 #endif 150 /* Invalid conversion to integer. */ 151 #ifndef FP_EX_INVALID_CVI 152 # define FP_EX_INVALID_CVI 0 153 #endif 154 /* Invalid comparison. */ 155 #ifndef FP_EX_INVALID_VC 156 # define FP_EX_INVALID_VC 0 157 #endif 158 159 /* _FP_STRUCT_LAYOUT may be defined as an attribute to determine the 160 struct layout variant used for structures where bit-fields are used 161 to access specific parts of binary floating-point numbers. This is 162 required for systems where the default ABI uses struct layout with 163 differences in how consecutive bit-fields are laid out from the 164 default expected by soft-fp. */ 165 #ifndef _FP_STRUCT_LAYOUT 166 # define _FP_STRUCT_LAYOUT 167 #endif 168 169 #ifdef _FP_DECL_EX 170 # define FP_DECL_EX \ 171 int _fex = 0; \ 172 _FP_DECL_EX 173 #else 174 # define FP_DECL_EX int _fex = 0 175 #endif 176 177 /* Initialize any machine-specific state used in FP_ROUNDMODE, 178 FP_TRAPPING_EXCEPTIONS or FP_HANDLE_EXCEPTIONS. */ 179 #ifndef FP_INIT_ROUNDMODE 180 # define FP_INIT_ROUNDMODE do {} while (0) 181 #endif 182 183 /* Initialize any machine-specific state used in 184 FP_TRAPPING_EXCEPTIONS or FP_HANDLE_EXCEPTIONS. */ 185 #ifndef FP_INIT_TRAPPING_EXCEPTIONS 186 # define FP_INIT_TRAPPING_EXCEPTIONS FP_INIT_ROUNDMODE 187 #endif 188 189 /* Initialize any machine-specific state used in 190 FP_HANDLE_EXCEPTIONS. */ 191 #ifndef FP_INIT_EXCEPTIONS 192 # define FP_INIT_EXCEPTIONS FP_INIT_TRAPPING_EXCEPTIONS 193 #endif 194 195 #ifndef FP_HANDLE_EXCEPTIONS 196 # define FP_HANDLE_EXCEPTIONS do {} while (0) 197 #endif 198 199 /* Whether to flush subnormal inputs to zero with the same sign. */ 200 #ifndef FP_DENORM_ZERO 201 # define FP_DENORM_ZERO 0 202 #endif 203 204 #ifndef FP_INHIBIT_RESULTS 205 /* By default we write the results always. 206 sfp-machine may override this and e.g. 207 check if some exceptions are unmasked 208 and inhibit it in such a case. */ 209 # define FP_INHIBIT_RESULTS 0 210 #endif 211 212 #define FP_SET_EXCEPTION(ex) \ 213 _fex |= (ex) 214 215 #define FP_CUR_EXCEPTIONS \ 216 (_fex) 217 218 #ifndef FP_TRAPPING_EXCEPTIONS 219 # define FP_TRAPPING_EXCEPTIONS 0 220 #endif 221 222 /* A file using soft-fp may define FP_NO_EXCEPTIONS before including 223 soft-fp.h to indicate that, although a macro used there could raise 224 exceptions, or do rounding and potentially thereby raise 225 exceptions, for some arguments, for the particular arguments used 226 in that file no exceptions or rounding can occur. Such a file 227 should not itself use macros relating to handling exceptions and 228 rounding modes; this is only for indirect uses (in particular, in 229 _FP_FROM_INT and the macros it calls). */ 230 #ifdef FP_NO_EXCEPTIONS 231 232 # undef FP_SET_EXCEPTION 233 # define FP_SET_EXCEPTION(ex) do {} while (0) 234 235 # undef FP_CUR_EXCEPTIONS 236 # define FP_CUR_EXCEPTIONS 0 237 238 # undef FP_TRAPPING_EXCEPTIONS 239 # define FP_TRAPPING_EXCEPTIONS 0 240 241 # undef FP_ROUNDMODE 242 # define FP_ROUNDMODE FP_RND_ZERO 243 244 # undef _FP_TININESS_AFTER_ROUNDING 245 # define _FP_TININESS_AFTER_ROUNDING 0 246 247 #endif 248 249 /* A file using soft-fp may define FP_NO_EXACT_UNDERFLOW before 250 including soft-fp.h to indicate that, although a macro used there 251 could allow for the case of exact underflow requiring the underflow 252 exception to be raised if traps are enabled, for the particular 253 arguments used in that file no exact underflow can occur. */ 254 #ifdef FP_NO_EXACT_UNDERFLOW 255 # undef FP_TRAPPING_EXCEPTIONS 256 # define FP_TRAPPING_EXCEPTIONS 0 257 #endif 258 259 #define _FP_ROUND_NEAREST(wc, X) \ 260 do \ 261 { \ 262 if ((_FP_FRAC_LOW_##wc (X) & 15) != _FP_WORK_ROUND) \ 263 _FP_FRAC_ADDI_##wc (X, _FP_WORK_ROUND); \ 264 } \ 265 while (0) 266 267 #define _FP_ROUND_ZERO(wc, X) (void) 0 268 269 #define _FP_ROUND_PINF(wc, X) \ 270 do \ 271 { \ 272 if (!X##_s && (_FP_FRAC_LOW_##wc (X) & 7)) \ 273 _FP_FRAC_ADDI_##wc (X, _FP_WORK_LSB); \ 274 } \ 275 while (0) 276 277 #define _FP_ROUND_MINF(wc, X) \ 278 do \ 279 { \ 280 if (X##_s && (_FP_FRAC_LOW_##wc (X) & 7)) \ 281 _FP_FRAC_ADDI_##wc (X, _FP_WORK_LSB); \ 282 } \ 283 while (0) 284 285 #define _FP_ROUND(wc, X) \ 286 do \ 287 { \ 288 if (_FP_FRAC_LOW_##wc (X) & 7) \ 289 { \ 290 FP_SET_EXCEPTION (FP_EX_INEXACT); \ 291 switch (FP_ROUNDMODE) \ 292 { \ 293 case FP_RND_NEAREST: \ 294 _FP_ROUND_NEAREST (wc, X); \ 295 break; \ 296 case FP_RND_ZERO: \ 297 _FP_ROUND_ZERO (wc, X); \ 298 break; \ 299 case FP_RND_PINF: \ 300 _FP_ROUND_PINF (wc, X); \ 301 break; \ 302 case FP_RND_MINF: \ 303 _FP_ROUND_MINF (wc, X); \ 304 break; \ 305 } \ 306 } \ 307 } \ 308 while (0) 309 310 #define FP_CLS_NORMAL 0 311 #define FP_CLS_ZERO 1 312 #define FP_CLS_INF 2 313 #define FP_CLS_NAN 3 314 315 #define _FP_CLS_COMBINE(x, y) (((x) << 2) | (y)) 316 317 #include "op-1.h" 318 #include "op-2.h" 319 #include "op-4.h" 320 #include "op-8.h" 321 #include "op-common.h" 322 323 /* Sigh. Silly things longlong.h needs. */ 324 #define UWtype _FP_W_TYPE 325 #define W_TYPE_SIZE _FP_W_TYPE_SIZE 326 327 typedef int QItype __attribute__ ((mode (QI))); 328 typedef int SItype __attribute__ ((mode (SI))); 329 typedef int DItype __attribute__ ((mode (DI))); 330 typedef unsigned int UQItype __attribute__ ((mode (QI))); 331 typedef unsigned int USItype __attribute__ ((mode (SI))); 332 typedef unsigned int UDItype __attribute__ ((mode (DI))); 333 #if _FP_W_TYPE_SIZE == 32 334 typedef unsigned int UHWtype __attribute__ ((mode (HI))); 335 #elif _FP_W_TYPE_SIZE == 64 336 typedef USItype UHWtype; 337 #endif 338 339 #ifndef CMPtype 340 # define CMPtype int 341 #endif 342 343 #define SI_BITS (__CHAR_BIT__ * (int) sizeof (SItype)) 344 #define DI_BITS (__CHAR_BIT__ * (int) sizeof (DItype)) 345 346 #ifndef umul_ppmm 347 # ifdef _LIBC 348 # include <stdlib/longlong.h> 349 # else 350 # include "longlong.h" 351 # endif 352 #endif 353 354 #endif /* !SOFT_FP_H */ 355