Skip to content

Commit

Permalink
simplifies cart quad export with face ordering
Browse files Browse the repository at this point in the history
  • Loading branch information
Mike Park authored and Mike Park committed Sep 29, 2022
1 parent 14168c6 commit a310dcf
Showing 1 changed file with 9 additions and 63 deletions.
72 changes: 9 additions & 63 deletions src/ref_oct.c
View file
Original file line number Diff line number Diff line change
Expand Up @@ -659,70 +659,16 @@ REF_FCN static REF_STATUS ref_oct_export_node(REF_OCT ref_oct, REF_INT node,
ref_cell = ref_grid_qua(ref_grid);
RSS(ref_oct_bbox_diag(bbox, &diag), "diag");

cell_face = 0; /* ymin */
if (ABS(bbox[2] - ref_oct->bbox[2]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
for (cell_face = 0; cell_face < 6; cell_face++) {
if (ABS(bbox[cell_face] - ref_oct->bbox[cell_face]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = cell_face + 1;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}
nodes[ref_cell_id_index(ref_cell)] = 3;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

cell_face = 1; /* xmax */
if (ABS(bbox[1] - ref_oct->bbox[1]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = 2;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

cell_face = 2; /* ymax */
if (ABS(bbox[3] - ref_oct->bbox[3]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = 4;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

cell_face = 3; /* xmin */
if (ABS(bbox[0] - ref_oct->bbox[0]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = 1;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

cell_face = 4; /* zmin */
if (ABS(bbox[4] - ref_oct->bbox[4]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = 5;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

cell_face = 5; /* zmax */
if (ABS(bbox[5] - ref_oct->bbox[5]) < tol * diag) {
for (face_node = 0; face_node < 4; face_node++) {
nodes[face_node] =
ref_cell_f2n_gen(ref_grid_hex(ref_grid), face_node, cell_face);
nodes[face_node] = ref_oct_c2n(ref_oct, nodes[face_node], node);
}
nodes[ref_cell_id_index(ref_cell)] = 6;
RSS(ref_cell_add(ref_cell, nodes, &new_cell), "add quad");
}

} else {
Expand Down

0 comments on commit a310dcf

Please sign in to comment.