OpenXM_contrib2/asir2000/include/interval.h - diff

Return to interval.h CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / include

Diff for /OpenXM_contrib2/asir2000/include/interval.h between version 1.12 and 1.16

version 1.12, 2015/08/06 10:01:52

version 1.16, 2019/06/04 07:11:22

Line 1

* $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.11 2014/01/07 06:22:08 saito Exp $

* $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.15 2018/03/29 01:32:53 noro Exp $

#ifndef _INTERVAL_H

#define _INTERVAL_H

#define PRINTF_G 0

#define INT_ASIR_VERSION 20181218

#define PRINTF_E 1

#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

#ifdef linux

#if 0

#include <fpu_control.h>

#ifdef linux

#include <fpu_control.h>

#if 1

#define LINUX_FPU_RC_MASK 0xf3ff

#define LINUX_FPU_SETCW(c) \

#if defined(__ARM_ARCH) || defined(ANDROID)

#ifndef _FPU_RC_NEAREST

#define _FPU_RC_NEAREST 0x0

#endif

#ifndef _FPU_RC_DOWN

#define _FPU_RC_DOWN 0x400

#endif

#ifndef _FPU_RC_UP

#define _FPU_RC_UP 0x800

#endif

#ifndef _FPU_RC_ZERO

#define _FPU_RC_ZERO 0xC00

#endif

#endif /* __ARM_ARCH */

#define LINUX_FPU_RC_MASK 0xf3ff

#define LINUX_FPU_SETCW(c) \

{fpu_control_t __tmp__; _FPU_GETCW(__fpu_control);\

__tmp__ = (__fpu_control & LINUX_FPU_RC_MASK | c);\

_FPU_SETCW(__tmp__);}

#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

#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

#if defined(__FreeBSD__) && defined(__GNUC__)

#include <floatingpoint.h>

#define FPNEAREST fpsetround(FP_RN);

#define FPPLUSINF fpsetround(FP_RP);

#define FPMINUSINF fpsetround(FP_RM);

#define FPTOZERO fpsetround(FP_RZ);

#endif

#if defined(VISUAL) || defined(__MINGW32__) || defined(__MINGW64__)

#if defined(VISUAL) || defined(__MINGW32__)

#include <float.h>

#define FPNEAREST _controlfp(_RC_NEAR,_MCW_RC);

#define FPPLUSINF _controlfp(_RC_UP,_MCW_RC);

#define FPMINUSINF _controlfp(_RC_DOWN,_MCW_RC);

#define FPTOZERO _controlfp(_RC_CHOP,_MCW_RC);

#endif

#include <fenv.h>

#define FPNEAREST fesetround(FE_TONEAREST);

#define FPPLUSINF fesetround(FE_UPWARD);

#define FPMINUSINF fesetround(FE_DOWNWARD);

#define FPTOZERO fesetround(FE_TOWARDZERO);

/* 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

#define DEFAULTPREC 0 /* for mpfr */

/* data structures */

struct oItv {

short id;

char nid;

char pad;

Num inf;

Num sup;

};

typedef struct oItv *Itv;

struct oIntervalDouble {

short id;

char nid;

char pad;

double inf;

double sup;

};

typedef struct oIntervalDouble *IntervalDouble;

struct oforth {

short sign;

short exp;

double u;

double l;

};

typedef struct oforth *forth;

struct oIntervalQuad {

short id;

char nid;

char pad;

forth inf;

forth sup;

};

typedef struct oIntervalQuad *IntervalQuad;

struct oIntervalBigFloat {

short id;

char nid;

char pad;

BF inf;

BF sup;

};

typedef struct oIntervalBigFloat *IntervalBigFloat;

extern int zerorewrite;

extern int zerorewriteCount;

/* 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 NEWIntervalDouble(q) ((q)=(IntervalDouble)MALLOC(sizeof(struct oIntervalDouble)),\

OID(q)=O_N,NID(q)=N_IntervalDouble)

#define NEWIntervalQuad(q) ((q)=(IntervalDouble)MALLOC(sizeof(struct oIntervalQuad)),\

OID(q)=O_N,NID(q)=N_IntervalQuad)

#define NEWIntervalBigFloat(q) ((q)=(IntervalBigFloat)MALLOC(sizeof(struct oIntervalBigFloat)),\

OID(q)=O_N,NID(q)=N_IntervalBigFloat)

#define MKItvP(a,b,c) (NEWItvP(c),(INF(c)=(a),SUP(c)=(b)))

#define MKIntervalDouble(a,b,c) if((zerorewrite) && ((a)<=0.0) && ((b)>=0.0)) (c)=0;\

#define MKIntervalDouble(a,b,c) if((zerorewrite) && ((a)<=0.0) && ((b)>=0.0)) {(c)=0;zerorewriteCount++;}\

else (NEWIntervalDouble(c),(INF(c)=(a),SUP(c)=(b)))

#define MKIntervalQuad(a,b,c) (NEWIntervalQuad(c),(INF(c)=(a),SUP(c)=(b)))

#define MKIntervalBigFloat(a,b,c) (NEWIntervalBigFloat(c),(INF(c)=(a),SUP(c)=(b)))

#define ToItvP(a,c) (NEWItvP(c),INF(c)=(a),SUP(c)=(a))

#define ToIntervalDouble(a,c) (NEWIntervalDouble(c),INF(c)=(ToReal(a)),SUP(c)=(ToReal(a)))

#define ToIntervalBigFloat(a,c) (NEWIntervalBigFloat(c),INF(c)=(a),SUP(c)=(a))

#define ITVP(a) (NID(a)==N_IP)

#define ITVD(a) (NID(a)==N_IntervalDouble)

#define ITVQ(a) (NID(a)==N_IntervalQuad)

#define ITVF(a) (NID(a)==N_IntervalBigFloat)

/*** engine/itv.c ***/

#if 0

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 *);

Line 212 void widthitv(Itv, Num *);

Line 239 void widthitv(Itv, Num *);

void distanceitv(Itv, Itv, Num *);

#endif

/*** engine/p-itv.c ***/

void itvtois(Itv, Num *, Num *);

void istoitv(Num, Num, Itv *);

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 chsgnitvp(Itv, Itv *);

int initvp(Num , Itv );

int itvinitvp(Itv, Itv);

int cmpitvp(Itv, Itv);

void miditvp(Itv, Num *);

void cupitvp(Itv, Itv, Itv *);

void capitvp(Itv, Itv, Itv *);

void widthitvp(Itv, Num *);

void absitvp(Itv, Num *);

void distanceitvp(Itv, Itv, Num *);

void ToBf(Num, BF *);

void double2bf(double, BF *);

double bf2double(BF);

void addulp(IntervalBigFloat, IntervalBigFloat *);

void getulp(BF, BF *);

/*** engine/f-itv.c ***/

//void ToBf(Num, BF *);

//void double2bf(double, BF *);

//double bf2double(BF);

//void addulp(IntervalBigFloat, IntervalBigFloat *);

//void getulp(BF, BF *);

void additvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);

void subitvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);

void mulitvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);

void divitvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);

void chsgnitvf(Itv, Itv *);

void chsgnitvf(IntervalBigFloat, IntervalBigFloat *);

int initvf(Num, Itv);

int itvinitvf(Itv, Itv);

int cmpitvf(Itv, Itv);

void pwritvf(Itv, Num, Itv *);

void pwritv0f(Itv, int, Itv *);

void miditvf(Itv, Num *);

void cupitvf(Itv, Itv, Itv *);

void capitvf(Itv, Itv, Itv *);

void widthitvf(Itv, Num *);

void absitvf(Itv, Num *);

void distanceitvf(Itv, Itv, Num *);

/*** engine/d-itv.c ***/

double ToRealDown(Num);

double ToRealUp(Num);

void Num2double(Num, double *, double *);

void additvd(Num, Num, IntervalDouble *);

void subitvd(Num, Num, IntervalDouble *);

void mulitvd(Num, Num, IntervalDouble *);

void divitvd(Num, Num, IntervalDouble *);

void chsgnitvd(IntervalDouble, IntervalDouble *);

int initvd(Num, IntervalDouble);

int cmpitvd(IntervalDouble, IntervalDouble);

void pwritvd(Num, Num, IntervalDouble *);

void pwritv0d(IntervalDouble, int, IntervalDouble *);

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>