From 5957a1bd9cb5f895bf00dbefc689c71c30b777d3 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Fri, 27 Mar 2020 13:22:02 +0100
Subject: [PATCH] shift evaluation for vtk output to correctly picture
discontinuities at the interface
---
dune/mmdg/dg.hh | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/dune/mmdg/dg.hh b/dune/mmdg/dg.hh
index 344c397..7ae4a82 100644
--- a/dune/mmdg/dg.hh
+++ b/dune/mmdg/dg.hh
@@ -524,8 +524,10 @@ protected:
for (int k = 0; k < geo.corners(); k++)
{
if (problem_.hasExactSolution())
- {
- exactPressure[elemIdxSLE + k] = problem_.exactSolution(geo.corner(k));
+ { //evaluation at the corners is slightly shifted towards the center of
+ //the element to correctly picture discontinuities of the solution
+ exactPressure[elemIdxSLE + k] = problem_.exactSolution(
+ 0.9999 * geo.corner(k) + 0.0001 * geo.center());
}
//contribution of the basis function
--
GitLab