From 77ac9ce0a5c55d4f79f8fb8f7daa59ddb53cb507 Mon Sep 17 00:00:00 2001
From: Florian Fischer <florian.fl.fischer@fau.de>
Date: Sat, 24 Aug 2019 17:57:51 +0200
Subject: add cfrac benchmark

---
 src/benchmarks/cfrac/pfloat.c | 61 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 61 insertions(+)
 create mode 100644 src/benchmarks/cfrac/pfloat.c

(limited to 'src/benchmarks/cfrac/pfloat.c')

diff --git a/src/benchmarks/cfrac/pfloat.c b/src/benchmarks/cfrac/pfloat.c
new file mode 100644
index 0000000..63f4344
--- /dev/null
+++ b/src/benchmarks/cfrac/pfloat.c
@@ -0,0 +1,61 @@
+/*
+ *	High Precision Math Library Supplement for floating point routines
+ */
+#include <stdio.h>
+#include <math.h>
+#include "pdefs.h"
+#include "pcvt.h"
+#include "precision.h"
+
+extern precision palloc();
+
+/*
+ * double to precision
+ */
+precision dtop(f)
+   register double f;
+{
+   register digitPtr	 uPtr;
+   register precision    u;
+
+   u = palloc(DOUBLESIZE);		/* pretty big */
+   if (u == pUndef) return u;
+
+   if (f < 0.0) {
+      f = -f;
+      u->sign = true;
+   } else {
+      u->sign = false;
+   }
+   uPtr	      = u->value;
+   do {
+      *uPtr++ = fmod(f, (double) BASE);
+      f	      = floor(f / (double) BASE);
+   } while (f != 0.0);
+
+   u->size = (uPtr - u->value);
+
+   return presult(u);
+}
+
+/*
+ *  precision to double (no overflow check)
+ */
+double ptod(u)
+   precision	  u;
+{
+   register digitPtr uPtr;
+   register double   f;
+
+   (void) pparm(u);
+   uPtr = u->value + u->size;
+   f    = 0.0;
+   do {
+      f = f * (double) BASE + (double) *--uPtr;
+   } while (uPtr > u->value);
+
+   if (u->sign) f = -f;
+
+   pdestroy(u);
+   return f;
+}
-- 
cgit v1.2.3