From a66e8e2d7a1e9d5eb51f4954435513f6ecb8383c Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Sat, 14 Jan 2017 11:59:39 -0500 Subject: tried to fix voltage boundary conditions --- src/fracture.c | 8 ++++++++ 1 file changed, 8 insertions(+) (limited to 'src/fracture.c') diff --git a/src/fracture.c b/src/fracture.c index bb7701b..af9cd0f 100644 --- a/src/fracture.c +++ b/src/fracture.c @@ -325,6 +325,14 @@ int main(int argc, char *argv[]) { } } + FILE *testout = fopen("test.txt", "w"); + double *tmp_voltage = net_voltages(net, &c); + for (uint_t j = 0; j < g->bound_inds[1]; j++) { + fprintf(testout, "%g ", tmp_voltage[g->bound_verts[j]]-tmp_voltage[g->nv+j]); + } + fclose(testout); + free(tmp_voltage); + if (save_damage) { uint_t would_break = 0; double *tmp_voltage = net_voltages(net, &c); -- cgit v1.2.3-70-g09d2