1 files changed, 48 insertions, 0 deletions
diff --git a/src/benchmarks/cfrac/pimod.c b/src/benchmarks/cfrac/pimod.c
new file mode 100644
index 0000000..b26536d
--- /dev/null
+++ b/src/benchmarks/cfrac/pimod.c
@@ -0,0 +1,48 @@
+#include "pdefs.h"
+#include "precision.h"
+#ifdef ASM_16BIT
+#include "asm16bit.h"
+#endif
+
+/*
+ * Single-digit remainder
+ */
+int pimod(u, v)
+   register precision u;
+   int 	    v;
+{
+#ifndef ASM_16BIT
+   register digitPtr	uPtr;
+   register accumulator temp;		       /* 0 <= temp < base^2 */
+#endif
+   register digit	r = 0, d;		/* 0 <= r,d < base */
+   register int res = 0;
+
+   (void) pparm(u);
+   if (v < 0) d = (digit) -v; else d = (digit) v;
+   if (d >= BASE) {
+      errorp(PDOMAIN, "pimod", "divisor too big for single digit");
+      goto done;
+   } 
+   if (d == 0) {
+      errorp(PDOMAIN, "pimod", "divide by zero");
+      goto done;
+   } 
+#ifndef ASM_16BIT
+   uPtr = u->value + u->size;
+   r    = 0;			  	/* r is current remainder */
+   do {
+       temp  = mulBase(r);		/* 0 <= temp <= (base-1)^2 */
+       temp += *--uPtr;			/* 0 <= temp <= base(base-1) */
+       r     = temp % d; 		/* 0 <= r < base */
+   } while (uPtr > u->value);
+#else
+   r = memmodw1(u->value, u->size, d);
+#endif
+
+   res = (int) r;
+   if (u->sign) res = -res;
+done:
+   pdestroy(u);
+   return res;
+}