simplify problem mmdg1

dune/mmdg/problems/mmdgproblem1.hh

@@ -19,22 +19,13 @@ public:
   //the exact bulk solution at position pos
   Scalar exactSolution (const Coordinate& pos) const
   {
-    return (pos[0] < 0.5) ? pos[0] / (1.0 + d_) : (pos[0] - 1.0) / (1.0 + d_) + 1.0;
-    Scalar solution = pos * Coordinate(1.0);
-
-    if (pos[0] > 0.5)
-    {
-      solution += 1.0;
-    }
-
-    return solution;
+    return (pos[0] < 0.5) ? pos[0] : pos[0] + 1;
   }
 
   //the exact solution on the interface at position pos
   Scalar exactInterfaceSolution (const Coordinate& pos) const
   {
-    return 0.5;
-    return pos * Coordinate(1.0) + 0.5;
+    return 1.0;
   }
 
   //indicates whether an exact solution is implemented for the problem
@@ -49,26 +40,6 @@ public:
     return d_;
   }
 
-  //tangential permeability tensor of the interface at position pos
-  Matrix Kparallel (const Coordinate& pos) const
-  {
-    Matrix permeability(0.0);
-
-    for (int i = 0; i < dim; i++)
-    {
-      permeability[i][i] = 1.0 / d_;
-    }
-
-    //return identity matrix
-    return permeability;
-  }
-
-  //permeability of the fracture in normal direction at position pos
-  Scalar Kperp (const Coordinate& pos) const
-  {
-    return Scalar(1.0) / d_;
-  }
-
   //returns the recommended quadrature order to compute an integral
   //over x * boundary(x)
   int quadratureOrderBoundary () const
@@ -79,7 +50,7 @@ public:
   //returns the recommended quadrature order to compute an integral
   //over x * interfaceBoundary(x)
   int quadratureOrderInterfaceBoundary () const
-  { //NOTE: only relevant for 3d, not implemented!
+  {
     return 2;
   }