Do not pass end_points to setup_boundary

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;