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/pgcd.c | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)
 create mode 100644 src/benchmarks/cfrac/pgcd.c

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

diff --git a/src/benchmarks/cfrac/pgcd.c b/src/benchmarks/cfrac/pgcd.c
new file mode 100644
index 0000000..a72a8a7
--- /dev/null
+++ b/src/benchmarks/cfrac/pgcd.c
@@ -0,0 +1,24 @@
+#include "precision.h"
+
+/*
+ * Euclid's Algorithm
+ *
+ * Given u and v, calculated and return their greatest common divisor.
+ */
+precision pgcd(u, v)
+   precision u, v;
+{
+   precision u3 = pnew(pabs(pparm(u))), v3 = pnew(pabs(pparm(v)));
+   precision q  = pUndef, r  = pUndef;
+
+   while (pnez(v3)) {
+      pdivmod(u3, v3, &q, &r);
+      pset(&u3, v3);
+      pset(&v3, r);
+   }
+
+   pdestroy(v3);
+   pdestroy(q);  pdestroy(r);  
+   pdestroy(u);  pdestroy(v);
+   return presult(u3);			/* result always positive */
+}
-- 
cgit v1.2.3