diff --git a/dune/mmdg/dg.hh b/dune/mmdg/dg.hh index 344c397e7575d645a5d17096a14080e5d0b9b3d2..7ae4a826f980642a79f9d4b6c9c1056e33406e4c 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