added support for hyperuniform lattices using existing hard-sphere jamming routine

1 files changed, 40 insertions, 0 deletions

diff --git a/src/spheres_poly/neighbor.cpp b/src/spheres_poly/neighbor.cpp
new file mode 100644
index 0000000..920c099
--- /dev/null
+++ b/src/spheres_poly/neighbor.cpp
@@ -0,0 +1,40 @@
+#include "vector.h"
+#include <iostream>
+#include "box.h"
+
+
+//==============================================================
+//==============================================================
+//  Class neighbor
+//==============================================================
+//==============================================================
+
+neighbor::neighbor(int i_i)
+ : i(i_i)
+{
+}
+
+
+//==============================================================
+//==============================================================
+//  Class collision
+//==============================================================
+//==============================================================
+
+collision::collision(int i_i, box *b_i)
+  : neighbor(i_i),
+    b(b_i)
+{
+  ctime = dblINF;
+  cpartner = i;
+}
+
+
+//==============================================================
+// Operation is finding the next collision from a given cell
+//==============================================================
+void collision::Operation(int j, vector<DIM, int>& pboffset)
+{
+  b->PredictCollision(i, j, pboffset, ctime, cpartner, cpartnerpboffset);
+}
+