From ee5b4765809a213cb11ffdb549d777e81b0b94af Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 18 Mar 2015 15:09:34 +0000
Subject: [PATCH] Port from the dune-fufem to the dune-functions mechanism for
writing VTK files
[[Imported from SVN: r10083]]
---
src/harmonicmaps.cc | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 69f64917..05a1bf07 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -21,11 +21,11 @@
#include <dune/grid/io/file/gmshreader.hh>
#include <dune/grid/io/file/vtk.hh>
+#include <dune/functions/gridfunctions/discretescalarglobalbasisfunction.hh>
#include <dune/functions/functionspacebases/pq1nodalbasis.hh>
#include <dune/functions/functionspacebases/pqknodalbasis.hh>
#include <dune/fufem/boundarypatch.hh>
-#include <dune/fufem/functions/vtkbasisgridfunction.hh>
#include <dune/fufem/functiontools/basisinterpolator.hh>
#include <dune/fufem/functiontools/boundarydofs.hh>
#include <dune/fufem/functionspacebases/dunefunctionsbasis.hh>
@@ -240,10 +240,11 @@ int main (int argc, char *argv[]) try
xEmbedded[i][2] = x[i].globalCoordinates()[2];
}
- VTKWriter<GridType::LeafGridView> vtkWriter(grid->leafGridView());
- auto vtkVectorField = std::make_shared<VTKBasisGridFunction<FufemFEBasis,EmbeddedVectorType> >(fufemFeBasis, xEmbedded, "orientation");
- vtkWriter.addVertexData(vtkVectorField);
+ Functions::DiscreteScalarGlobalBasisFunction<decltype(feBasis),decltype(xEmbedded)> xFunction(feBasis,xEmbedded);
+ auto localXFunction = localFunction(xFunction);
+ VTKWriter<GridType::LeafGridView> vtkWriter(grid->leafGridView());
+ vtkWriter.addVertexData(localXFunction, VTK::FieldInfo("orientation", VTK::FieldInfo::Type::vector, 3));
vtkWriter.write(resultPath + "_" + energy + "_result");
/////////////////////////////////////////////////////////////////
--
GitLab