Comments

src/one-body-sample.cc

@@ -95,13 +95,13 @@ void setup_boundary(GridView const &gridView,
     assert(it->geometry().corners() == 1);
     Dune::FieldVector<double, dim> const coordinates = it->geometry().corner(0);
     size_t const id = myVertexMapper.map(*it);
-    if (coordinates[1] == upperRight[1])
+    if (coordinates[1] == upperRight[1]) // upper face
       ignoreNodes[id] = true;
-    else if (coordinates[1] == lowerLeft[1]) {
+    else if (coordinates[1] == lowerLeft[1]) { // lower face
       frictionalNodes[id] = true;
-      ignoreNodes[id][1] = true; // Zero displacement in direction y
-    } else if (coordinates[0] == lowerLeft[0] ||
-               coordinates[0] == upperRight[0])
+      ignoreNodes[id][1] = true;
+    } else if (coordinates[0] == lowerLeft[0]      // left face
+               || coordinates[0] == upperRight[0]) // right face
       neumannNodes[id] = true;
   }
 }