From ae6ad8569615d81fd4b4a8f13318c8f90a768a37 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 27 May 2019 22:26:21 -0400 Subject: refactored much of the fracture library and fixed some percolation measurements --- src/analysis_tools.hpp | 2 ++ src/perc_meas.cpp | 4 +++- src/percolation.cpp | 5 ++++- 3 files changed, 9 insertions(+), 2 deletions(-) (limited to 'src') diff --git a/src/analysis_tools.hpp b/src/analysis_tools.hpp index 3beb1a8..659cacf 100644 --- a/src/analysis_tools.hpp +++ b/src/analysis_tools.hpp @@ -18,6 +18,8 @@ struct EdgeProperties { }; typedef boost::adjacency_list Graph; +typedef boost::graph_traits::vertex_descriptor Vertex; +typedef boost::graph_traits::vertices_size_type VertexIndex; template diff --git a/src/perc_meas.cpp b/src/perc_meas.cpp index 4220a93..ec7746e 100644 --- a/src/perc_meas.cpp +++ b/src/perc_meas.cpp @@ -159,7 +159,9 @@ void pm::post_fracture(network &n) { } else { ss[components[i].size() - 1]++; } - sn[p][components[i].size() - 1]++; + for (unsigned j = p; j < sn.size(); j++) { + sn[j][components[i].size() - 1]++; + } } sp[p - 1]++; diff --git a/src/percolation.cpp b/src/percolation.cpp index 03da117..757ab44 100644 --- a/src/percolation.cpp +++ b/src/percolation.cpp @@ -76,20 +76,23 @@ int main(int argc, char* argv[]) { pm meas(n, a); for (unsigned trial = 0; trial < N; trial++) { + /* while (true) { try { + */ graph G(n, a, rng); percolation_network fuse_network(G, &c); fuse_network.set_thresholds(beta, rng); fuse_network.fracture(meas); + /* break; } catch (std::exception &e) { std::cout << e.what() << '\n'; } } - if (quit.load()) break; + */ } } -- cgit v1.2.3-70-g09d2