From 45cbed8d1fe61bd704a0c2f4d67cf830b7e44945 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 22 May 2014 19:43:55 +0000
Subject: [PATCH] Temporary: Add the z-coordinate of the deformation as a
 scalar field

This makes visualizing wrinkles much easier.

[[Imported from SVN: r9753]]
---
 dune/gfe/cosseratvtkwriter.hh | 9 +++++++++
 1 file changed, 9 insertions(+)

diff --git a/dune/gfe/cosseratvtkwriter.hh b/dune/gfe/cosseratvtkwriter.hh
index 577b0db8..ef03a8c5 100644
--- a/dune/gfe/cosseratvtkwriter.hh
+++ b/dune/gfe/cosseratvtkwriter.hh
@@ -6,6 +6,7 @@
 
 #include <dune/fufem/functionspacebases/p1nodalbasis.hh>
 #include <dune/fufem/functions/vtkbasisgridfunction.hh>
+#include <dune/fufem/functiontools/basisinterpolator.hh>
 #include <dune/gfe/rigidbodymotion.hh>
 
 
@@ -186,6 +187,14 @@ public:
             vtkWriter.addVertexData(vtkDirector);
         }
 
+        // For easier visualization of wrinkles: add z-coordinate as scalar field
+        std::vector<double> zCoord(downsampledConfig.size());
+        for (size_t i=0; i<zCoord.size(); i++)
+          zCoord[i] = downsampledConfig[i].r[2];
+
+        vtkWriter.addVertexData(zCoord, "zCoord");
+
+        // Write the file to disk
         vtkWriter.write(filename);
 
     }
-- 
GitLab