diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index f7dedf7750f30841c45b1b79b00fbdba2e11c4b7..a815fb4b38573cc977d005ba4e83595796a216d9 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -84,24 +84,18 @@ void setup_boundary(GridView const &gridView, Dune::MultipleCodimMultipleGeomTypeMapper< GridView, Dune::MCMGVertexLayout> const myVertexMapper(gridView); - size_t dirichlet_nodes = 0; - size_t neumann_nodes = 0; - size_t frictional_nodes = 0; for (auto it = gridView.template begin<dim>(); it != gridView.template end<dim>(); ++it) { assert(it->geometry().corners() == 1); Dune::FieldVector<double, dim> const coordinates = it->geometry().corner(0); if (coordinates[1] == 1) { - ++dirichlet_nodes; size_t const id = myVertexMapper.map(*it); ignoreNodes[id] = true; } else if (coordinates[1] == 0) { - ++frictional_nodes; size_t const id = myVertexMapper.map(*it); frictionalNodes[id] = true; ignoreNodes[id][1] = true; // Zero displacement in direction y } else if (coordinates[0] == 0 || coordinates[0] == 1) { - ++neumann_nodes; size_t const id = myVertexMapper.map(*it); neumannNodes[id] = true; }