From 97b9dcd11bc64cfd31c5621749e121a52898e417 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Sun, 31 May 2020 14:46:10 +0200
Subject: [PATCH] add vtk output for aperture
---
dune/mmdg/mmdg.hh | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh
index a83541d..be04ec7 100644
--- a/dune/mmdg/mmdg.hh
+++ b/dune/mmdg/mmdg.hh
@@ -894,6 +894,7 @@ private:
//create the output using a linear interpolation for each interface element
Dune::DynamicVector<Scalar> iPressure(interfaceDOF, 0.0);
Dune::DynamicVector<Scalar> exactIPressure(interfaceDOF, 0.0);
+ Dune::DynamicVector<Scalar> aperture(interfaceDOF, 0.0);
for (const auto& iElem : elements(iGridView_))
{
@@ -914,6 +915,9 @@ private:
Base::problem_.exactInterfaceSolution(iGeo.corner(k));
}
+ aperture[iElemIdxOutput + k] =
+ Base::problem_.aperture(iGeo.corner(k));
+
//contribution of the basis function
// phi_iElem,0 (x) = indicator(iElem);
//at the kth corner of iElem
@@ -934,6 +938,8 @@ private:
vtkWriter(iGridView_, Dune::VTK::nonconforming);
vtkWriter.addVertexData( std::shared_ptr<const VTKFunction>(
new P1Function(iGridView_, iPressure, "interfacePressure")));
+ vtkWriter.addVertexData( std::shared_ptr<const VTKFunction>(
+ new P1Function(iGridView_, aperture, "aperture")));
if (Base::problem_.hasExactSolution())
{
--
GitLab