Drop Neumann boundary

src/one-body-sample.cc

@@ -114,14 +114,23 @@ void setup_boundary(GridView const &gridView,
         assert(false);
     }
 
-    if (coordinates[1] == upperRight[1]) // upper face
-      ignoreNodes[id] = true;
-    else if (coordinates[1] == lowerLeft[1]) { // lower face
+    // lower face
+    if (coordinates[1] == lowerLeft[1]) {
       frictionalNodes[id] = true;
       ignoreNodes[id][1] = true;
-    } else if (coordinates[0] == lowerLeft[0]      // left face
-               || coordinates[0] == upperRight[0]) // right face
-      neumannNodes[id] = true;
+    }
+
+    // upper face
+    else if (coordinates[1] == upperRight[1])
+      ignoreNodes[id] = true;
+
+    // right face (except for both corners)
+    else if (coordinates[0] == upperRight[0])
+      ;
+
+    // left face (except for both corners)
+    else if (coordinates[0] == lowerLeft[0])
+      ;
   }
 }