libgcc/fp-bit.h

   1 /* Header file for fp-bit.c.  */
   2 /* Copyright (C) 2000-2017 Free Software Foundation, Inc.
   3
   4 This file is part of GCC.
   5
   6 GCC is free software; you can redistribute it and/or modify it under
   7 the terms of the GNU General Public License as published by the Free
   8 Software Foundation; either version 3, or (at your option) any later
   9 version.
  10
  11 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
  12 WARRANTY; without even the implied warranty of MERCHANTABILITY or
  13 FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  14 for more details.
  15
  16 Under Section 7 of GPL version 3, you are granted additional
  17 permissions described in the GCC Runtime Library Exception, version
  18 3.1, as published by the Free Software Foundation.
  19
  20 You should have received a copy of the GNU General Public License and
  21 a copy of the GCC Runtime Library Exception along with this program;
  22 see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
  23 <http://www.gnu.org/licenses/>.  */
  24
  25 #ifndef GCC_FP_BIT_H
  26 #define GCC_FP_BIT_H
  27
  28 /* Defining FINE_GRAINED_LIBRARIES allows one to select which routines
  29    from this file are compiled via additional -D options.
  30
  31    This avoids the need to pull in the entire fp emulation library
  32    when only a small number of functions are needed.
  33
  34    If FINE_GRAINED_LIBRARIES is not defined, then compile every
  35    suitable routine.  */
  36 #ifndef FINE_GRAINED_LIBRARIES
  37 #define L_pack_df
  38 #define L_unpack_df
  39 #define L_pack_sf
  40 #define L_unpack_sf
  41 #define L_addsub_sf
  42 #define L_addsub_df
  43 #define L_mul_sf
  44 #define L_mul_df
  45 #define L_div_sf
  46 #define L_div_df
  47 #define L_fpcmp_parts_sf
  48 #define L_fpcmp_parts_df
  49 #define L_compare_sf
  50 #define L_compare_df
  51 #define L_eq_sf
  52 #define L_eq_df
  53 #define L_ne_sf
  54 #define L_ne_df
  55 #define L_gt_sf
  56 #define L_gt_df
  57 #define L_ge_sf
  58 #define L_ge_df
  59 #define L_lt_sf
  60 #define L_lt_df
  61 #define L_le_sf
  62 #define L_le_df
  63 #define L_unord_sf
  64 #define L_unord_df
  65 #define L_usi_to_sf
  66 #define L_usi_to_df
  67 #define L_si_to_sf
  68 #define L_si_to_df
  69 #define L_sf_to_si
  70 #define L_df_to_si
  71 #define L_f_to_usi
  72 #define L_df_to_usi
  73 #define L_negate_sf
  74 #define L_negate_df
  75 #define L_make_sf
  76 #define L_make_df
  77 #define L_sf_to_df
  78 #define L_df_to_sf
  79 #ifdef FLOAT
  80 #define L_thenan_sf
  81 #else
  82 #define L_thenan_df
  83 #endif
  84 #endif /* ! FINE_GRAINED_LIBRARIES */
  85
  86 #if __LDBL_MANT_DIG__ == 113 || __LDBL_MANT_DIG__ == 106
  87 # if defined(TFLOAT) || defined(L_sf_to_tf) || defined(L_df_to_tf)
  88 #  define TMODES
  89 # endif
  90 #endif
  91
  92 typedef float SFtype __attribute__ ((mode (SF)));
  93 typedef float DFtype __attribute__ ((mode (DF)));
  94 #ifdef TMODES
  95 typedef float TFtype __attribute__ ((mode (TF)));
  96 #endif
  97
  98 typedef int HItype __attribute__ ((mode (HI)));
  99 typedef int SItype __attribute__ ((mode (SI)));
 100 typedef int DItype __attribute__ ((mode (DI)));
 101 #ifdef TMODES
 102 typedef int TItype __attribute__ ((mode (TI)));
 103 #endif
 104
 105 /* The type of the result of a floating point comparison.  This must
 106    match `__libgcc_cmp_return__' in GCC for the target.  */
 107 #ifndef CMPtype
 108 typedef int CMPtype __attribute__ ((mode (__libgcc_cmp_return__)));
 109 #endif
 110
 111 typedef unsigned int UHItype __attribute__ ((mode (HI)));
 112 typedef unsigned int USItype __attribute__ ((mode (SI)));
 113 typedef unsigned int UDItype __attribute__ ((mode (DI)));
 114 #ifdef TMODES
 115 typedef unsigned int UTItype __attribute__ ((mode (TI)));
 116 #endif
 117
 118 #define MAX_USI_INT  (~(USItype)0)
 119 #define MAX_SI_INT   ((SItype) (MAX_USI_INT >> 1))
 120 #define BITS_PER_SI  (4 * __CHAR_BIT__)
 121 #ifdef TMODES
 122 #define MAX_UDI_INT  (~(UDItype)0)
 123 #define MAX_DI_INT   ((DItype) (MAX_UDI_INT >> 1))
 124 #define BITS_PER_DI  (8 * __CHAR_BIT__)
 125 #endif
 126
 127 #ifdef FLOAT_ONLY
 128 #define NO_DI_MODE
 129 #endif
 130
 131 #if __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__
 132 #define FLOAT_BIT_ORDER_MISMATCH
 133 #endif
 134
 135 #if __BYTE_ORDER__ != __FLOAT_WORD_ORDER__
 136 #define FLOAT_WORD_ORDER_MISMATCH
 137 #endif
 138
 139 #ifdef TFLOAT
 140 # ifndef TMODES
 141 #  error "TFLOAT requires long double to have 113 bits of mantissa"
 142 # endif
 143
 144 #       define PREFIXFPDP tp
 145 #       define PREFIXSFDF tf
 146 #       define NGARDS 10L /* Is this right? */
 147 #       define GARDROUND 0x1ff
 148 #       define GARDMASK  0x3ff
 149 #       define GARDMSB   0x200
 150 #       define FRAC_NBITS 128
 151
 152 # if __LDBL_MANT_DIG__ == 113 /* IEEE quad */
 153 #       define EXPBITS 15
 154 #       define EXPBIAS 16383
 155 #       define EXPMAX (0x7fff)
 156 #       define QUIET_NAN ((TItype)0x8 << 108)
 157 #       define FRACHIGH  ((TItype)0x8 << 124)
 158 #       define FRACHIGH2 ((TItype)0xc << 124)
 159 #       define FRACBITS 112
 160 # endif
 161
 162 # if __LDBL_MANT_DIG__ == 106 /* IBM extended (double+double) */
 163 #       define EXPBITS 11
 164 #       define EXPBIAS 1023
 165 #       define EXPMAX (0x7ff)
 166 #       define QUIET_NAN ((TItype)0x8 << (48 + 64))
 167 #       define FRACHIGH  ((TItype)0x8 << 124)
 168 #       define FRACHIGH2 ((TItype)0xc << 124)
 169 #       define FRACBITS 105
 170 #       define HALFFRACBITS 52
 171 #       define HALFSHIFT 64
 172 # endif
 173
 174 #       define pack_d __pack_t
 175 #       define unpack_d __unpack_t
 176 #       define __fpcmp_parts __fpcmp_parts_t
 177         typedef UTItype fractype;
 178         typedef UDItype halffractype;
 179         typedef USItype qrtrfractype;
 180 #define qrtrfractype qrtrfractype
 181         typedef TFtype FLO_type;
 182         typedef TItype intfrac;
 183 #elif defined FLOAT
 184 #       define NGARDS    7L
 185 #       define GARDROUND 0x3f
 186 #       define GARDMASK  0x7f
 187 #       define GARDMSB   0x40
 188 #       define EXPBITS 8
 189 #       define EXPBIAS 127
 190 #       define FRACBITS 23
 191 #       define EXPMAX (0xff)
 192 #       define QUIET_NAN 0x400000L
 193 #       define FRAC_NBITS 32
 194 #       define FRACHIGH  0x80000000L
 195 #       define FRACHIGH2 0xc0000000L
 196 #       define pack_d __pack_f
 197 #       define unpack_d __unpack_f
 198 #       define __fpcmp_parts __fpcmp_parts_f
 199         typedef USItype fractype;
 200         typedef UHItype halffractype;
 201         typedef SFtype FLO_type;
 202         typedef SItype intfrac;
 203
 204 #else
 205 #       define PREFIXFPDP dp
 206 #       define PREFIXSFDF df
 207 #       define NGARDS 8L
 208 #       define GARDROUND 0x7f
 209 #       define GARDMASK  0xff
 210 #       define GARDMSB   0x80
 211 #       define EXPBITS 11
 212 #       define EXPBIAS 1023
 213 #       define FRACBITS 52
 214 #       define EXPMAX (0x7ff)
 215 #       define QUIET_NAN 0x8000000000000LL
 216 #       define FRAC_NBITS 64
 217 #       define FRACHIGH  0x8000000000000000LL
 218 #       define FRACHIGH2 0xc000000000000000LL
 219 #       define pack_d __pack_d
 220 #       define unpack_d __unpack_d
 221 #       define __fpcmp_parts __fpcmp_parts_d
 222         typedef UDItype fractype;
 223         typedef USItype halffractype;
 224         typedef DFtype FLO_type;
 225         typedef DItype intfrac;
 226 #endif /* FLOAT */
 227
 228 #ifdef TFLOAT
 229 #       define add              __addtf3
 230 #       define sub              __subtf3
 231 #       define multiply         __multf3
 232 #       define divide           __divtf3
 233 #       define compare          __cmptf2
 234 #       define _eq_f2           __eqtf2
 235 #       define _ne_f2           __netf2
 236 #       define _gt_f2           __gttf2
 237 #       define _ge_f2           __getf2
 238 #       define _lt_f2           __lttf2
 239 #       define _le_f2           __letf2
 240 #       define _unord_f2        __unordtf2
 241 #       define usi_to_float     __floatunsitf
 242 #       define si_to_float      __floatsitf
 243 #       define float_to_si      __fixtfsi
 244 #       define float_to_usi     __fixunstfsi
 245 #       define negate           __negtf2
 246 #       define tf_to_sf         __trunctfsf2
 247 #       define tf_to_df         __trunctfdf2
 248 #elif defined FLOAT
 249 #       define add              __addsf3
 250 #       define sub              __subsf3
 251 #       define multiply         __mulsf3
 252 #       define divide           __divsf3
 253 #       define compare          __cmpsf2
 254 #       define _eq_f2           __eqsf2
 255 #       define _ne_f2           __nesf2
 256 #       define _gt_f2           __gtsf2
 257 #       define _ge_f2           __gesf2
 258 #       define _lt_f2           __ltsf2
 259 #       define _le_f2           __lesf2
 260 #       define _unord_f2        __unordsf2
 261 #       define usi_to_float     __floatunsisf
 262 #       define si_to_float      __floatsisf
 263 #       define float_to_si      __fixsfsi
 264 #       define float_to_usi     __fixunssfsi
 265 #       define negate           __negsf2
 266 #       define sf_to_df         __extendsfdf2
 267 #       define sf_to_tf         __extendsftf2
 268 #else
 269 #       define add              __adddf3
 270 #       define sub              __subdf3
 271 #       define multiply         __muldf3
 272 #       define divide           __divdf3
 273 #       define compare          __cmpdf2
 274 #       define _eq_f2           __eqdf2
 275 #       define _ne_f2           __nedf2
 276 #       define _gt_f2           __gtdf2
 277 #       define _ge_f2           __gedf2
 278 #       define _lt_f2           __ltdf2
 279 #       define _le_f2           __ledf2
 280 #       define _unord_f2        __unorddf2
 281 #       define usi_to_float     __floatunsidf
 282 #       define si_to_float      __floatsidf
 283 #       define float_to_si      __fixdfsi
 284 #       define float_to_usi     __fixunsdfsi
 285 #       define negate           __negdf2
 286 #       define df_to_sf         __truncdfsf2
 287 #       define df_to_tf         __extenddftf2
 288 #endif /* FLOAT */
 289
 290 #ifndef INLINE
 291 #define INLINE __inline__
 292 #endif
 293
 294 /* Preserve the sticky-bit when shifting fractions to the right.  */
 295 #define LSHIFT(a, s) { a = (a >> s) | !!(a & (((fractype) 1 << s) - 1)); }
 296
 297 /* numeric parameters */
 298 /* F_D_BITOFF is the number of bits offset between the MSB of the mantissa
 299    of a float and of a double. Assumes there are only two float types.
 300    (double::FRAC_BITS+double::NGARDS-(float::FRAC_BITS+float::NGARDS))
 301  */
 302 #define F_D_BITOFF (52+8-(23+7))
 303
 304 #ifdef TMODES
 305 # define F_T_BITOFF (__LDBL_MANT_DIG__-1+10-(23+7))
 306 # define D_T_BITOFF (__LDBL_MANT_DIG__-1+10-(52+8))
 307 #endif
 308
 309
 310 #define NORMAL_EXPMIN (-(EXPBIAS)+1)
 311 #define IMPLICIT_1 ((fractype)1<<(FRACBITS+NGARDS))
 312 #define IMPLICIT_2 ((fractype)1<<(FRACBITS+1+NGARDS))
 313
 314 /* common types */
 315
 316 typedef enum
 317 {
 318   CLASS_SNAN,
 319   CLASS_QNAN,
 320   CLASS_ZERO,
 321   CLASS_NUMBER,
 322   CLASS_INFINITY
 323 } fp_class_type;
 324
 325 typedef struct
 326 {
 327 #ifdef SMALL_MACHINE
 328   char class;
 329   unsigned char sign;
 330   short normal_exp;
 331 #else
 332   fp_class_type class;
 333   unsigned int sign;
 334   int normal_exp;
 335 #endif
 336
 337   union
 338     {
 339       fractype ll;
 340       halffractype l[2];
 341     } fraction;
 342 } fp_number_type;
 343
 344 typedef union
 345 {
 346   FLO_type value;
 347   fractype value_raw;
 348
 349 #ifndef FLOAT
 350 # ifdef qrtrfractype
 351   qrtrfractype qwords[4];
 352 # else
 353   halffractype words[2];
 354 # endif
 355 #endif
 356
 357 #ifdef FLOAT_BIT_ORDER_MISMATCH
 358   struct
 359     {
 360       fractype fraction:FRACBITS __attribute__ ((packed));
 361       unsigned int exp:EXPBITS __attribute__ ((packed));
 362       unsigned int sign:1 __attribute__ ((packed));
 363     }
 364   bits;
 365 #endif
 366
 367 #ifdef _DEBUG_BITFLOAT
 368   struct
 369     {
 370       unsigned int sign:1 __attribute__ ((packed));
 371       unsigned int exp:EXPBITS __attribute__ ((packed));
 372       fractype fraction:FRACBITS __attribute__ ((packed));
 373     }
 374   bits_big_endian;
 375
 376   struct
 377     {
 378       fractype fraction:FRACBITS __attribute__ ((packed));
 379       unsigned int exp:EXPBITS __attribute__ ((packed));
 380       unsigned int sign:1 __attribute__ ((packed));
 381     }
 382   bits_little_endian;
 383 #endif
 384 }
 385 FLO_union_type;
 386
 387 /* Prototypes.  */
 388
 389 #if defined(L_pack_df) || defined(L_pack_sf) || defined(L_pack_tf)
 390 extern FLO_type pack_d (const fp_number_type *);
 391 #endif
 392
 393 extern void unpack_d (FLO_union_type *, fp_number_type *);
 394
 395 #if defined(L_addsub_sf) || defined(L_addsub_df) || defined(L_addsub_tf)
 396 extern FLO_type add (FLO_type, FLO_type);
 397 extern FLO_type sub (FLO_type, FLO_type);
 398 #endif
 399
 400 #if defined(L_mul_sf) || defined(L_mul_df) || defined(L_mul_tf)
 401 extern FLO_type multiply (FLO_type, FLO_type);
 402 #endif
 403
 404 #if defined(L_div_sf) || defined(L_div_df) || defined(L_div_tf)
 405 extern FLO_type divide (FLO_type, FLO_type);
 406 #endif
 407
 408 extern int __fpcmp_parts (fp_number_type *, fp_number_type *);
 409
 410 #if defined(L_compare_sf) || defined(L_compare_df) || defined(L_compare_tf)
 411 extern CMPtype compare (FLO_type, FLO_type);
 412 #endif
 413
 414 #if defined(L_eq_sf) || defined(L_eq_df) || defined(L_eq_tf)
 415 extern CMPtype _eq_f2 (FLO_type, FLO_type);
 416 #endif
 417
 418 #if defined(L_ne_sf) || defined(L_ne_df) || defined(L_ne_tf)
 419 extern CMPtype _ne_f2 (FLO_type, FLO_type);
 420 #endif
 421
 422 #if defined(L_gt_sf) || defined(L_gt_df) || defined(L_gt_tf)
 423 extern CMPtype _gt_f2 (FLO_type, FLO_type);
 424 #endif
 425
 426 #if defined(L_ge_sf) || defined(L_ge_df) || defined(L_ge_tf)
 427 extern CMPtype _ge_f2 (FLO_type, FLO_type);
 428 #endif
 429
 430 #if defined(L_lt_sf) || defined(L_lt_df) || defined(L_lt_tf)
 431 extern CMPtype _lt_f2 (FLO_type, FLO_type);
 432 #endif
 433
 434 #if defined(L_le_sf) || defined(L_le_df) || defined(L_le_tf)
 435 extern CMPtype _le_f2 (FLO_type, FLO_type);
 436 #endif
 437
 438 #if defined(L_unord_sf) || defined(L_unord_df) || defined(L_unord_tf)
 439 extern CMPtype _unord_f2 (FLO_type, FLO_type);
 440 #endif
 441
 442 #if defined(L_si_to_sf) || defined(L_si_to_df) || defined(L_si_to_tf)
 443 extern FLO_type si_to_float (SItype);
 444 #endif
 445
 446 #if defined(L_sf_to_si) || defined(L_df_to_si) || defined(L_tf_to_si)
 447 extern SItype float_to_si (FLO_type);
 448 #endif
 449
 450 #if defined(L_tf_to_usi)
 451 extern USItype float_to_usi (FLO_type);
 452 #endif
 453
 454 #if defined(L_usi_to_sf) || defined(L_usi_to_df) || defined(L_usi_to_tf)
 455 extern FLO_type usi_to_float (USItype);
 456 #endif
 457
 458 #if defined(L_negate_sf) || defined(L_negate_df) || defined(L_negate_tf)
 459 extern FLO_type negate (FLO_type);
 460 #endif
 461
 462 #ifdef FLOAT
 463 #if defined(L_make_sf)
 464 extern SFtype __make_fp (fp_class_type, unsigned int, int, USItype);
 465 #endif
 466 #ifndef FLOAT_ONLY
 467 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
 468 #if defined(L_sf_to_df)
 469 extern DFtype sf_to_df (SFtype);
 470 #endif
 471 #if defined(L_sf_to_tf) && defined(TMODES)
 472 extern TFtype sf_to_tf (SFtype);
 473 #endif
 474 #endif /* ! FLOAT_ONLY */
 475 #endif /* FLOAT */
 476
 477 #ifndef FLOAT
 478 extern SFtype __make_fp (fp_class_type, unsigned int, int, USItype);
 479 #if defined(L_make_df)
 480 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
 481 #endif
 482 #if defined(L_df_to_sf)
 483 extern SFtype df_to_sf (DFtype);
 484 #endif
 485 #if defined(L_df_to_tf) && defined(TMODES)
 486 extern TFtype df_to_tf (DFtype);
 487 #endif
 488 #endif /* ! FLOAT */
 489
 490 #ifdef TMODES
 491 extern DFtype __make_dp (fp_class_type, unsigned int, int, UDItype);
 492 extern TFtype __make_tp (fp_class_type, unsigned int, int, UTItype);
 493 #ifdef TFLOAT
 494 #if defined(L_tf_to_sf)
 495 extern SFtype tf_to_sf (TFtype);
 496 #endif
 497 #if defined(L_tf_to_df)
 498 extern DFtype tf_to_df (TFtype);
 499 #endif
 500 #if defined(L_di_to_tf)
 501 extern TFtype di_to_df (DItype);
 502 #endif
 503 #endif /* TFLOAT */
 504 #endif /* TMODES */
 505
 506 #endif /* ! GCC_FP_BIT_H */