0,0 → 1,256 |
#ifndef _IEEE_FP_H_ |
#define _IEEE_FP_H_ |
|
#include "_ansi.h" |
|
#include <machine/ieeefp.h> |
|
_BEGIN_STD_C |
|
/* FIXME FIXME FIXME: |
Neither of __ieee_{float,double}_shape_tape 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; |
|
#endif |
|
#ifdef __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 |
|
#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; |
|
#endif |
|
#ifdef __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 |
|
|
|
|
|
/* 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 _EXFUN(fpgetround,(void)); |
fp_rnd _EXFUN(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 _EXFUN(fpgetmask,(void)); |
fp_except _EXFUN(fpsetmask,(fp_except)); |
fp_except _EXFUN(fpgetsticky,(void)); |
fp_except _EXFUN(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 _EXFUN(fpgetroundtoi,(void)); |
fp_rdi _EXFUN(fpsetroundtoi,(fp_rdi)); |
|
#undef isnan |
#undef isinf |
|
int _EXFUN(isnan, (double)); |
int _EXFUN(isinf, (double)); |
int _EXFUN(finite, (double)); |
|
|
|
int _EXFUN(isnanf, (float)); |
int _EXFUN(isinff, (float)); |
int _EXFUN(finitef, (float)); |
|
#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 |
|
#ifndef __ieeefp_isnanf |
#define __ieeefp_isnanf(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \ |
((*(long *)&(x) & 0x007fffffL)!=0000000000L)) |
#endif |
#define isnanf(x) __ieeefp_isnanf(x) |
|
#ifndef __ieeefp_isinff |
#define __ieeefp_isinff(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \ |
((*(long *)&(x) & 0x007fffffL)==0000000000L)) |
#endif |
#define isinff(x) __ieeefp_isinff(x) |
|
#ifndef __ieeefp_finitef |
#define __ieeefp_finitef(x) (((*(long *)&(x) & 0x7f800000L)!=0x7f800000L)) |
#endif |
#define finitef(x) __ieeefp_finitef(x) |
|
#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_ */ |