From 77043e249cf49fd533a1ffd6f53c0b6d6fcaaba8 Mon Sep 17 00:00:00 2001
From: Ian Jauslin <ian@jauslin.org>
Date: Thu, 26 May 2022 20:54:30 -0400
Subject: Make N be the smallest power of 2 larger than 3*K+1

---
 src/int_tools.c | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)
 create mode 100644 src/int_tools.c

(limited to 'src/int_tools.c')

diff --git a/src/int_tools.c b/src/int_tools.c
new file mode 100644
index 0000000..13e2cf9
--- /dev/null
+++ b/src/int_tools.c
@@ -0,0 +1,26 @@
+#include <math.h>
+#include "int_tools.h"
+
+// return smallest power of 2 that is > x
+int smallest_pow2(
+  int x
+){
+  return ipow(2,((int)log2(x)+1));
+}
+
+// integer power
+int ipow(
+  int x,
+  int n
+){
+  int out=1;
+  while (n>0)
+  {
+    if (n%2==1){
+      out*=x; 
+    }
+    n/=2;
+    x*=x;
+  }
+  return out;
+}
-- 
cgit v1.2.3-54-g00ecf