nkeynes@185: #ifndef _IEEE_FP_H_
nkeynes@185: #define _IEEE_FP_H_
nkeynes@185: 
nkeynes@185: #include "_ansi.h"
nkeynes@185: 
nkeynes@185: #include <machine/ieeefp.h>
nkeynes@185: 
nkeynes@185: /* FIXME FIXME FIXME:
nkeynes@185:    Neither of __ieee_{float,double}_shape_tape seem to be used anywhere
nkeynes@185:    except in libm/test.  If that is the case, please delete these from here.
nkeynes@185:    If that is not the case, please insert documentation here describing why
nkeynes@185:    they're needed.  */
nkeynes@185: 
nkeynes@185: #ifdef __IEEE_BIG_ENDIAN
nkeynes@185: 
nkeynes@185: typedef union 
nkeynes@185: {
nkeynes@185:   double value;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int sign : 1;
nkeynes@185:     unsigned int exponent: 11;
nkeynes@185:     unsigned int fraction0:4;
nkeynes@185:     unsigned int fraction1:16;
nkeynes@185:     unsigned int fraction2:16;
nkeynes@185:     unsigned int fraction3:16;
nkeynes@185:     
nkeynes@185:   } number;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int sign : 1;
nkeynes@185:     unsigned int exponent: 11;
nkeynes@185:     unsigned int quiet:1;
nkeynes@185:     unsigned int function0:3;
nkeynes@185:     unsigned int function1:16;
nkeynes@185:     unsigned int function2:16;
nkeynes@185:     unsigned int function3:16;
nkeynes@185:   } nan;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned long msw;
nkeynes@185:     unsigned long lsw;
nkeynes@185:   } parts;
nkeynes@185:     long aslong[2];
nkeynes@185: } __ieee_double_shape_type;
nkeynes@185: 
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: #ifdef __IEEE_LITTLE_ENDIAN
nkeynes@185: 
nkeynes@185: typedef union 
nkeynes@185: {
nkeynes@185:   double value;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185: #ifdef __SMALL_BITFIELDS
nkeynes@185:     unsigned int fraction3:16;
nkeynes@185:     unsigned int fraction2:16;
nkeynes@185:     unsigned int fraction1:16;
nkeynes@185:     unsigned int fraction0: 4;
nkeynes@185: #else
nkeynes@185:     unsigned int fraction1:32;
nkeynes@185:     unsigned int fraction0:20;
nkeynes@185: #endif
nkeynes@185:     unsigned int exponent :11;
nkeynes@185:     unsigned int sign     : 1;
nkeynes@185:   } number;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185: #ifdef __SMALL_BITFIELDS
nkeynes@185:     unsigned int function3:16;
nkeynes@185:     unsigned int function2:16;
nkeynes@185:     unsigned int function1:16;
nkeynes@185:     unsigned int function0:3;
nkeynes@185: #else
nkeynes@185:     unsigned int function1:32;
nkeynes@185:     unsigned int function0:19;
nkeynes@185: #endif
nkeynes@185:     unsigned int quiet:1;
nkeynes@185:     unsigned int exponent: 11;
nkeynes@185:     unsigned int sign : 1;
nkeynes@185:   } nan;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned long lsw;
nkeynes@185:     unsigned long msw;
nkeynes@185:   } parts;
nkeynes@185: 
nkeynes@185:   long aslong[2];
nkeynes@185: 
nkeynes@185: } __ieee_double_shape_type;
nkeynes@185: 
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: #ifdef __IEEE_BIG_ENDIAN
nkeynes@185: 
nkeynes@185: typedef union
nkeynes@185: {
nkeynes@185:   float value;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int sign : 1;
nkeynes@185:     unsigned int exponent: 8;
nkeynes@185:     unsigned int fraction0: 7;
nkeynes@185:     unsigned int fraction1: 16;
nkeynes@185:   } number;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int sign:1;
nkeynes@185:     unsigned int exponent:8;
nkeynes@185:     unsigned int quiet:1;
nkeynes@185:     unsigned int function0:6;
nkeynes@185:     unsigned int function1:16;
nkeynes@185:   } nan;
nkeynes@185:   long p1;
nkeynes@185:   
nkeynes@185: } __ieee_float_shape_type;
nkeynes@185: 
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: #ifdef __IEEE_LITTLE_ENDIAN
nkeynes@185: 
nkeynes@185: typedef union
nkeynes@185: {
nkeynes@185:   float value;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int fraction0: 7;
nkeynes@185:     unsigned int fraction1: 16;
nkeynes@185:     unsigned int exponent: 8;
nkeynes@185:     unsigned int sign : 1;
nkeynes@185:   } number;
nkeynes@185:   struct 
nkeynes@185:   {
nkeynes@185:     unsigned int function1:16;
nkeynes@185:     unsigned int function0:6;
nkeynes@185:     unsigned int quiet:1;
nkeynes@185:     unsigned int exponent:8;
nkeynes@185:     unsigned int sign:1;
nkeynes@185:   } nan;
nkeynes@185:   long p1;
nkeynes@185:   
nkeynes@185: } __ieee_float_shape_type;
nkeynes@185: 
nkeynes@185: #endif
nkeynes@185: 
nkeynes@185: 
nkeynes@185: 
nkeynes@185: 
nkeynes@185: 
nkeynes@185: /* FLOATING ROUNDING */
nkeynes@185: 
nkeynes@185: typedef int fp_rnd;
nkeynes@185: #define FP_RN 0 	/* Round to nearest 		*/
nkeynes@185: #define FP_RM 1		/* Round down 			*/
nkeynes@185: #define FP_RP 2		/* Round up 			*/
nkeynes@185: #define FP_RZ 3		/* Round to zero (trunate) 	*/
nkeynes@185: 
nkeynes@185: fp_rnd _EXFUN(fpgetround,(void));
nkeynes@185: fp_rnd _EXFUN(fpsetround, (fp_rnd));
nkeynes@185: 
nkeynes@185: /* EXCEPTIONS */
nkeynes@185: 
nkeynes@185: typedef int fp_except;
nkeynes@185: #define FP_X_INV 0x10	/* Invalid operation 		*/
nkeynes@185: #define FP_X_DX  0x80	/* Divide by zero		*/
nkeynes@185: #define FP_X_OFL 0x04	/* Overflow exception		*/
nkeynes@185: #define FP_X_UFL 0x02	/* Underflow exception		*/
nkeynes@185: #define FP_X_IMP 0x01	/* imprecise exception		*/
nkeynes@185: 
nkeynes@185: fp_except _EXFUN(fpgetmask,(void));
nkeynes@185: fp_except _EXFUN(fpsetmask,(fp_except));
nkeynes@185: fp_except _EXFUN(fpgetsticky,(void));
nkeynes@185: fp_except _EXFUN(fpsetsticky, (fp_except));
nkeynes@185: 
nkeynes@185: /* INTEGER ROUNDING */
nkeynes@185: 
nkeynes@185: typedef int fp_rdi;
nkeynes@185: #define FP_RDI_TOZ 0	/* Round to Zero 		*/
nkeynes@185: #define FP_RDI_RD  1	/* Follow float mode		*/
nkeynes@185: 
nkeynes@185: fp_rdi _EXFUN(fpgetroundtoi,(void));
nkeynes@185: fp_rdi _EXFUN(fpsetroundtoi,(fp_rdi));
nkeynes@185: 
nkeynes@185: int _EXFUN(isnan, (double));
nkeynes@185: int _EXFUN(isinf, (double));
nkeynes@185: int _EXFUN(finite, (double));
nkeynes@185: 
nkeynes@185: 
nkeynes@185: 
nkeynes@185: int _EXFUN(isnanf, (float));
nkeynes@185: int _EXFUN(isinff, (float));
nkeynes@185: int _EXFUN(finitef, (float));
nkeynes@185: 
nkeynes@185: #define __IEEE_DBL_EXPBIAS 1023
nkeynes@185: #define __IEEE_FLT_EXPBIAS 127
nkeynes@185: 
nkeynes@185: #define __IEEE_DBL_EXPLEN 11
nkeynes@185: #define __IEEE_FLT_EXPLEN 8
nkeynes@185: 
nkeynes@185: 
nkeynes@185: #define __IEEE_DBL_FRACLEN (64 - (__IEEE_DBL_EXPLEN + 1))
nkeynes@185: #define __IEEE_FLT_FRACLEN (32 - (__IEEE_FLT_EXPLEN + 1))
nkeynes@185: 
nkeynes@185: #define __IEEE_DBL_MAXPOWTWO	((double)(1L << 32 - 2) * (1L << (32-11) - 32 + 1))
nkeynes@185: #define __IEEE_FLT_MAXPOWTWO	((float)(1L << (32-8) - 1))
nkeynes@185: 
nkeynes@185: #define __IEEE_DBL_NAN_EXP 0x7ff
nkeynes@185: #define __IEEE_FLT_NAN_EXP 0xff
nkeynes@185: 
nkeynes@185: 
nkeynes@185: #define isnanf(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \
nkeynes@185: 		   ((*(long *)&(x) & 0x007fffffL)!=0000000000L))
nkeynes@185: 
nkeynes@185: #define isinff(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \
nkeynes@185: 		   ((*(long *)&(x) & 0x007fffffL)==0000000000L))
nkeynes@185: 
nkeynes@185: #define finitef(x) (((*(long *)&(x) & 0x7f800000L)!=0x7f800000L))
nkeynes@185: 
nkeynes@185: #ifdef _DOUBLE_IS_32BITS
nkeynes@185: #undef __IEEE_DBL_EXPBIAS
nkeynes@185: #define __IEEE_DBL_EXPBIAS __IEEE_FLT_EXPBIAS
nkeynes@185: 
nkeynes@185: #undef __IEEE_DBL_EXPLEN
nkeynes@185: #define __IEEE_DBL_EXPLEN __IEEE_FLT_EXPLEN
nkeynes@185: 
nkeynes@185: #undef __IEEE_DBL_FRACLEN
nkeynes@185: #define __IEEE_DBL_FRACLEN __IEEE_FLT_FRACLEN
nkeynes@185: 
nkeynes@185: #undef __IEEE_DBL_MAXPOWTWO
nkeynes@185: #define __IEEE_DBL_MAXPOWTWO __IEEE_FLT_MAXPOWTWO
nkeynes@185: 
nkeynes@185: #undef __IEEE_DBL_NAN_EXP
nkeynes@185: #define __IEEE_DBL_NAN_EXP __IEEE_FLT_NAN_EXP
nkeynes@185: 
nkeynes@185: #undef __ieee_double_shape_type
nkeynes@185: #define __ieee_double_shape_type __ieee_float_shape_type
nkeynes@185: 
nkeynes@185: #endif /* _DOUBLE_IS_32BITS */
nkeynes@185: 
nkeynes@185: #endif /* _IEEE_FP_H_ */