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

1 files changed, 44 insertions, 0 deletions

diff --git a/src/spheres_poly/sphere.h b/src/spheres_poly/sphere.h
new file mode 100644
index 0000000..28c64b4
--- /dev/null
+++ b/src/spheres_poly/sphere.h
@@ -0,0 +1,44 @@
+#ifndef  SPHERE_H
+#define  SPHERE_H
+
+
+#include "vector.h"
+
+
+class sphere {
+
+ public:
+
+  // constructor and destructor
+
+  sphere();
+  sphere(const sphere& s);
+  sphere(int i_i, vector<DIM> x, vector<DIM, int> cell_i, double lutime_i, 
+	 double r_i, double gr_i, double m_i, int species_i);
+  ~sphere();
+
+ //variables
+  
+  int i;                          // sphere ID
+
+  // impending event
+  event nextevent;                // next event...can be collision or transfer
+  event nextcollision;            // next collision if next event is transfer
+  // maybe nextnext event
+  
+  // past information
+  double lutime;                  // last update time
+  vector<DIM, int> cell;          // cell that it belongs to
+  vector<DIM, double> x;          // position
+  vector<DIM, double> v;          // velocity
+  double r;                       // sphere radius
+  double gr;                      // sphere growth rate
+  double m;                       // sphere mass
+  int species;                    // species number (not used during the MD)
+  // make sure efficent in memory
+
+ 
+
+};
+
+#endif