Skip to content

Commit

Permalink
adds iso interp test
Browse files Browse the repository at this point in the history
  • Loading branch information
Mike Park authored and Mike Park committed Aug 30, 2021
1 parent b999dfa commit 089dea5
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 1 deletion.
1 change: 1 addition & 0 deletions src/ref_iso.c
Expand Up @@ -133,6 +133,7 @@ static REF_STATUS ref_iso_interp(REF_GRID iso_grid, REF_EDGE ref_edge,
}

ref_free(edge_real);

return REF_SUCCESS;
}

Expand Down
8 changes: 7 additions & 1 deletion src/ref_iso_test.c
Expand Up @@ -150,6 +150,8 @@ int main(int argc, char *argv[]) {
REF_GRID ref_grid, iso_grid;
REF_NODE ref_node;
REF_DBL *field;
REF_DBL *out;
REF_INT node;

RSS(ref_fixture_tet_grid(&ref_grid, ref_mpi), "tri");
ref_node = ref_grid_node(ref_grid);
Expand All @@ -158,12 +160,16 @@ int main(int argc, char *argv[]) {
field[1] = -1;
field[2] = -1;
field[3] = 1;
RSS(ref_iso_insert(&iso_grid, ref_grid, field, 0, NULL, NULL), "iso");
RSS(ref_iso_insert(&iso_grid, ref_grid, field, 1, field, &out), "iso");
if (!ref_mpi_para(ref_mpi)) {
REIS(4, ref_node_n(ref_grid_node(iso_grid)), "three nodes");
REIS(2, ref_cell_n(ref_grid_tri(iso_grid)), "one tri");
REIS(0, ref_cell_n(ref_grid_edg(iso_grid)), "no edg");
}
each_ref_node_valid_node(ref_grid_node(iso_grid), node) {
RWDS(0.0, out[node], -1, "interp");
}
ref_free(out);
ref_grid_free(iso_grid);
ref_free(field);
ref_grid_free(ref_grid);
Expand Down

0 comments on commit 089dea5

Please sign in to comment.