Skip to content
Snippets Groups Projects
Commit e775bd9f authored by Stéphane Del Pino's avatar Stéphane Del Pino
Browse files

Use m_cell_type to compute faces

parent 12805fd6
No related branches found
No related tags found
1 merge request!6Feature/crs
...@@ -14,8 +14,8 @@ void Connectivity<3>::_computeFaceCellConnectivities() ...@@ -14,8 +14,8 @@ void Connectivity<3>::_computeFaceCellConnectivities()
for (unsigned int j=0; j<this->numberOfCells(); ++j) { for (unsigned int j=0; j<this->numberOfCells(); ++j) {
const auto& cell_nodes = cell_to_node_matrix.rowConst(j); const auto& cell_nodes = cell_to_node_matrix.rowConst(j);
switch (cell_nodes.length) { switch (m_cell_type[j]) {
case 4: { // tetrahedron case CellType::Tetrahedron: {
cell_nb_faces[j] = 4; cell_nb_faces[j] = 4;
// face 0 // face 0
Face f0({cell_nodes(1), Face f0({cell_nodes(1),
...@@ -42,7 +42,7 @@ void Connectivity<3>::_computeFaceCellConnectivities() ...@@ -42,7 +42,7 @@ void Connectivity<3>::_computeFaceCellConnectivities()
face_cells_map[f3].emplace_back(std::make_tuple(j, 3, f3.reversed())); face_cells_map[f3].emplace_back(std::make_tuple(j, 3, f3.reversed()));
break; break;
} }
case 8: { // hexahedron case CellType::Hexahedron: {
// face 0 // face 0
Face f0({cell_nodes(3), Face f0({cell_nodes(3),
cell_nodes(2), cell_nodes(2),
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment