1 files changed, 24 insertions, 0 deletions
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 */
+}