From 3a13fa257a7589befafced123a4bfd02320b1e48 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Mon, 30 Mar 2020 16:46:35 +0200
Subject: [PATCH] simplify problem mmdg1
---
dune/mmdg/problems/mmdgproblem1.hh | 35 +++---------------------------
1 file changed, 3 insertions(+), 32 deletions(-)
diff --git a/dune/mmdg/problems/mmdgproblem1.hh b/dune/mmdg/problems/mmdgproblem1.hh
index c5efe4b..6ea0cd8 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;
}
--
GitLab