[Cleanup] Move vertexCoordinates

src/one-body-sample.cc

@@ -154,17 +154,6 @@ int main(int argc, char *argv[]) {
       });
     }
 
-    Vector vertexCoordinates(fineVertexCount);
-    {
-      Dune::MultipleCodimMultipleGeomTypeMapper<
-          GridView, Dune::MCMGVertexLayout> const vertexMapper(leafView);
-      for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) {
-        auto const geometry = it->geometry();
-        assert(geometry.corners() == 1);
-        vertexCoordinates[vertexMapper.map(*it)] = geometry.corner(0);
-      }
-    }
-
     // Neumann boundary
     BoundaryPatch<GridView> const neumannBoundary(leafView);
 
@@ -340,6 +329,17 @@ int main(int argc, char *argv[]) {
     auto multigridStep = factory.getSolver();
 
     {
+      Vector vertexCoordinates(fineVertexCount);
+      {
+        Dune::MultipleCodimMultipleGeomTypeMapper<
+            GridView, Dune::MCMGVertexLayout> const vertexMapper(leafView);
+        for (auto it = leafView.begin<dims>(); it != leafView.end<dims>();
+             ++it) {
+          auto const geometry = it->geometry();
+          assert(geometry.corners() == 1);
+          vertexCoordinates[vertexMapper.map(*it)] = geometry.corner(0);
+        }
+      }
       std::fstream vertexCoordinateWriter("coordinates", std::fstream::out);
       for (size_t i = 0; i < fineVertexCount; ++i)
         if (frictionalNodes[i][0])