Record coordinates

src/one-body-sample.cc

@@ -192,15 +192,17 @@ int main(int argc, char *argv[]) {
     Dune::BitSetVector<dims> ignoreNodes(finestSize, false);
     Dune::BitSetVector<1> neumannNodes(finestSize, false);
     Dune::BitSetVector<1> frictionalNodes(finestSize, false);
+
+    VectorType coordinates(finestSize);
     {
       Dune::MultipleCodimMultipleGeomTypeMapper<
           GridView, Dune::MCMGVertexLayout> const myVertexMapper(leafView);
 
       for (auto it = leafView.begin<dims>(); it != leafView.end<dims>(); ++it) {
         assert(it->geometry().corners() == 1);
-        Dune::FieldVector<double, dims> const localCoordinates =
-            it->geometry().corner(0);
         size_t const id = myVertexMapper.map(*it);
+        coordinates[id] = it->geometry().corner(0);
+        auto const &localCoordinates = coordinates[id];
 
         // Find the center of the lower face
         switch (dims) {