1 files changed, 28 insertions, 0 deletions
diff --git a/src/benchmarks/cfrac/ppowmod.c b/src/benchmarks/cfrac/ppowmod.c
new file mode 100644
index 0000000..4528db9
--- /dev/null
+++ b/src/benchmarks/cfrac/ppowmod.c
@@ -0,0 +1,28 @@
+#include "precision.h"
+
+/*
+ * Raise to precision power mod m
+ */
+precision ppowmod(u, v, m)
+   precision 	u, v, m;
+{
+   precision j = pUndef, i = pUndef, n = pUndef;
+
+   (void) pparm(m);
+   pset(&i,   pparm(u));
+   pset(&n,   pparm(v));
+   pset(&j,   pone);
+
+   do {
+      if (podd(n)) {
+	 pset(&j, pmod(pmul(i, j), m));
+      }
+      pset(&n, phalf(n));
+      if (peqz(n)) break;
+      pset(&i, pmod(pmul(i, i), m));
+   } while (1);
+
+   pdestroy(i); pdestroy(n);
+   pdestroy(u); pdestroy(v); pdestroy(m);
+   return presult(j);
+}