diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-12-05 19:16:28 -0500 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-12-05 19:16:28 -0500 |
commit | 4f4cf365eae07e04298459bf8f9e27ad0cfcc834 (patch) | |
tree | f9823c732c03d1ed0166c765a39c2a77d6f1afa9 /lib/src/network.cpp | |
parent | 66aad9c05bce441e3f74049b6c055312287e6a4a (diff) | |
download | fuse_networks-4f4cf365eae07e04298459bf8f9e27ad0cfcc834.tar.gz fuse_networks-4f4cf365eae07e04298459bf8f9e27ad0cfcc834.tar.bz2 fuse_networks-4f4cf365eae07e04298459bf8f9e27ad0cfcc834.zip |
now takes Lx and Ly for anisotropic fracture
Diffstat (limited to 'lib/src/network.cpp')
-rw-r--r-- | lib/src/network.cpp | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/lib/src/network.cpp b/lib/src/network.cpp index 9f6cb0d..bc7a0c6 100644 --- a/lib/src/network.cpp +++ b/lib/src/network.cpp @@ -7,7 +7,7 @@ network::network(const graph& G, cholmod_common *c) : c(c), G(G), fuses(G.edges. double v0y = G.vertices[G.edges[i][0]].r.y; double v1y = G.vertices[G.edges[i][1]].r.y; - if (fabs(v0y - v1y) > 0.5) { + if (fabs(v0y - v1y) > G.L.y / 2) { bool ind = v0y < v1y ? 0 : 1; ((double *)b->x)[G.edges[i][ind]] += 1.0; @@ -138,7 +138,7 @@ void network::break_edge(unsigned int e, bool unbreak) { double v0y = G.vertices[v0].r.y; double v1y = G.vertices[v1].r.y; - if (fabs(v0y - v1y) > 0.5) { + if (fabs(v0y - v1y) > G.L.y / 2) { bool ind = v0y < v1y ? unbreak : !unbreak; ((double *)b->x)[G.edges[e][ind]] -= 1.0; @@ -172,7 +172,7 @@ current_info network::get_current_info() { double v0y = G.vertices[G.edges[i][0]].r.y; double v1y = G.vertices[G.edges[i][1]].r.y; - if (fabs(v0y - v1y) > 0.5) { + if (fabs(v0y - v1y) > G.L.y / 2) { if (v0y > v1y) { currents[i] += 1.0; total_current += currents[i]; |