Finding the lower left corner is easy

src/one-body-sample.cc

@@ -204,20 +204,9 @@ int main(int argc, char *argv[]) {
         coordinates[id] = it->geometry().corner(0);
         auto const &localCoordinates = coordinates[id];
 
-        // Find the center of the lower face
-        switch (dims) {
-          case 3:
-            if (localCoordinates[2] != lowerLeft[2])
-              break;
-          // fallthrough
-          case 2:
-            if (localCoordinates[1] == lowerLeft[1] &&
-                std::abs(localCoordinates[0] - lowerLeft[0]) < 1e-8)
-              specialNode = id;
-            break;
-          default:
-            assert(false);
-        }
+        // Find the lower left corner
+        if (localCoordinates == lowerLeft)
+          specialNode = id;
 
         // lower face
         if (localCoordinates[1] == lowerLeft[1]) {