1 files changed, 54 insertions, 0 deletions
diff --git a/print_raw.c b/print_raw.c
new file mode 100644
index 000000000..afc2590d6
--- /dev/null
+++ b/print_raw.c
@@ -0,0 +1,54 @@
+#include <stdio.h>
+#include "gmp.h"
+#include "gmp-impl.h"
+#include "mpfr.h"
+
+void
+mpfr_get_str_raw(char *digit_ptr, mpfr_srcptr x)
+{
+  mp_limb_t *mx, wd, t; long ex, sx, k, l, p;
+
+  mx = MANT(x); 
+  ex = EXP(x); 
+  p = PREC(x); 
+
+  if (SIGN(x) < 0) { *digit_ptr = '-'; digit_ptr++; }
+  sprintf(digit_ptr, "0."); digit_ptr += 2; 
+
+  sx = 1+(p-1)/mp_bits_per_limb; /* number of significant limbs */
+  for (k = sx - 1; k >= 0 ; k--)
+    { 
+      wd = mx[k]; 
+      t = 1UL << (BITS_PER_MP_LIMB - 1); 
+      for (l = BITS_PER_MP_LIMB - 1; l>=0; l--)
+	{
+	  if (wd & t) 
+	    { *digit_ptr = '1'; digit_ptr++; } 
+	  else 
+	    { *digit_ptr = '0'; digit_ptr++; }
+	  t >>= 1; 
+	  if (--p==0) { *digit_ptr = '['; digit_ptr++; }
+	}
+    }
+  sprintf(digit_ptr, "]E%ld", ex); 
+}
+ 
+void
+mpfr_print_raw(mpfr_srcptr x)
+{
+  char *str; 
+
+  /* 3 char for sign + 0 + binary point
+     + ABSSIZE(x) * BITS_PER_MP_LIMB for mantissa
+     + 2 for brackets in mantissa
+     + 1 for 'E'
+     + 11 for exponent (including sign)
+     = 17 + ABSSIZE(x) * BITS_PER_MP_LIMB
+   */
+  str = (char *) malloc((17 + ABSSIZE(x) * BITS_PER_MP_LIMB)*sizeof(char)); 
+  mpfr_get_str_raw(str, x);
+
+  printf("%s", str); 
+  free(str); 
+}
+