added support for rotating a lattice

2 files changed, 11 insertions, 4 deletions

diff --git a/src/sample.cpp b/src/sample.cpp
index b8af937..488d20d 100644
--- a/src/sample.cpp
+++ b/src/sample.cpp
@@ -16,25 +16,25 @@ void sample::pre_fracture(const network& n) {
   sample_file << "<|";
   sample_file << "\"vr\"->{";
   for (const graph::vertex &v : n.G.vertices) {
-    sample_file << "{" << v.r.x << "," << v.r.y << "},";
+    sample_file <<std::fixed<< "{" << v.r.x << "," << v.r.y << "},";
   }
   sample_file << "},\"vp\"->{";
   for (const graph::vertex &v : n.G.vertices) {
     sample_file << "{";
     for (const graph::coordinate &r : v.polygon) {
-      sample_file << "{" << r.x << "," << r.y << "},";
+      sample_file <<std::fixed<< "{" << r.x << "," << r.y << "},";
     }
     sample_file << "},";
   }
   sample_file << "},\"ur\"->{";
   for (const graph::vertex &v : n.G.dual_vertices) {
-    sample_file << "{" << v.r.x << "," << v.r.y << "},";
+    sample_file <<std::fixed<< "{" << v.r.x << "," << v.r.y << "},";
   }
   sample_file << "},\"up\"->{";
   for (const graph::vertex &v : n.G.dual_vertices) {
     sample_file << "{";
     for (const graph::coordinate &r : v.polygon) {
-      sample_file << "{" << r.x << "," << r.y << "},";
+      sample_file <<std::fixed<< "{" << r.x << "," << r.y << "},";
     }
     sample_file << "},";
   }
diff --git a/src/sample_fracture.cpp b/src/sample_fracture.cpp
index 42a6a5a..782b3e4 100644
--- a/src/sample_fracture.cpp
+++ b/src/sample_fracture.cpp
@@ -47,6 +47,7 @@ int main(int argc, char* argv[]) {
   CHOL_F(start)(&c);
 
   sample meas(Lx, Ly, beta);
+  //sample meas2(Ly, Lx, beta);
 
   randutils::auto_seed_128 seeds;
   std::mt19937 rng{seeds};
@@ -56,6 +57,12 @@ int main(int argc, char* argv[]) {
     network network(G, &c);
     network.set_thresholds(beta, rng);
     network.fracture(meas);
+
+    /*graph G2 = G.rotate();
+    class network network2(G2, &c);
+    network2.thresholds = network.thresholds;
+    network2.fracture(meas2);
+    */
   }
 
   CHOL_F(finish)(&c);