[BACK]Return to interval.h CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / include

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

version 1.16, 2019/06/04 07:11:22 version 1.17, 2019/11/12 10:52:04
Line 1 
Line 1 
 /*  /*
  * $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.15 2018/03/29 01:32:53 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/include/interval.h,v 1.16 2019/06/04 07:11:22 kondoh Exp $
 */  */
 #ifndef  _INTERVAL_H  #ifndef  _INTERVAL_H
 #define  _INTERVAL_H  #define  _INTERVAL_H
   
 #define       INT_ASIR_VERSION        20181218  #define       INT_ASIR_VERSION        20191111
   
 #define  PRINTF_G  0  #define  PRINTF_G  0
 #define  PRINTF_E  1  #define  PRINTF_E  1
   
 #if defined(INTERVAL)  #if defined(INTERVAL)
   #if INTERVAL == 0
   #undef INTERVAL
   #endif
   #endif
   #if defined(INTERVAL)
   
 #include <math.h>  #include <math.h>
   
Line 129  static char  *Interval_dummy;
Line 134  static char  *Interval_dummy;
 #define  MID_PRINTF_G  2  #define  MID_PRINTF_G  2
 #define  MID_PRINTF_E  3  #define  MID_PRINTF_E  3
   
 #define DEFAULTPREC     0       /* for mpfr */  #define DEFAULTPREC     0 //mpfr_get_default_prec()     /* for mpfr */
   
 /* data structures */  /* data structures */
 struct oItv {  struct oItv {
Line 184  typedef struct oIntervalBigFloat *IntervalBigFloat;
Line 189  typedef struct oIntervalBigFloat *IntervalBigFloat;
 extern int zerorewrite;  extern int zerorewrite;
 extern int zerorewriteCount;  extern int zerorewriteCount;
   
   #define ZEROREWRITE       if (zerorewrite && initvp(0,c) ) { *rp = 0; zerorewriteCount++; }
   
 /* general macros */  /* general macros */
 #define INF(p)  ((p)->inf)  #define INF(p)  ((p)->inf)
 #define SUP(p)  ((p)->sup)  #define SUP(p)  ((p)->sup)
Line 203  extern int zerorewriteCount;
Line 210  extern int zerorewriteCount;
 #define MKIntervalBigFloat(a,b,c)  (NEWIntervalBigFloat(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 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 ToIntervalDouble(a,c)  (NEWIntervalDouble(c),INF(c)=(toRealDown(a)),SUP(c)=(toRealUp(a)))
 #define ToIntervalBigFloat(a,c)  (NEWIntervalBigFloat(c),INF(c)=(a),SUP(c)=(a))  #define ToIntervalBigFloat(a,c)  (NEWIntervalBigFloat(c),INF(c)=(a),SUP(c)=(a))
   
 #define ITVP(a) (NID(a)==N_IP)  #define ITVP(a) (NID(a)==N_IP)
Line 211  extern int zerorewriteCount;
Line 218  extern int zerorewriteCount;
 #define ITVQ(a) (NID(a)==N_IntervalQuad)  #define ITVQ(a) (NID(a)==N_IntervalQuad)
 #define ITVF(a) (NID(a)==N_IntervalBigFloat)  #define ITVF(a) (NID(a)==N_IntervalBigFloat)
   
   #define EvalIntervalDouble              0
   #define EvalIntervalQuad                1
   #define EvalIntervalBigFloat    2
   
 #if 0  #if 0
 double  ToRealSup(Num);  double  ToRealSup(Num);
 double  ToRealInf(Num);  double  ToRealInf(Num);
Line 257  void miditvp(Itv, Num *);
Line 268  void miditvp(Itv, Num *);
 void    cupitvp(Itv, Itv, Itv *);  void    cupitvp(Itv, Itv, Itv *);
 void    capitvp(Itv, Itv, Itv *);  void    capitvp(Itv, Itv, Itv *);
 void    widthitvp(Itv, Num *);  void    widthitvp(Itv, Num *);
 void absitvp(Itv, Num *);  void    absitvp(Itv, Num *);
   void    absintvalp(Itv, Itv *);
 void    distanceitvp(Itv, Itv, Num *);  void    distanceitvp(Itv, Itv, Num *);
   
   
Line 268  void    distanceitvp(Itv, Itv, Num *);
Line 280  void    distanceitvp(Itv, Itv, Num *);
 //double  bf2double(BF);  //double  bf2double(BF);
 //void  addulp(IntervalBigFloat, IntervalBigFloat *);  //void  addulp(IntervalBigFloat, IntervalBigFloat *);
 //void  getulp(BF, BF *);  //void  getulp(BF, BF *);
   double mpfr2dblDown(mpfr_t);
   double mpfr2dblUp(mpfr_t);
   void toInterval(Num, int, int, Num *);
   
 void    additvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);  void    additvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);
 void    subitvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);  void    subitvf(IntervalBigFloat, IntervalBigFloat, IntervalBigFloat *);
Line 287  void absitvf(Itv, Num *);
Line 302  void absitvf(Itv, Num *);
 void distanceitvf(Itv, Itv, Num *);  void distanceitvf(Itv, Itv, Num *);
   
 /***    engine/d-itv.c    ***/  /***    engine/d-itv.c    ***/
 double  ToRealDown(Num);  double  toRealDown(Num);
 double  ToRealUp(Num);  double  toRealUp(Num);
 void  Num2double(Num, double *, double *);  void  Num2double(Num, double *, double *);
   
 void    additvd(Num, Num, IntervalDouble *);  void    additvd(Num, Num, IntervalDouble *);
Line 305  void    cupitvd(IntervalDouble, IntervalDouble, Interv
Line 320  void    cupitvd(IntervalDouble, IntervalDouble, Interv
 void    capitvd(IntervalDouble, IntervalDouble, IntervalDouble *);  void    capitvd(IntervalDouble, IntervalDouble, IntervalDouble *);
 void    widthitvd(IntervalDouble, Num *);  void    widthitvd(IntervalDouble, Num *);
 void    absitvd(IntervalDouble, Num *);  void    absitvd(IntervalDouble, Num *);
   void    absintvald(IntervalDouble, IntervalDouble *);
 void    distanceitvd(IntervalDouble, IntervalDouble, Num *);  void    distanceitvd(IntervalDouble, IntervalDouble, Num *);
   
   /***    builtin/itvnum.c    ***/
   
   void evalitvr(VL ,Obj ,int , int , Obj *);
   void evalitvp(VL ,P ,int , int , P *);
   void evalitvv(VL ,V ,int , int , Obj *);
   //void evalitvins(PFINS ,int , int , Obj *);
   extern void (*pi_itv_ft[])();
   extern void (*e_itv_ft[])();
   extern void (*sin_itv_ft[])();
   extern void (*cos_itv_ft[])();
   extern void (*tan_itv_ft[])();
   extern void (*asin_itv_ft[])();
   extern void (*acos_itv_ft[])();
   extern void (*atan_itv_ft[])();
   extern void (*sinh_itv_ft[])();
   extern void (*cosh_itv_ft[])();
   extern void (*tanh_itv_ft[])();
   extern void (*asinh_itv_ft[])();
   extern void (*acosh_itv_ft[])();
   extern void (*atanh_itv_ft[])();
   extern void (*exp_itv_ft[])();
   extern void (*log_itv_ft[])();
   extern void (*abs_itv_ft[])();
   extern void (*pow_itv_ft[])();
   //void devalr(VL ,Obj ,Obj *);
   //void devalp(VL ,P ,P *);
   //void devalv(VL ,V ,Obj *);
   //void devalins(PFINS ,Obj *);
   
   void Pitvbf_pi(NODE ,Obj *);
   void Pitvbf_e(NODE ,Obj *);
   void Pitvbf_sin(NODE ,Obj *);
   void Pitvbf_cos(NODE ,Obj *);
   void Pitvbf_tan(NODE ,Obj *);
   void Pitvbf_asin(NODE ,Obj *);
   void Pitvbf_acos(NODE ,Obj *);
   void Pitvbf_atan(NODE ,Obj *);
   void Pitvbf_sinh(NODE ,Obj *);
   void Pitvbf_cosh(NODE ,Obj *);
   void Pitvbf_tanh(NODE ,Obj *);
   void Pitvbf_asinh(NODE ,Obj *);
   void Pitvbf_acosh(NODE ,Obj *);
   void Pitvbf_atanh(NODE ,Obj *);
   void Pitvbf_exp(NODE ,Obj *);
   void Pitvbf_log(NODE ,Obj *);
   void Pitvbf_abs(NODE ,Obj *);
   //void mp_factorial(NODE ,Num *);
   void Pitvbf_pow(NODE ,Num *);
   
   void Pitvd_pi(NODE ,Obj *);
   void Pitvd_e(NODE ,Obj *);
   void Pitvd_sin(NODE ,Obj *);
   void Pitvd_cos(NODE ,Obj *);
   void Pitvd_tan(NODE ,Obj *);
   void Pitvd_asin(NODE ,Obj *);
   void Pitvd_acos(NODE ,Obj *);
   void Pitvd_atan(NODE ,Obj *);
   void Pitvd_sinh(NODE ,Obj *);
   void Pitvd_cosh(NODE ,Obj *);
   void Pitvd_tanh(NODE ,Obj *);
   void Pitvd_asinh(NODE ,Obj *);
   void Pitvd_acosh(NODE ,Obj *);
   void Pitvd_atanh(NODE ,Obj *);
   void Pitvd_exp(NODE ,Obj *);
   void Pitvd_log(NODE ,Obj *);
   void Pitvd_abs(NODE ,Obj *);
   void Pitvd_pow(NODE ,Num *);
   
 #endif /* end of INTERVAL */  #endif /* end of INTERVAL */
 #endif /* end of _INTERVAL_H */  #endif /* end of _INTERVAL_H */

Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17

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