Support MPFR floats in numkondov1.4

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

1 files changed, 12 insertions, 20 deletions

diff --git a/src/types.h b/src/types.h
index a5d47b9..d30840e 100644
--- a/src/types.h
+++ b/src/types.h
@@ -21,6 +21,8 @@ limitations under the License.
 #ifndef TYPES_H
 #define TYPES_H
 
+#include <mpfr.h>
+
 
 // rational number
 typedef struct Q{
@@ -114,6 +116,12 @@ typedef struct RCC{
   int* indices;
   int length;
 } RCC;
+// rcc using mpfr floats
+typedef struct RCC_mpfr{
+  mpfr_t* values;
+  int* indices;
+  int length;
+} RCC_mpfr;
 
 // identities between fields
 typedef struct Identities{
@@ -174,25 +182,6 @@ typedef struct Id_Table{
   int memory;
 } Id_Table;
 
-/*
-// polynomial scalar and vectors
-typedef struct Polynomial_Scalar{
-  Coefficient coef;
-  int* indices;
-  int length;
-} Polynomial_Scalar;
-typedef struct Polynomial_Vector{
-  Coefficient* coefv;
-  int* indices;
-  int length;
-} Polynomial_Vector;
-typedef struct Polynomial_Matrix{
-  Coefficient** coefm;
-  int* indices;
-  int length;
-} Polynomial_Matrix;
-*/
-
 
 // command line options
 typedef struct Meankondo_Options{
@@ -203,8 +192,9 @@ typedef struct Meankondo_Options{
 typedef struct Numkondo_Options{
   int display_mode;
   int niter;
-  long double tol;
   Char_Array eval_rccstring;
+  mpfr_prec_t mpfr_prec;
+  mpfr_exp_t mpfr_emax;
 } Numkondo_Options;
 
 typedef struct Meantools_Options{
@@ -213,6 +203,8 @@ typedef struct Meantools_Options{
   Int_Array deriv_vars;
   Char_Array eval_rccstring;
   int chain;
+  mpfr_prec_t mpfr_prec;
+  mpfr_exp_t mpfr_emax;
 } Meantools_Options;
 
 typedef struct Kondopp_Options{