/* * $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.4 2002/01/04 17:01:40 saito Exp $ */ #ifndef _INTERVAL_H #define _INTERVAL_H #define PRINTF_G 0 #define PRINTF_E 1 #if defined(INTERVAL) #include #ifdef sun #if OSMajorVersion < 4 || (OSMajorVersion == 4 && OSMinorVersion < 1) #include #endif #if defined(__svr4__) #include #define FPNEAREST fpsetround(FP_RN); #define FPPLUSINF fpsetround(FP_RP); #define FPMINUSINF fpsetround(FP_RM); #define FPTOZERO fpsetround(FP_RZ); #else static char *Interval_dummy; #define FPNEAREST ieee_flags("clear", "direction", Interval_dummy, &Interval_dummy); #define FPPLUSINF ieee_flags("set", "direction", "positive", &Interval_dummy); #define FPMINUSINF ieee_flags("set", "direction", "negative", &Interval_dummy); #define FPTOZERO ieee_flags("set", "direction", "tozero", &Interval_dummy); #endif #endif #ifdef linux #include #if 1 #define LINUX_FPU_RC_MASK 0xf3ff #define LINUX_FPU_SETCW(c) {_FPU_GETCW(__fpu_control);\ _FPU_SETCW(__fpu_control & LINUX_FPU_RC_MASK | c);} #define FPNEAREST LINUX_FPU_SETCW(_FPU_RC_NEAREST); #define FPPLUSINF LINUX_FPU_SETCW(_FPU_RC_UP); #define FPMINUSINF LINUX_FPU_SETCW(_FPU_RC_DOWN); #define FPTOZERO LINUX_FPU_SETCW(_FPU_RC_ZERO); #else #define _FPU_DEFAULT_p_FPU_RC_UP 0x1b72 #define _FPU_DEFAULT_p_FPU_RC_DOWN 0x1772 #define _FPU_DEFAULT_p_FPU_RC_ZERO 0x1f72 #define FPNEAREST __setfpucw(_FPU_DEFAULT); #define FPPLUSINF __setfpucw(_FPU_DEFAULT_p_FPU_RC_UP); #define FPMINUSINF __setfpucw(_FPU_DEFAULT_p_FPU_RC_DOWN); #define FPTOZERO __setfpucw(_FPU_DEFAULT_p_FPU_RC_ZERO); #endif #endif #if defined(__osf__) #if 0 #include #define FPNEAREST write_rnd(FP_RND_RN); #define FPPLUSINF write_rnd(FP_RND_RP); #define FPMINUSINF write_rnd(FP_RND_RM); #define FPTOZERO write_rnd(FP_RND_RZ); #else #define FPNEAREST #define FPPLUSINF #define FPMINUSINF #define FPTOZERO #endif #endif #if defined(__FreeBSD__) #include #define FPNEAREST fpsetround(FP_RN); #define FPPLUSINF fpsetround(FP_RP); #define FPMINUSINF fpsetround(FP_RM); #define FPTOZERO fpsetround(FP_RZ); #endif /* no control function of floating point rounding */ #ifndef FPNEAREST #define FPNEAREST fprintf(stderr, "Fpu control FPNEAREST is not supported in this machine yet.\n"); #endif #ifndef FPMINUSINF #define FPMINUSINF fprintf(stderr, "Fpu control FPMINUSINF is not supported in this machine yet.\n"); #endif #ifndef FPPLUSINF #define FPPLUSINF fprintf(stderr, "Fpu control FPPLUSINF is not supported in this machine yet.\n"); #endif #ifndef FPTOZERO #define FPTOZERO fprintf(stderr, "Fpu control FPZERO is not supported in this machine yet.\n"); #endif #define FPCLEAR FPNEAREST #define MID_PRINTF_G 2 #define MID_PRINTF_E 3 #undef N_C #undef N_M #undef N_LM #undef N_GF2N #undef N_GFPN #undef N_GFS #undef N_GFSN #define N_IP (N_B+1) #define N_ID (N_B+2) #define N_IT (N_B+3) #define N_IF (N_B+4) #define N_PRE_C N_IF #define N_C (N_B+5) #define N_M (N_B+6) #define N_LM (N_B+7) #define N_GF2N (N_B+8) #define N_GFPN (N_B+9) #define N_GFS (N_B+10) #define N_GFSN (N_B+11) /* data structures */ struct oItv { short id; char nid; char pad; Num inf; Num sup; }; typedef struct oItv *Itv; struct oItvD { short id; char nid; char pad; double inf; double sup; }; typedef struct oItvD *ItvD; struct oforth { short sign; short exp; double u; double l; }; typedef struct oforth *forth; struct oItvT { short id; char nid; char pad; forth inf; forth sup; }; typedef struct oItvT *ItvT; struct oItvF { short id; char nid; char pad; BF inf; BF sup; }; typedef struct oItvF *ItvF; extern int zerorewrite; /* general macros */ #define INF(p) ((p)->inf) #define SUP(p) ((p)->sup) #define NEWItvP(q) ((q)=(Itv)MALLOC(sizeof(struct oItv)),\ OID(q)=O_N,NID(q)=N_IP) #define NEWItvD(q) ((q)=(ItvD)MALLOC(sizeof(struct oItvD)),\ OID(q)=O_N,NID(q)=N_ID) #define NEWItvT(q) ((q)=(ItvD)MALLOC(sizeof(struct oItvT)),\ OID(q)=O_N,NID(q)=N_IT) #define NEWItvF(q) ((q)=(ItvF)MALLOC(sizeof(struct oItvF)),\ OID(q)=O_N,NID(q)=N_IF) #define MKItvP(a,b,c) (NEWItvP(c),(INF(c)=(a),SUP(c)=(b))) #define MKItvD(a,b,c) if((zerorewrite) && ((a)<=0.0) && ((b)>=0.0)) (c)=0;\ else (NEWItvD(c),(INF(c)=(a),SUP(c)=(b))) #define MKItvT(a,b,c) (NEWItvT(c),(INF(c)=(a),SUP(c)=(b))) #define MKItvF(a,b,c) (NEWItvF(c),(INF(c)=(a),SUP(c)=(b))) #define ToItvP(a,c) (NEWItvP(c),INF(c)=(a),SUP(c)=(a)) #define ToItvD(a,c) (NEWItvD(c),INF(c)=(ToReal(a)),SUP(c)=(ToReal(a))) #define ToItvF(a,c) (NEWItvF(c),INF(c)=(a),SUP(c)=(a)) #define ITVP(a) (NID(a)==N_IP) #define ITVD(a) (NID(a)==N_ID) #define ITVT(a) (NID(a)==N_IT) #define ITVF(a) (NID(a)==N_IF) /*** engine/itv.c ***/ double ToRealSup(Num); double ToRealInf(Num); double RatnToRealUp(Q); double NatToRealUp(N, int *); void double2bf(double, BF *); double bf2double(BF); void itvtois(Itv, Num *, Num *); void istoitv(Num, Num, Itv *); #if 0 void additv(Num, Num, Num *); void subitv(Num, Num, Num *); void mulitv(Num, Num, Num *); int initv(Num, double); int itvinitv(Itv, Itv); void divitv(Num, Num, Num *); void sqrtitv(Num, Num *); void pwritv(Num, Num, Num *); void pwritv0(Itv, int, Itv *); void chsgnitv(Num, Num *); void miditv(Itv, Num *); int cmpitv(Itv, Itv); void cupitv(Itv, Itv, Itv *); void capitv(Itv, Itv, Itv *); void widthitv(Itv, Num *); void distanceitv(Itv, Itv, Num *); #endif void additvp(Itv, Itv, Itv *); void subitvp(Itv, Itv, Itv *); void mulitvp(Itv, Itv, Itv *); void divitvp(Itv, Itv, Itv *); void chsgnitvp(Itv, Itv *); int cmpitvp(Itv, Itv); void pwritvp(Itv, Num, Itv *); void pwritv0p(Itv, int, Itv *); void cupitvp(Itv, Itv, Itv *); void capitvp(Itv, Itv, Itv *); void widthitvp(Itv, Num *); void distanceitvp(Itv, Itv, Num *); void ToBf(Num, BF *); void double2bf(double, BF *); double bf2double(BF); void addulp(ItvF, ItvF *); void getulp(BF, BF *); void additvf(ItvF, ItvF, ItvF *); void subitvf(ItvF, ItvF, ItvF *); void mulitvf(ItvF, ItvF, ItvF *); void divitvf(ItvF, ItvF, ItvF *); void chsgnitvf(Itv, Itv *); int cmpitvf(Itv, Itv); void pwritvf(Itv, Num, Itv *); void pwritv0f(Itv, int, Itv *); /*** engine/d-itv.c ***/ double ToRealDown(Num); double ToRealUp(Num); void Num2double(Num, double *, double *); void additvd(Num, Num, ItvD *); void subitvd(Num, Num, ItvD *); void mulitvd(Num, Num, ItvD *); void divitvd(Num, Num, ItvD *); void chsgnitvd(ItvD, ItvD *); int cmpitvd(ItvD, ItvD); void pwritvd(Num, Num, ItvD *); void pwritv0d(ItvD, int, ItvD *); void miditvd(ItvD, Num *); void cupitvd(ItvD, ItvD, ItvD *); void capitvd(ItvD, ItvD, ItvD *); void widthitvd(ItvD, Num *); void absitvd(ItvD, Num *); void distanceitvd(ItvD, ItvD, Num *); #endif /* end of INTERVAL */ #endif /* end of _INTERVAL_H */