From a43ff1f98e9b9814f858bccb11c174b418458491 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 10 Oct 2018 21:45:32 -0400
Subject: big rearrangement of files to make libraries and example (research)
 files clearer, and changed to c++ std lib random numbers

---
 lib/src/graph.cpp | 44 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 44 insertions(+)
 create mode 100644 lib/src/graph.cpp

(limited to 'lib/src/graph.cpp')

diff --git a/lib/src/graph.cpp b/lib/src/graph.cpp
new file mode 100644
index 0000000..4043413
--- /dev/null
+++ b/lib/src/graph.cpp
@@ -0,0 +1,44 @@
+
+#include <wolff/graph.hpp>
+
+graph_t::graph_t(D_t D, L_t L) {
+  nv = pow(L, D);
+  ne = D * nv;
+
+  v_adj.resize(nv);
+  coordinate.resize(nv);
+
+  for (std::vector<v_t> v_adj_i : v_adj) {
+    v_adj_i.reserve(2 * D);
+  }
+
+  for (v_t i = 0; i < nv; i++) {
+    coordinate[i].resize(D);
+    for (D_t j = 0; j < D; j++) {
+      coordinate[i][j] = (i / (v_t)pow(L, D - j - 1)) % L;
+
+      v_adj[i].push_back(pow(L, j + 1) * (i / ((v_t)pow(L, j + 1))) + fmod(i + pow(L, j), pow(L, j + 1)));
+      v_adj[i].push_back(pow(L, j + 1) * (i / ((v_t)pow(L, j + 1))) + fmod(pow(L, j+1) + i - pow(L, j), pow(L, j + 1)));
+    }
+  }
+}
+
+void graph_t::add_ext() {
+  for (std::vector<v_t>& v_adj_i : v_adj) {
+    v_adj_i.push_back(nv);
+  }
+
+  v_adj.resize(nv + 1);
+  coordinate.resize(nv + 1);
+  v_adj[nv].reserve(nv);
+
+  for (v_t i = 0; i < nv; i++) {
+    v_adj[nv].push_back(i);
+  }
+
+  coordinate[nv].resize(coordinate[0].size());
+
+  ne += nv;
+  nv += 1;
+}
+
-- 
cgit v1.2.3-70-g09d2