summaryrefslogtreecommitdiff
path: root/src/fracture.c
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2017-01-14 11:59:39 -0500
committerJaron Kent-Dobias <jaron@kent-dobias.com>2017-01-14 11:59:39 -0500
commita66e8e2d7a1e9d5eb51f4954435513f6ecb8383c (patch)
treefd9fcbee58f5e3b22f718e856ebc2b1670a6a650 /src/fracture.c
parent3d6f2bbb756c7fed4e627da196e99047ec7406cc (diff)
downloadfuse_networks-a66e8e2d7a1e9d5eb51f4954435513f6ecb8383c.tar.gz
fuse_networks-a66e8e2d7a1e9d5eb51f4954435513f6ecb8383c.tar.bz2
fuse_networks-a66e8e2d7a1e9d5eb51f4954435513f6ecb8383c.zip
tried to fix voltage boundary conditions
Diffstat (limited to 'src/fracture.c')
-rw-r--r--src/fracture.c8
1 files changed, 8 insertions, 0 deletions
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);