From 5c9909485667385951cf8294d7dad77c0944d8f9 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 22 May 2014 19:21:15 +0000
Subject: [PATCH] Write to VTK instead of AmiraMesh
[[Imported from SVN: r9751]]
---
harmonicmaps.cc | 17 +++++++++++------
1 file changed, 11 insertions(+), 6 deletions(-)
diff --git a/harmonicmaps.cc b/harmonicmaps.cc
index a4008adf..d37a1c3a 100644
--- a/harmonicmaps.cc
+++ b/harmonicmaps.cc
@@ -21,9 +21,10 @@
#include <dune/grid/utility/structuredgridfactory.hh>
#include <dune/grid/io/file/amirameshreader.hh>
-#include <dune/grid/io/file/amirameshwriter.hh>
+#include <dune/grid/io/file/vtk.hh>
#include <dune/fufem/boundarypatch.hh>
+#include <dune/fufem/functions/vtkbasisgridfunction.hh>
#include <dune/solvers/solvers/iterativesolver.hh>
#include <dune/solvers/norms/energynorm.hh>
@@ -268,7 +269,8 @@ int main (int argc, char *argv[]) try
// Output result
// //////////////////////////////
- BlockVector<FieldVector<double,3> > xEmbedded(x.size());
+ typedef BlockVector<FieldVector<double,3> > EmbeddedVectorType;
+ EmbeddedVectorType xEmbedded(x.size());
for (size_t i=0; i<x.size(); i++) {
#ifdef UNITVECTOR2
xEmbedded[i][0] = x[i].globalCoordinates()[0];
@@ -292,10 +294,13 @@ int main (int argc, char *argv[]) try
#endif
}
- LeafAmiraMeshWriter<GridType> amiramesh;
- amiramesh.addGrid(grid.leafGridView());
- amiramesh.addVertexData(xEmbedded, grid.leafGridView());
- amiramesh.write("resultGrid", 1);
+ VTKWriter<GridType::LeafGridView> vtkWriter(grid.leafGridView());
+ Dune::shared_ptr<VTKBasisGridFunction<FEBasis,EmbeddedVectorType> > vtkVectorField
+ = Dune::shared_ptr<VTKBasisGridFunction<FEBasis,EmbeddedVectorType> >
+ (new VTKBasisGridFunction<FEBasis,EmbeddedVectorType>(feBasis, xEmbedded, "orientation"));
+ vtkWriter.addVertexData(vtkVectorField);
+
+ vtkWriter.write("resultGrid");
// //////////////////////////////////////////////////////////
// Recompute and compare against exact solution
--
GitLab