diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 75bf68df8da2caccfad86aed1c2822b97e63942b..41cf79d9a1cb7a05dbb3c88ad87cadcb3aaf90ae 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -253,7 +253,7 @@ int main (int argc, char *argv[]) try
 
     auto xFunction = Dune::Functions::makeDiscreteGlobalBasisFunction<TargetSpace::CoordinateType>(feBasis,TypeTree::hybridTreePath(),xEmbedded);
 
-    VTKWriter<GridType::LeafGridView> vtkWriter(grid->leafGridView());
+    SubsamplingVTKWriter<GridType::LeafGridView> vtkWriter(grid->leafGridView(),order-1);
     vtkWriter.addVertexData(xFunction, VTK::FieldInfo("orientation", VTK::FieldInfo::Type::vector, xEmbedded[0].size()));
     vtkWriter.write(resultPath + "_" + energy + "_result");