From e76333883d3620622ac616b7ca68a8baf459c640 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 12 Jun 2019 10:34:51 -0400
Subject: fixed the current cutoff is the percolation measurements
---
src/perc_meas.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
(limited to 'src')
diff --git a/src/perc_meas.cpp b/src/perc_meas.cpp
index 23b46b3..e400982 100644
--- a/src/perc_meas.cpp
+++ b/src/perc_meas.cpp
@@ -147,7 +147,7 @@ void pm::bond_broken(const network& net, const current_info& cur, unsigned i) {
std::vector<bool> vertex_in(net.G.vertices.size());
for (unsigned i = 0; i < net.G.edges.size(); i++) {
- if (cur.currents[i] > 1.0 / pow(net.G.edges.size(), 2)) {
+ if (cur.currents[i] > CURRENT_CUTOFF) {
vertex_in[net.G.edges[i].v[0]] = true;
vertex_in[net.G.edges[i].v[1]] = true;
}
@@ -198,7 +198,7 @@ void pm::post_fracture(network &n) {
std::vector<bool> vertex_in(n.G.vertices.size());
for (unsigned i = 0; i < n.G.edges.size(); i++) {
- if (last_cur.currents[i] > 1.0 / n.G.edges.size()) {
+ if (last_cur.currents[i] > CURRENT_CUTOFF) {
vertex_in[n.G.edges[i].v[0]] = true;
vertex_in[n.G.edges[i].v[1]] = true;
}
--
cgit v1.2.3-70-g09d2