diff options
author | pants <jaron@kent-dobias.com> | 2016-09-08 11:54:06 -0400 |
---|---|---|
committer | pants <jaron@kent-dobias.com> | 2016-09-08 11:54:06 -0400 |
commit | b74bf87fd7eb2e75e71f6b667bd12b0ca7a31c11 (patch) | |
tree | c3fa46f10f6ae8604c4e784febc3eaa49b51b821 /src/update_boundary.c | |
parent | a264723b23d9ad6c3cb13450a6cf4b131f8d7338 (diff) | |
download | fuse_networks-b74bf87fd7eb2e75e71f6b667bd12b0ca7a31c11.tar.gz fuse_networks-b74bf87fd7eb2e75e71f6b667bd12b0ca7a31c11.tar.bz2 fuse_networks-b74bf87fd7eb2e75e71f6b667bd12b0ca7a31c11.zip |
cleaned up old files that don't matter
Diffstat (limited to 'src/update_boundary.c')
-rw-r--r-- | src/update_boundary.c | 86 |
1 files changed, 0 insertions, 86 deletions
diff --git a/src/update_boundary.c b/src/update_boundary.c deleted file mode 100644 index 6bc4fca..0000000 --- a/src/update_boundary.c +++ /dev/null @@ -1,86 +0,0 @@ - -#include "fracture.h" - -void update_boundary(net_t *instance, const double *avg_field) { - int size = instance->graph->ne; - int width = sqrt(size); - int num_verts = instance->graph->nv; - double *boundary = (double *)instance->boundary_cond->x; - - boundary[num_verts] = 0; - boundary[num_verts + 1] = 0; - - if (instance->voltage_bound) { - for (int i = 0; i < width / 4; i++) { - int t_ind = (width + 1) * (width / 8) + width / 4 - 1 + i; - double t_val = avg_field[t_ind]; - boundary[2 * i] = t_val; - boundary[2 * i + 1] = t_val; - - int b_ind = num_verts - (width + 1) * (width / 8) - width / 4 - i; - double b_val = avg_field[b_ind]; - boundary[num_verts - 1 - 2 * i] = b_val; - boundary[num_verts - 2 - 2 * i] = b_val; - - int l_ind = (width + 1) * (width / 8) + width / 2 + width / 4 - 1 + - i * (width + 1); - double l_val = avg_field[l_ind]; - boundary[width / 2 + 2 * i * (width + 1)] = l_val; - boundary[width / 2 + (2 * i + 1) * (width + 1)] = l_val; - - int r_ind = (width + 1) * (width / 8) + width - 1 + i * (width + 1); - double r_val = avg_field[r_ind]; - boundary[width + 2 * i * (width + 1)] = r_val; - boundary[width + (2 * i + 1) * (width + 1)] = r_val; - } - } else { - const double *stress = avg_field; - - for (int i = 0; i < width / 4; i++) { - boundary[2 * i] = stress[(width / 4 - 1) * width + width / 4 + 2 * i] + - stress[(width / 4 - 1) * width + width / 4 + 1 + 2 * i]; - boundary[2 * i + 1] = - stress[(width / 4 - 1) * width + width / 4 + 2 * i] + - stress[(width / 4 - 1) * width + width / 4 + 1 + 2 * i]; - - boundary[num_verts - 2 * i - 1] = - -stress[size - 1 - ((width / 4 - 1) * width + width / 4 + 2 * i)] - - stress[size - 1 - ((width / 4 - 1) * width + width / 4 + 1 + 2 * i)]; - boundary[num_verts - 2 * i - 2] = - -stress[size - 1 - ((width / 4 - 1) * width + width / 4 + 2 * i)] - - stress[size - 1 - ((width / 4 - 1) * width + width / 4 + 1 + 2 * i)]; - - boundary[(width + 1) / 2 + 2 * i * (width + 1)] = - stress[((width / 4) + 2 * i) * width + width / 4 - 1] - - stress[((width / 4) + 2 * i + 1) * width + width / 4 - 1]; - boundary[(width + 1) / 2 + (2 * i + 1) * (width + 1)] = - stress[((width / 4) + 2 * i) * width + width / 4 - 1] - - stress[((width / 4) + 2 * i + 1) * width + width / 4 - 1]; - - boundary[(width + 1) / 2 + 2 * i * (width + 1) + width / 2] = - stress[((width / 4) + 2 * i) * width + width / 4 + width / 2] - - stress[((width / 4) + 2 * i + 1) * width + width / 4 + width / 2]; - boundary[(width + 1) / 2 + (2 * i + 1) * (width + 1) + width / 2] = - stress[((width / 4) + 2 * i) * width + width / 4 + width / 2] - - stress[((width / 4) + 2 * i + 1) * width + width / 4 + width / 2]; - } - - double s_total = 0; - double total = 0; - for (int i = 0; i < num_verts; i++) { - if (boundary[i] != boundary[i]) { - boundary[i] = 0; - } else { - total += fabs(boundary[i]); - s_total += boundary[i]; - } - } - - // assert(fabs(s_total) < inf); - boundary[num_verts] -= s_total; - - for (int i = 0; i < num_verts + 1; i++) { - boundary[i] /= total; - } - } -} |