diff --git a/dune/mmdg/problems/dgproblem.hh b/dune/mmdg/problems/dgproblem.hh
index b59b0a6fe05c457efd60ebb3a383b16ed45c3332..5f7e6307c6b7823e4107807ff636b73bc3d09a3f 100644
--- a/dune/mmdg/problems/dgproblem.hh
+++ b/dune/mmdg/problems/dgproblem.hh
@@ -56,7 +56,7 @@ class DGProblem
     }
 
     //returns the recommended quadrature order to compute an integral
-    //over x * boundary(x)
+    //over x * boundary(x) and K(x) * boundary(x)
     virtual int quadratureOrderBoundary () const
     {
       return 1;
diff --git a/dune/mmdg/problems/poisson.hh b/dune/mmdg/problems/poisson.hh
index 571d17613b2b6038adeb3df48fd884aa3cb845dc..549fb82c8610faeace62d4756158187aa09b73c8 100644
--- a/dune/mmdg/problems/poisson.hh
+++ b/dune/mmdg/problems/poisson.hh
@@ -8,6 +8,8 @@ template<class Coordinate, class Scalar = double>
 class Poisson : public DGProblem<Coordinate, Scalar>
 {
   public:
+    static constexpr int dim = Coordinate::dimension;
+    
     //the exact solution at position pos
     //1d: f(x)     = x,
     //2d: f(x,y)   = x*y,
@@ -29,6 +31,13 @@ class Poisson : public DGProblem<Coordinate, Scalar>
     {
       return true;
     };
+
+    //returns the recommended quadrature order to compute an integral
+    //over x * boundary(x) and K(x) * boundary(x)
+    int quadratureOrderBoundary () const
+    {
+      return dim;
+    }
 };
 
 #endif