From a73445a1235047fa81c001a1cd642445b5df4e1f Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 4 Nov 2011 12:21:24 +0100
Subject: [PATCH] Write displacement to a vtk file
---
src/one-body-sample.cc | 12 ++++++++++++
1 file changed, 12 insertions(+)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index cfe1c64f..faf33e22 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -7,6 +7,8 @@
#include <exception>
#include <iostream>
+#include <dune/common/shared_ptr.hh>
+
#include <dune/common/exceptions.hh>
#include <dune/common/stdstreams.hh>
@@ -26,6 +28,7 @@
#include <dune/fufem/boundarypatch.hh>
#include <dune/fufem/functions/constantfunction.hh>
#include <dune/fufem/functionspacebases/p1nodalbasis.hh>
+#include <dune/grid/io/file/vtk/vtkwriter.hh>
#include <dune/tnnmg/iterationsteps/genericnonlineargs.hh>
@@ -34,6 +37,8 @@
#include <dune/solvers/norms/energynorm.hh>
#include <dune/solvers/solvers/loopsolver.hh>
+#include <dune/vtkgridfunction.hh> // currently homeless
+
#include <dune/tectonic/globalnonlinearity.hh>
#include <dune/tectonic/myconvexproblem.hh>
#include <dune/tectonic/myblockproblem.hh>
@@ -237,6 +242,13 @@ int main() {
}
}
+ Dune::shared_ptr<Dune::VTKBasisGridFunction<P1Basis, VectorType>> ptr(
+ new Dune::VTKBasisGridFunction<P1Basis, VectorType>(p1Basis, u1,
+ "displacement"));
+ Dune::VTKWriter<GridType::LeafGridView> writer(leafView);
+ writer.addVertexData(ptr);
+ writer.write("weeee");
+
VectorType diff = u2;
diff -= u1;
std::cout << "Infinity norm of the difference of the two solutions: "
--
GitLab