diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 2ecfbe7d7d58a12bb84e5065b93efd2ca00ae381..52142ad4431229c1ad38377c7dfd764747423e0f 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -252,10 +252,9 @@ int main (int argc, char *argv[]) try
         xEmbedded[i] = x[i].globalCoordinates();
 
     auto xFunction = Dune::Functions::makeDiscreteGlobalBasisFunction<TargetSpace::CoordinateType>(feBasis,TypeTree::hybridTreePath(),xEmbedded);
-    auto localXFunction = localFunction(xFunction);
 
     VTKWriter<GridType::LeafGridView> vtkWriter(grid->leafGridView());
-    vtkWriter.addVertexData(localXFunction, VTK::FieldInfo("orientation", VTK::FieldInfo::Type::scalar, xEmbedded[0].size()));
+    vtkWriter.addVertexData(xFunction, VTK::FieldInfo("orientation", VTK::FieldInfo::Type::scalar, xEmbedded[0].size()));
     vtkWriter.write(resultPath + "_" + energy + "_result");
 
     // Write the corresponding coefficient vector: verbatim in binary, to be completely lossless