[407-fix-deformed-dg-plotting] Fixed deforming grid

src/lib4_fem/plot/VtuPlot.cpp

@@ -614,7 +614,7 @@ void VtuPlot::deformDGGrid(const Vector &displacement) {
       DGVectorFieldElement E(displacement, c);
       for (int i = 0; i < c.Corners(); ++i) {
         exBuffer.Send(0) << c.Corner(i);
-        VectorField u = E.VectorValue(c.Corner(i), displacement);
+        VectorField u = E.VectorValue(c.LocalCorner(i), displacement);
         for (int j = 0; j < dim; ++j)
           exBuffer.Send(0) << u[j];
       }
@@ -643,7 +643,7 @@ void VtuPlot::deformDGGrid(const Vector &displacement) {
     }
 
 //  for (auto &p: points) {
-  for (int i = 0; i < points.size(); i++) {
+  for (int i = 0; i < points[0.0].size(); i++) {
     displacements[0.0][i] =
       deformationData[points[0.0][i]] / deformationCount[points[0.0][i]];
   }
@@ -678,7 +678,7 @@ void VtuPlot::deformLagrangeGrid(const Vector &displacement) {
       deformationData->insert({z, d});
     }
 
-  for (int i = 0; i < points.size(); i++) {
+  for (int i = 0; i < points[0.0].size(); i++) {
     displacements[0.0][i] = deformationData->at(points[0.0].at(i));
   }