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

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

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;
+}
-- 
cgit v1.2.3