diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 0cc55e51e09f54c94d8670967a46987d514bb30c..8002b890e22fd5594534ccc5cf976957d430ae5d 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -58,7 +58,6 @@ int const dim = 3; template <class GridView> void setup_boundary(GridView const &gridView, - Dune::FieldVector<double, dim> const &end_points, Dune::BitSetVector<dim> &ignoreNodes, Dune::BitSetVector<1> &neumannNodes, Dune::BitSetVector<1> &frictionalNodes) { @@ -75,7 +74,7 @@ void setup_boundary(GridView const &gridView, it != gridView.template end<dim>(); ++it) { assert(it->geometry().corners() == 1); Dune::FieldVector<double, dim> const coordinates = it->geometry().corner(0); - if (coordinates[1] == end_points[0]) { + if (coordinates[1] == 1) { ++dirichlet_nodes; size_t const id = myVertexMapper.map(*it); ignoreNodes[id] = true; @@ -84,7 +83,7 @@ void setup_boundary(GridView const &gridView, 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] == end_points[0]) { + } else if (coordinates[0] == 0 || coordinates[0] == 1) { ++neumann_nodes; size_t const id = myVertexMapper.map(*it); neumannNodes[id] = true; @@ -238,8 +237,7 @@ int main(int argc, char *argv[]) { false); Dune::BitSetVector<1> frictionalNodes(grid->size(grid->maxLevel(), dim), false); - setup_boundary(leafView, end_points, ignoreNodes, neumannNodes, - frictionalNodes); + setup_boundary(leafView, ignoreNodes, neumannNodes, frictionalNodes); typedef MyConvexProblem<OperatorType, VectorType> MyConvexProblemType; typedef MyBlockProblem<MyConvexProblemType> MyBlockProblemType;