added preliminary support for gathering statistics at crack growth, not critical stress

1 files changed, 4 insertions, 2 deletions

diff --git a/src/net.c b/src/net.c
index 9f6965b..f3aeda9 100644
--- a/src/net.c
+++ b/src/net.c
@@ -57,8 +57,8 @@ net_t *net_create(const graph_t *g, double inf, double beta, double notch_len, b
 	net->voltage_bound = vb;
 	net->boundary_cond = bound_set(g, vb, notch_len, c);
 
-	if (g->boundary != TORUS_BOUND) net->adjacency = gen_adjacency(net, false, false, 0, c);
-	else net->adjacency = gen_adjacency(net, true, false, 0, c);
+	net->adjacency = gen_adjacency(net, false, false, 0, c);
+	net->dual_adjacency = gen_adjacency(net, true, false, 0, c);
 
 	net->marks = (uint_t *)malloc((net->graph->break_dim) * sizeof(uint_t));
 	assert(net->marks != NULL);
@@ -112,6 +112,7 @@ net_t *net_copy(const net_t *net, cholmod_common *c) {
 	memcpy(net_copy->dual_marks, net->dual_marks, dual_marks_size);
 
 	net_copy->adjacency = CHOL_F(copy_sparse)(net->adjacency, c);
+	net_copy->dual_adjacency = CHOL_F(copy_sparse)(net->dual_adjacency, c);
 	net_copy->boundary_cond = CHOL_F(copy_dense)(net->boundary_cond, c);
 	net_copy->factor = CHOL_F(copy_factor)(net->factor, c);
 
@@ -123,6 +124,7 @@ void net_free(net_t *net, cholmod_common *c) {
 	free(net->thres);
 	CHOL_F(free_dense)(&(net->boundary_cond), c);
 	CHOL_F(free_sparse)(&(net->adjacency), c);
+	CHOL_F(free_sparse)(&(net->dual_adjacency), c);
 	CHOL_F(free_factor)(&(net->factor), c);
 	free(net->marks);
 	free(net->dual_marks);