Support MPFR floats in numkondov1.4

Remove '-D' option (error tolerance) in numkondo

1 files changed, 14 insertions, 11 deletions

diff --git a/src/rational_float.c b/src/rational_float.c
index eebc4f4..d6a6a4c 100644
--- a/src/rational_float.c
+++ b/src/rational_float.c
@@ -19,21 +19,14 @@ limitations under the License.
 #include "rational_float.h"
 #include <stdio.h>
 #include <stdlib.h>
+#include <stdarg.h>
+// define MPFR_USE_VA_LIST to enable the use of mpfr_inits and mpfr_clears
+#define MPFR_USE_VA_LIST
+#include <mpfr.h>
 #include "istring.h"
 #include "array.h"
 #include "math.h"
 
-Q quot(long double p, long double q){
-  Q ret;
-  if(q==0){
-    fprintf(stderr,"error: %Lf/%Lf is ill defined\n",p,q);
-    exit(-1);
-  }
-  ret.numerator=p;
-  ret.denominator=q;
-  return(ret);
-}
-
 // add
 Q Q_add(Q x1,Q x2){
   Q ret;
@@ -141,6 +134,16 @@ long double lcm(long double x,long double y){
 double Q_double_value(Q q){
   return(1.0*q.numerator/q.denominator);
 }
+// approximate value as mpfr float
+int Q_mpfr_value(mpfr_t out, Q q){
+  mpfr_t x;
+  mpfr_init(out);
+  mpfr_init(x);
+  mpfr_set_ld(x, q.denominator, MPFR_RNDN);
+  mpfr_ld_div(out, q.numerator, x, MPFR_RNDN);
+  mpfr_clear(x);
+  return(0);
+}
 
 
 // print to string