Use the SubsamplingVTKWriter

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");