1#ifndef _IEEE_FP_H_ 2#define _IEEE_FP_H_ 3 4#include "_ansi.h" 5 6#include <machine/ieeefp.h> 7 8_BEGIN_STD_C 9 10/* FIXME FIXME FIXME: 11 Neither of __ieee_{float,double}_shape_tape seem to be used anywhere 12 except in libm/test. If that is the case, please delete these from here. 13 If that is not the case, please insert documentation here describing why 14 they're needed. */ 15 16#ifdef __IEEE_BIG_ENDIAN 17 18typedef union 19{ 20 double value; 21 struct 22 { 23 unsigned int sign : 1; 24 unsigned int exponent: 11; 25 unsigned int fraction0:4; 26 unsigned int fraction1:16; 27 unsigned int fraction2:16; 28 unsigned int fraction3:16; 29 30 } number; 31 struct 32 { 33 unsigned int sign : 1; 34 unsigned int exponent: 11; 35 unsigned int quiet:1; 36 unsigned int function0:3; 37 unsigned int function1:16; 38 unsigned int function2:16; 39 unsigned int function3:16; 40 } nan; 41 struct 42 { 43 unsigned long msw; 44 unsigned long lsw; 45 } parts; 46 long aslong[2]; 47} __ieee_double_shape_type; 48 49#endif 50 51#ifdef __IEEE_LITTLE_ENDIAN 52 53typedef union 54{ 55 double value; 56 struct 57 { 58#ifdef __SMALL_BITFIELDS 59 unsigned int fraction3:16; 60 unsigned int fraction2:16; 61 unsigned int fraction1:16; 62 unsigned int fraction0: 4; 63#else 64 unsigned int fraction1:32; 65 unsigned int fraction0:20; 66#endif 67 unsigned int exponent :11; 68 unsigned int sign : 1; 69 } number; 70 struct 71 { 72#ifdef __SMALL_BITFIELDS 73 unsigned int function3:16; 74 unsigned int function2:16; 75 unsigned int function1:16; 76 unsigned int function0:3; 77#else 78 unsigned int function1:32; 79 unsigned int function0:19; 80#endif 81 unsigned int quiet:1; 82 unsigned int exponent: 11; 83 unsigned int sign : 1; 84 } nan; 85 struct 86 { 87 unsigned long lsw; 88 unsigned long msw; 89 } parts; 90 91 long aslong[2]; 92 93} __ieee_double_shape_type; 94 95#endif 96 97#ifdef __IEEE_BIG_ENDIAN 98 99typedef union 100{ 101 float value; 102 struct 103 { 104 unsigned int sign : 1; 105 unsigned int exponent: 8; 106 unsigned int fraction0: 7; 107 unsigned int fraction1: 16; 108 } number; 109 struct 110 { 111 unsigned int sign:1; 112 unsigned int exponent:8; 113 unsigned int quiet:1; 114 unsigned int function0:6; 115 unsigned int function1:16; 116 } nan; 117 long p1; 118 119} __ieee_float_shape_type; 120 121#endif 122 123#ifdef __IEEE_LITTLE_ENDIAN 124 125typedef union 126{ 127 float value; 128 struct 129 { 130 unsigned int fraction0: 7; 131 unsigned int fraction1: 16; 132 unsigned int exponent: 8; 133 unsigned int sign : 1; 134 } number; 135 struct 136 { 137 unsigned int function1:16; 138 unsigned int function0:6; 139 unsigned int quiet:1; 140 unsigned int exponent:8; 141 unsigned int sign:1; 142 } nan; 143 long p1; 144 145} __ieee_float_shape_type; 146 147#endif 148 149 150 151 152 153/* FLOATING ROUNDING */ 154 155typedef int fp_rnd; 156#define FP_RN 0 /* Round to nearest */ 157#define FP_RM 1 /* Round down */ 158#define FP_RP 2 /* Round up */ 159#define FP_RZ 3 /* Round to zero (trunate) */ 160 161fp_rnd _EXFUN(fpgetround,(void)); 162fp_rnd _EXFUN(fpsetround, (fp_rnd)); 163 164/* EXCEPTIONS */ 165 166typedef int fp_except; 167#define FP_X_INV 0x10 /* Invalid operation */ 168#define FP_X_DX 0x80 /* Divide by zero */ 169#define FP_X_OFL 0x04 /* Overflow exception */ 170#define FP_X_UFL 0x02 /* Underflow exception */ 171#define FP_X_IMP 0x01 /* imprecise exception */ 172 173fp_except _EXFUN(fpgetmask,(void)); 174fp_except _EXFUN(fpsetmask,(fp_except)); 175fp_except _EXFUN(fpgetsticky,(void)); 176fp_except _EXFUN(fpsetsticky, (fp_except)); 177 178/* INTEGER ROUNDING */ 179 180typedef int fp_rdi; 181#define FP_RDI_TOZ 0 /* Round to Zero */ 182#define FP_RDI_RD 1 /* Follow float mode */ 183 184fp_rdi _EXFUN(fpgetroundtoi,(void)); 185fp_rdi _EXFUN(fpsetroundtoi,(fp_rdi)); 186 187#undef isnan 188#undef isinf 189 190int _EXFUN(isnan, (double)); 191int _EXFUN(isinf, (double)); 192int _EXFUN(finite, (double)); 193 194 195 196int _EXFUN(isnanf, (float)); 197int _EXFUN(isinff, (float)); 198int _EXFUN(finitef, (float)); 199 200#define __IEEE_DBL_EXPBIAS 1023 201#define __IEEE_FLT_EXPBIAS 127 202 203#define __IEEE_DBL_EXPLEN 11 204#define __IEEE_FLT_EXPLEN 8 205 206 207#define __IEEE_DBL_FRACLEN (64 - (__IEEE_DBL_EXPLEN + 1)) 208#define __IEEE_FLT_FRACLEN (32 - (__IEEE_FLT_EXPLEN + 1)) 209 210#define __IEEE_DBL_MAXPOWTWO ((double)(1L << 32 - 2) * (1L << (32-11) - 32 + 1)) 211#define __IEEE_FLT_MAXPOWTWO ((float)(1L << (32-8) - 1)) 212 213#define __IEEE_DBL_NAN_EXP 0x7ff 214#define __IEEE_FLT_NAN_EXP 0xff 215 216#ifndef __ieeefp_isnanf 217#define __ieeefp_isnanf(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \ 218 ((*(long *)&(x) & 0x007fffffL)!=0000000000L)) 219#endif 220#define isnanf(x) __ieeefp_isnanf(x) 221 222#ifndef __ieeefp_isinff 223#define __ieeefp_isinff(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \ 224 ((*(long *)&(x) & 0x007fffffL)==0000000000L)) 225#endif 226#define isinff(x) __ieeefp_isinff(x) 227 228#ifndef __ieeefp_finitef 229#define __ieeefp_finitef(x) (((*(long *)&(x) & 0x7f800000L)!=0x7f800000L)) 230#endif 231#define finitef(x) __ieeefp_finitef(x) 232 233#ifdef _DOUBLE_IS_32BITS 234#undef __IEEE_DBL_EXPBIAS 235#define __IEEE_DBL_EXPBIAS __IEEE_FLT_EXPBIAS 236 237#undef __IEEE_DBL_EXPLEN 238#define __IEEE_DBL_EXPLEN __IEEE_FLT_EXPLEN 239 240#undef __IEEE_DBL_FRACLEN 241#define __IEEE_DBL_FRACLEN __IEEE_FLT_FRACLEN 242 243#undef __IEEE_DBL_MAXPOWTWO 244#define __IEEE_DBL_MAXPOWTWO __IEEE_FLT_MAXPOWTWO 245 246#undef __IEEE_DBL_NAN_EXP 247#define __IEEE_DBL_NAN_EXP __IEEE_FLT_NAN_EXP 248 249#undef __ieee_double_shape_type 250#define __ieee_double_shape_type __ieee_float_shape_type 251 252#endif /* _DOUBLE_IS_32BITS */ 253 254_END_STD_C 255 256#endif /* _IEEE_FP_H_ */ 257