123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295 |
- #ifndef _IEEE_FP_H_
- #define _IEEE_FP_H_
- #include "_ansi.h"
- #include <machine/ieeefp.h>
- #include <float.h>
- _BEGIN_STD_C
- /* FIXME FIXME FIXME:
- Neither of __ieee_{float,double}_shape_type seem to be used anywhere
- except in libm/test. If that is the case, please delete these from here.
- If that is not the case, please insert documentation here describing why
- they're needed. */
- #ifdef __IEEE_BIG_ENDIAN
- typedef union
- {
- double value;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 11;
- unsigned int fraction0:4;
- unsigned int fraction1:16;
- unsigned int fraction2:16;
- unsigned int fraction3:16;
-
- } number;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 11;
- unsigned int quiet:1;
- unsigned int function0:3;
- unsigned int function1:16;
- unsigned int function2:16;
- unsigned int function3:16;
- } nan;
- struct
- {
- unsigned long msw;
- unsigned long lsw;
- } parts;
- long aslong[2];
- } __ieee_double_shape_type;
- #elif defined __IEEE_LITTLE_ENDIAN
- typedef union
- {
- double value;
- struct
- {
- #ifdef __SMALL_BITFIELDS
- unsigned int fraction3:16;
- unsigned int fraction2:16;
- unsigned int fraction1:16;
- unsigned int fraction0: 4;
- #else
- unsigned int fraction1:32;
- unsigned int fraction0:20;
- #endif
- unsigned int exponent :11;
- unsigned int sign : 1;
- } number;
- struct
- {
- #ifdef __SMALL_BITFIELDS
- unsigned int function3:16;
- unsigned int function2:16;
- unsigned int function1:16;
- unsigned int function0:3;
- #else
- unsigned int function1:32;
- unsigned int function0:19;
- #endif
- unsigned int quiet:1;
- unsigned int exponent: 11;
- unsigned int sign : 1;
- } nan;
- struct
- {
- unsigned long lsw;
- unsigned long msw;
- } parts;
- long aslong[2];
- } __ieee_double_shape_type;
- #endif /* __IEEE_LITTLE_ENDIAN */
- #ifdef __IEEE_BIG_ENDIAN
- typedef union
- {
- float value;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 8;
- unsigned int fraction0: 7;
- unsigned int fraction1: 16;
- } number;
- struct
- {
- unsigned int sign:1;
- unsigned int exponent:8;
- unsigned int quiet:1;
- unsigned int function0:6;
- unsigned int function1:16;
- } nan;
- long p1;
-
- } __ieee_float_shape_type;
- #elif defined __IEEE_LITTLE_ENDIAN
- typedef union
- {
- float value;
- struct
- {
- unsigned int fraction0: 7;
- unsigned int fraction1: 16;
- unsigned int exponent: 8;
- unsigned int sign : 1;
- } number;
- struct
- {
- unsigned int function1:16;
- unsigned int function0:6;
- unsigned int quiet:1;
- unsigned int exponent:8;
- unsigned int sign:1;
- } nan;
- long p1;
-
- } __ieee_float_shape_type;
- #endif /* __IEEE_LITTLE_ENDIAN */
- #ifndef _LDBL_EQ_DBL
- #ifndef LDBL_MANT_DIG
- #error "LDBL_MANT_DIG not defined - should be found in float.h"
- #elif LDBL_MANT_DIG == DBL_MANT_DIG
- #error "double and long double are the same size but LDBL_EQ_DBL is not defined"
- #elif LDBL_MANT_DIG == 53
- /* This happens when doubles are 32-bits and long doubles are 64-bits. */
- #define EXT_EXPBITS 11
- #define EXT_FRACHBITS 20
- #define EXT_FRACLBITS 32
- #define __ieee_ext_field_type unsigned long
- #elif LDBL_MANT_DIG == 64
- #define EXT_EXPBITS 15
- #define EXT_FRACHBITS 32
- #define EXT_FRACLBITS 32
- #define __ieee_ext_field_type unsigned int
- #elif LDBL_MANT_DIG == 65
- #define EXT_EXPBITS 15
- #define EXT_FRACHBITS 32
- #define EXT_FRACLBITS 32
- #define __ieee_ext_field_type unsigned int
- #elif LDBL_MANT_DIG == 112
- #define EXT_EXPBITS 15
- #define EXT_FRACHBITS 48
- #define EXT_FRACLBITS 64
- #define __ieee_ext_field_type unsigned long long
- #elif LDBL_MANT_DIG == 113
- #define EXT_EXPBITS 15
- #define EXT_FRACHBITS 48
- #define EXT_FRACLBITS 64
- #define __ieee_ext_field_type unsigned long long
- #else
- #error Unsupported value for LDBL_MANT_DIG
- #endif
- #define EXT_EXP_INFNAN ((1 << EXT_EXPBITS) - 1) /* 32767 */
- #define EXT_EXP_BIAS ((1 << (EXT_EXPBITS - 1)) - 1) /* 16383 */
- #define EXT_FRACBITS (EXT_FRACLBITS + EXT_FRACHBITS)
- typedef struct ieee_ext
- {
- __ieee_ext_field_type ext_fracl : EXT_FRACLBITS;
- __ieee_ext_field_type ext_frach : EXT_FRACHBITS;
- __ieee_ext_field_type ext_exp : EXT_EXPBITS;
- __ieee_ext_field_type ext_sign : 1;
- } ieee_ext;
- typedef union ieee_ext_u
- {
- long double extu_ld;
- struct ieee_ext extu_ext;
- } ieee_ext_u;
- #endif /* ! _LDBL_EQ_DBL */
- /* FLOATING ROUNDING */
- typedef int fp_rnd;
- #define FP_RN 0 /* Round to nearest */
- #define FP_RM 1 /* Round down */
- #define FP_RP 2 /* Round up */
- #define FP_RZ 3 /* Round to zero (trunate) */
- fp_rnd fpgetround (void);
- fp_rnd fpsetround (fp_rnd);
- /* EXCEPTIONS */
- typedef int fp_except;
- #define FP_X_INV 0x10 /* Invalid operation */
- #define FP_X_DX 0x80 /* Divide by zero */
- #define FP_X_OFL 0x04 /* Overflow exception */
- #define FP_X_UFL 0x02 /* Underflow exception */
- #define FP_X_IMP 0x01 /* imprecise exception */
- fp_except fpgetmask (void);
- fp_except fpsetmask (fp_except);
- fp_except fpgetsticky (void);
- fp_except fpsetsticky (fp_except);
- /* INTEGER ROUNDING */
- typedef int fp_rdi;
- #define FP_RDI_TOZ 0 /* Round to Zero */
- #define FP_RDI_RD 1 /* Follow float mode */
- fp_rdi fpgetroundtoi (void);
- fp_rdi fpsetroundtoi (fp_rdi);
- #define __IEEE_DBL_EXPBIAS 1023
- #define __IEEE_FLT_EXPBIAS 127
- #define __IEEE_DBL_EXPLEN 11
- #define __IEEE_FLT_EXPLEN 8
- #define __IEEE_DBL_FRACLEN (64 - (__IEEE_DBL_EXPLEN + 1))
- #define __IEEE_FLT_FRACLEN (32 - (__IEEE_FLT_EXPLEN + 1))
- #define __IEEE_DBL_MAXPOWTWO ((double)(1L << 32 - 2) * (1L << (32-11) - 32 + 1))
- #define __IEEE_FLT_MAXPOWTWO ((float)(1L << (32-8) - 1))
- #define __IEEE_DBL_NAN_EXP 0x7ff
- #define __IEEE_FLT_NAN_EXP 0xff
- #ifdef __ieeefp_isnanf
- #define isnanf(x) __ieeefp_isnanf(x)
- #endif
- #ifdef __ieeefp_isinff
- #define isinff(x) __ieeefp_isinff(x)
- #endif
- #ifdef __ieeefp_finitef
- #define finitef(x) __ieeefp_finitef(x)
- #endif
- #ifdef _DOUBLE_IS_32BITS
- #undef __IEEE_DBL_EXPBIAS
- #define __IEEE_DBL_EXPBIAS __IEEE_FLT_EXPBIAS
- #undef __IEEE_DBL_EXPLEN
- #define __IEEE_DBL_EXPLEN __IEEE_FLT_EXPLEN
- #undef __IEEE_DBL_FRACLEN
- #define __IEEE_DBL_FRACLEN __IEEE_FLT_FRACLEN
- #undef __IEEE_DBL_MAXPOWTWO
- #define __IEEE_DBL_MAXPOWTWO __IEEE_FLT_MAXPOWTWO
- #undef __IEEE_DBL_NAN_EXP
- #define __IEEE_DBL_NAN_EXP __IEEE_FLT_NAN_EXP
- #undef __ieee_double_shape_type
- #define __ieee_double_shape_type __ieee_float_shape_type
- #endif /* _DOUBLE_IS_32BITS */
- _END_STD_C
- #endif /* _IEEE_FP_H_ */
|