diff --git a/dune/mmdg/problems/mmdgproblem1.hh b/dune/mmdg/problems/mmdgproblem1.hh index c5efe4ba3d24aaa5107fb1e43dad3d1cce499197..6ea0cd89881315afa5d63d895939cc4ff5b7ba17 100644 --- a/dune/mmdg/problems/mmdgproblem1.hh +++ b/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; }