nsz Git - libm/blob - src/internal/ldhack.h

   1 #ifndef _LDHACK_H
   2 #define _LDHACK_H
   3
   4 #include <float.h>
   5 #include <stdint.h>
   6
   7 // FIXME: conveniance macros
   8
   9 #undef LD64
  10 #undef LD80
  11 #undef LD128
  12 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
  13 #define LD64 1
  14 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
  15 #define LD80 1
  16 union ldshape {
  17         long double value;
  18         struct {
  19                 uint64_t m;
  20                 uint16_t exp:15;
  21                 uint16_t sign:1;
  22                 uint16_t pad;
  23         } bits;
  24 };
  25 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
  26 #define LD128 1
  27 union ldshape {
  28         long double value;
  29         struct {
  30                 uint64_t mlo;
  31                 uint64_t mhi:48;
  32                 uint16_t exp:15;
  33                 uint16_t sign:1;
  34         } bits;
  35 };
  36 #else
  37 #error Unsupported long double representation
  38 #endif
  39
  40
  41 // FIXME: hacks to make freebsd+openbsd long double code happy
  42
  43 // union and macros for freebsd
  44
  45 #if LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
  46
  47 union IEEEl2bits {
  48         long double e;
  49         struct {
  50                 uint32_t manl:32;
  51                 uint32_t manh:32;
  52                 uint32_t exp:15;
  53                 uint32_t sign:1;
  54                 uint32_t pad:16;
  55         } bits;
  56         struct {
  57                 uint64_t man:64;
  58                 uint32_t expsign:16;
  59                 uint32_t pad:16;
  60         } xbits;
  61 };
  62
  63 #define LDBL_MANL_SIZE 32
  64 #define LDBL_MANH_SIZE 32
  65 #define LDBL_NBIT (1ull << LDBL_MANH_SIZE-1)
  66 #undef LDBL_IMPLICIT_NBIT
  67 #define mask_nbit_l(u) ((u).bits.manh &= ~LDBL_NBIT)
  68
  69 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
  70 /*
  71 // ld128 float.h
  72 //#define LDBL_MAX 1.189731495357231765085759326628007016E+4932L
  73 #define LDBL_MAX 0x1.ffffffffffffffffffffffffffffp+16383
  74 #define LDBL_MAX_EXP 16384
  75 #define LDBL_HAS_INFINITY 1
  76 //#define LDBL_MIN 3.362103143112093506262677817321752603E-4932L
  77 #define LDBL_MIN 0x1p-16382
  78 #define LDBL_HAS_QUIET_NAN 1
  79 #define LDBL_HAS_DENORM 1
  80 //#define LDBL_EPSILON 1.925929944387235853055977942584927319E-34L
  81 #define LDBL_EPSILON 0x1p-112
  82 #define LDBL_MANT_DIG 113
  83 #define LDBL_MIN_EXP (-16381)
  84 #define LDBL_MAX_10_EXP 4932
  85 #define LDBL_DENORM_MIN 0x0.0000000000000000000000000001p-16381
  86 #define LDBL_MIN_10_EXP (-4931)
  87 #define LDBL_DIG 33
  88 */
  89
  90 union IEEEl2bits {
  91         long double e;
  92         struct {
  93                 uint64_t manl:64;
  94                 uint64_t manh:48;
  95                 uint32_t exp:15;
  96                 uint32_t sign:1;
  97         } bits;
  98         struct {
  99                 uint64_t unused0:64;
 100                 uint64_t unused1:48;
 101                 uint32_t expsign:16;
 102         } xbits;
 103 };
 104
 105 #define LDBL_MANL_SIZE 64
 106 #define LDBL_MANH_SIZE 48
 107 #define LDBL_NBIT (1ull << LDBL_MANH_SIZE)
 108 #define LDBL_IMPLICIT_NBIT 1
 109 #define mask_nbit_l(u)
 110
 111 #endif
 112
 113
 114 // macros for openbsd
 115
 116 #define GET_LDOUBLE_WORDS(se,mh,ml, f) do{ \
 117         union IEEEl2bits u; \
 118         u.e = (f); \
 119         (se) = u.xbits.expsign; \
 120         (mh) = u.bits.manh; \
 121         (ml) = u.bits.manl; \
 122 }while(0)
 123
 124 #define SET_LDOUBLE_WORDS(f,  se,mh,ml) do{ \
 125         union IEEEl2bits u; \
 126         u.xbits.expsign = (se); \
 127         u.bits.manh = (mh); \
 128         u.bits.manl = (ml); \
 129         (f) = u.e; \
 130 }while(0)
 131
 132 #define GET_LDOUBLE_EXP(se, f) do{ \
 133         union IEEEl2bits u; \
 134         u.e = (f); \
 135         (se) = u.xbits.expsign; \
 136 }while(0)
 137
 138 #define SET_LDOUBLE_EXP(f, se) do{ \
 139         union IEEEl2bits u; \
 140         u.e = (f); \
 141         u.xbits.expsign = (se); \
 142         (f) = u.e; \
 143 }while(0)
 144
 145 #endif