===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/bfaux.c,v
retrieving revision 1.3
retrieving revision 1.6
diff -u -p -r1.3 -r1.6
--- OpenXM_contrib2/asir2000/builtin/bfaux.c	2015/08/05 01:23:02	1.3
+++ OpenXM_contrib2/asir2000/builtin/bfaux.c	2015/08/07 06:15:00	1.6
@@ -1,17 +1,23 @@
-/* $OpenXM: OpenXM_contrib2/asir2000/builtin/bfaux.c,v 1.2 2015/08/05 01:10:38 noro Exp $ */
+/* $OpenXM: OpenXM_contrib2/asir2000/builtin/bfaux.c,v 1.5 2015/08/07 05:30:35 noro Exp $ */
 #include "ca.h"
 #include "parse.h"
 
-void Peval(), Psetprec(), Psetbprec(), Ptodouble();
+void Peval(), Psetprec(), Psetbprec(), Ptodouble(), Psetround();
+void Pmpfr_gamma(), Pmpfr_floor();
 
 struct ftab bf_tab[] = {
 	{"eval",Peval,-2},
 	{"setprec",Psetprec,-1},
 	{"setbprec",Psetbprec,-1},
+	{"setround",Psetround,-1},
 	{"todouble",Ptodouble,1},
+	{"mpfr_gamma",Pmpfr_gamma,-2},
+	{"mpfr_floor",Pmpfr_floor,-1},
 	{0,0,0},
 };
 
+int mpfr_roundmode = MPFR_RNDN;
+
 void Ptodouble(NODE arg,Num *rp)
 {
 	double r,i;
@@ -43,9 +49,7 @@ void Ptodouble(NODE arg,Num *rp)
 	}
 }
 
-void Peval(arg,rp)
-NODE arg;
-Obj *rp;
+void Peval(NODE arg,Obj *rp)
 {
   int prec;
 
@@ -103,6 +107,43 @@ void Psetbprec(NODE arg,Obj *rp)
 	mpfr_set_default_prec(prec);
 }
 
+void Psetround(NODE arg,Q *rp)
+{
+  int round;
+
+  STOQ(mpfr_roundmode,*rp);
+  if ( arg ) {
+		asir_assert(ARG0(arg),O_N,"setround");
+    round = QTOS((Q)ARG0(arg)); 
+    switch ( round ) {
+    case 0:
+      mpfr_roundmode = MPFR_RNDN;
+      break;
+    case 1:
+      mpfr_roundmode = MPFR_RNDZ;
+      break;
+    case 2:
+      mpfr_roundmode = MPFR_RNDU;
+      break;
+    case 3:
+      mpfr_roundmode = MPFR_RNDD;
+      break;
+    case 4:
+      mpfr_roundmode = MPFR_RNDA;
+      break;
+    case 5:
+      mpfr_roundmode = MPFR_RNDF;
+      break;
+    case 6:
+      mpfr_roundmode = MPFR_RNDNA;
+      break;
+    default:
+      error("setround : invalid rounding mode");
+      break;
+    }
+  }
+}
+
 Num tobf(Num a,int prec);
 
 void mp_pi(NODE arg,BF *rp)
@@ -113,7 +154,7 @@ void mp_pi(NODE arg,BF *rp)
 	prec = arg ? QTOS((Q)ARG0(arg)) : 0;
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_const_pi(r->body,MPFR_RNDN);
+	mpfr_const_pi(r->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -127,8 +168,8 @@ void mp_e(NODE arg,BF *rp)
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
 	mpfr_init(one);
-	mpfr_set_ui(one,1,MPFR_RNDN);
-	mpfr_exp(r->body,one,MPFR_RNDN);
+	mpfr_set_ui(one,1,mpfr_roundmode);
+	mpfr_exp(r->body,one,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -142,7 +183,7 @@ void mp_sin(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_sin(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_sin(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -156,7 +197,7 @@ void mp_cos(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_cos(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_cos(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -170,7 +211,7 @@ void mp_tan(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_tan(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_tan(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -184,7 +225,7 @@ void mp_asin(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_asin(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_asin(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 void mp_acos(NODE arg,BF *rp)
@@ -197,7 +238,7 @@ void mp_acos(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_acos(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_acos(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 void mp_atan(NODE arg,BF *rp)
@@ -210,7 +251,7 @@ void mp_atan(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_atan(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_atan(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -224,7 +265,7 @@ void mp_sinh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_sinh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_sinh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -238,7 +279,7 @@ void mp_cosh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_cosh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_cosh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -252,7 +293,7 @@ void mp_tanh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_tanh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_tanh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -266,7 +307,7 @@ void mp_asinh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_asinh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_asinh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 void mp_acosh(NODE arg,BF *rp)
@@ -279,7 +320,7 @@ void mp_acosh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_acosh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_acosh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 void mp_atanh(NODE arg,BF *rp)
@@ -292,7 +333,7 @@ void mp_atanh(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_atanh(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_atanh(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -306,7 +347,7 @@ void mp_exp(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_exp(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_exp(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -320,7 +361,7 @@ void mp_log(NODE arg,BF *rp)
 	a = tobf(ARG0(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_log(r->body,((BF)a)->body,MPFR_RNDN);
+	mpfr_log(r->body,((BF)a)->body,mpfr_roundmode);
     *rp = r; 
 }
 
@@ -335,6 +376,36 @@ void mp_pow(NODE arg,BF *rp)
 	e = tobf(ARG1(arg),prec);
 	NEWBF(r);
 	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
-	mpfr_pow(r->body,((BF)a)->body,((BF)e)->body,MPFR_RNDN);
+	mpfr_pow(r->body,((BF)a)->body,((BF)e)->body,mpfr_roundmode);
     *rp = r; 
+}
+
+void Pmpfr_gamma(NODE arg,BF *rp)
+{
+	Num a;
+  int prec;
+	BF r;
+
+	prec = NEXT(arg) ? QTOS((Q)ARG1(arg)) : 0;
+  prec *= 3.32193;
+	a = tobf(ARG0(arg),prec);
+	NEWBF(r);
+	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
+	mpfr_gamma(r->body,((BF)a)->body,mpfr_roundmode);
+  *rp = r; 
+}
+
+void Pmpfr_floor(NODE arg,BF *rp)
+{
+	Num a;
+  int prec;
+	BF r;
+
+	prec = NEXT(arg) ? QTOS((Q)ARG1(arg)) : 0;
+  prec *= 3.32193;
+	a = tobf(ARG0(arg),prec);
+	NEWBF(r);
+	prec ? mpfr_init2(r->body,prec) : mpfr_init(r->body);
+	mpfr_floor(r->body,((BF)a)->body);
+  *rp = r; 
 }