/*
* $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.3 2001/03/19 01:13:50 noro Exp $
*/
#ifndef _INTERVAL_H
#define _INTERVAL_H
#define PRINTF_G 0
#define PRINTF_E 1
#if defined(INTERVAL)
#include <math.h>
#ifdef sun
#if OSMajorVersion < 4 || (OSMajorVersion == 4 && OSMinorVersion < 1)
#include <sys/ieeefp.h>
#endif
#if defined(__svr4__)
#include <ieeefp.h>
#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 <fpu_control.h>
#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 <float.h>
#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 <floatingpoint.h>
#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
#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)
/* 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 */