diff --git a/src/finite-strain-elasticity.cc b/src/finite-strain-elasticity.cc
index 3a15baf2b8270ce26430e4a9fb606a09ad7fed29..d40c095f48ce2bb6a1e011d3dce028f482076638 100644
--- a/src/finite-strain-elasticity.cc
+++ b/src/finite-strain-elasticity.cc
@@ -363,11 +363,10 @@ int main (int argc, char *argv[]) try
     displacement -= identity;
 
     auto displacementFunction = Dune::Functions::makeDiscreteGlobalBasisFunction<FieldVector<double,dim>>(feBasis, displacement);
-    auto localDisplacementFunction = localFunction(displacementFunction);
 
     //  We need to subsample, because VTK cannot natively display real second-order functions
     SubsamplingVTKWriter<GridView> vtkWriter(gridView, Dune::refinementLevels(2));
-    vtkWriter.addVertexData(localDisplacementFunction, VTK::FieldInfo("displacement", VTK::FieldInfo::Type::vector, dim));
+    vtkWriter.addVertexData(displacementFunction, VTK::FieldInfo("displacement", VTK::FieldInfo::Type::vector, dim));
     vtkWriter.write(resultPath + "finite-strain_homotopy_" + std::to_string(i+1));
   }