From b0117b1d5ae73916c6f0d289be1f693bb5f46824 Mon Sep 17 00:00:00 2001
From: mbkma <johannes.unruh@fau.de>
Date: Thu, 5 Mar 2020 13:06:45 +0100
Subject: Port to GNU MPFR/MPC Library

For further information please visit:
https://www.mpfr.org/
http://www.multiprecision.org/mpc
---
 src/mp.h | 63 ++++++++++++++++++++++++++++++++++-----------------------------
 1 file changed, 34 insertions(+), 29 deletions(-)

(limited to 'src/mp.h')

diff --git a/src/mp.h b/src/mp.h
index 53943b8..1748c03 100644
--- a/src/mp.h
+++ b/src/mp.h
@@ -33,27 +33,21 @@
 #include <stdbool.h>
 #include <stdint.h>
 #include <glib.h>
+#include <glib/gi18n.h>
+#include <mpfr.h>
+#include <mpc.h>
 
-/* Size of the multiple precision values */
-#define MP_SIZE 1000
+/* If we're not using GNU C, elide __attribute__ */
+#ifndef __GNUC__
+#  define  __attribute__(x)  /*NOTHING*/
+#endif
 
-/* Base for numbers */
-#define MP_BASE 10000
+/* Precision of mpfr_t and mpc_t type objects */
+#define PRECISION 1000
 
-/* Object for a high precision floating point number representation
- *
- * x = sign * (MP_BASE^(exponent-1) + MP_BASE^(exponent-2) + ...)
- */
 typedef struct
 {
-   /* Sign (+1, -1) or 0 for the value zero */
-   int sign, im_sign;
-
-   /* Exponent (to base MP_BASE) */
-   int exponent, im_exponent;
-
-   /* Normalized fraction */
-   int fraction[MP_SIZE], im_fraction[MP_SIZE];
+    mpc_t num;
 } MPNumber;
 
 typedef enum
@@ -65,17 +59,28 @@ typedef enum
 
 /* Returns error string or NULL if no error */
 // FIXME: Global variable
-const char *mp_get_error(void);
+const char  *mp_get_error(void);
 
 /* Clear any current error */
-void mp_clear_error(void);
+void        mp_clear_error(void);
+
+void        mperr(const char *format, ...) __attribute__((format(printf, 1, 2)));
+
+/* Returns initialized MPNumber object */
+MPNumber    mp_new(void);
+
+MPNumber*   mp_new_ptr(void);
+
+void        mp_clear(MPNumber *z);
+
+void        mp_free(MPNumber *z);
 
 /* Returns:
  *  0 if x == y
  * <0 if x < y
  * >0 if x > y
  */
-int    mp_compare_mp_to_mp(const MPNumber *x, const MPNumber *y);
+int    mp_compare(const MPNumber *x, const MPNumber *y);
 
 /* Return true if the value is x == 0 */
 bool   mp_is_zero(const MPNumber *x);
@@ -134,9 +139,6 @@ void   mp_add(const MPNumber *x, const MPNumber *y, MPNumber *z);
 /* Sets z = x + y */
 void   mp_add_integer(const MPNumber *x, int64_t y, MPNumber *z);
 
-/* Sets z = x + numerator ÷ denominator */
-void   mp_add_fraction(const MPNumber *x, int64_t numerator, int64_t denominator, MPNumber *z);
-
 /* Sets z = x − y */
 void   mp_subtract(const MPNumber *x, const MPNumber *y, MPNumber *z);
 
@@ -146,9 +148,6 @@ void   mp_multiply(const MPNumber *x, const MPNumber *y, MPNumber *z);
 /* Sets z = x × y */
 void   mp_multiply_integer(const MPNumber *x, int64_t y, MPNumber *z);
 
-/* Sets z = x × numerator ÷ denominator */
-void   mp_multiply_fraction(const MPNumber *x, int64_t numerator, int64_t denominator, MPNumber *z);
-
 /* Sets z = x ÷ y */
 void   mp_divide(const MPNumber *x, const MPNumber *y, MPNumber *z);
 
@@ -217,6 +216,8 @@ void   mp_epowy(const MPNumber *x, MPNumber *z);
 /* Returns a list of all prime factors in x as MPNumbers */
 GList* mp_factorize(const MPNumber *x);
 
+GList* mp_factorize_unit64 (uint64_t n);
+
 /* Sets z = x */
 void   mp_set_from_mp(const MPNumber *x, MPNumber *z);
 
@@ -250,16 +251,16 @@ void   mp_set_from_random(MPNumber *z);
 bool   mp_set_from_string(const char *text, int default_base, MPNumber *z);
 
 /* Returns x as a native single-precision floating point number */
-float  mp_cast_to_float(const MPNumber *x);
+float  mp_to_float(const MPNumber *x);
 
 /* Returns x as a native double-precision floating point number */
-double mp_cast_to_double(const MPNumber *x);
+double mp_to_double(const MPNumber *x);
 
 /* Returns x as a native integer */
-int64_t mp_cast_to_int(const MPNumber *x);
+int64_t mp_to_integer(const MPNumber *x);
 
 /* Returns x as a native unsigned integer */
-uint64_t mp_cast_to_unsigned_int(const MPNumber *x);
+uint64_t mp_to_unsigned_integer(const MPNumber *x);
 
 /* Sets z = sin x */
 void   mp_sin(const MPNumber *x, MPAngleUnit unit, MPNumber *z);
@@ -324,4 +325,8 @@ void   mp_ones_complement(const MPNumber *x, int wordlen, MPNumber *z);
 /* Sets z to be the twos complement of x for word of length 'wordlen' */
 void   mp_twos_complement(const MPNumber *x, int wordlen, MPNumber *z);
 
+void convert_to_radians(const MPNumber *x, MPAngleUnit unit, MPNumber *z);
+
+void convert_from_radians(const MPNumber *x, MPAngleUnit unit, MPNumber *z);
+
 #endif /* MP_H */
-- 
cgit v1.2.1