From 8afe36687544eae2bfaea3b9ec177bc589180034 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Maximilian=20H=C3=B6rl?=
<maximilian.hoerl@mathematik.uni-stuttgart.de>
Date: Sun, 22 Mar 2020 12:49:01 +0100
Subject: [PATCH] [bugfix] fix coupling index error
---
dune/mmdg/mmdg.hh | 24 ++++++++++++------------
1 file changed, 12 insertions(+), 12 deletions(-)
diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh
index c830287..1076ec9 100644
--- a/dune/mmdg/mmdg.hh
+++ b/dune/mmdg/mmdg.hh
@@ -350,10 +350,10 @@ private:
const Scalar couplingUpdate1 = quadratureFactor * betaEvaluation2;
- (*Base::A)[elemInIdxSLE][iElemIdx] -= couplingUpdate1;
- (*Base::A)[iElemIdx][elemInIdxSLE] -= couplingUpdate1;
- (*Base::A)[elemOutIdxSLE][iElemIdx] -= couplingUpdate1;
- (*Base::A)[iElemIdx][elemOutIdxSLE] -= couplingUpdate1;
+ (*Base::A)[elemInIdxSLE][iElemIdxSLE] -= couplingUpdate1;
+ (*Base::A)[iElemIdxSLE][elemInIdxSLE] -= couplingUpdate1;
+ (*Base::A)[elemOutIdxSLE][iElemIdxSLE] -= couplingUpdate1;
+ (*Base::A)[iElemIdxSLE][elemOutIdxSLE] -= couplingUpdate1;
for (int i = 0; i < dim - 1; i++)
{
@@ -361,10 +361,10 @@ private:
const Scalar couplingUpdate3 =
(iFrame[i] * qpGlobal) * couplingUpdate1;
- (*Base::A)[elemInIdxSLE + i + 1][iElemIdx] -= couplingUpdate2;
- (*Base::A)[iElemIdx][elemInIdxSLE + i + 1] -= couplingUpdate2;
- (*Base::A)[elemOutIdxSLE + i + 1][iElemIdx] -= couplingUpdate2;
- (*Base::A)[iElemIdx][elemOutIdxSLE + i + 1] -= couplingUpdate2;
+ (*Base::A)[elemInIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
+ (*Base::A)[iElemIdxSLE][elemInIdxSLE + i + 1] -= couplingUpdate2;
+ (*Base::A)[elemOutIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
+ (*Base::A)[iElemIdxSLE][elemOutIdxSLE + i + 1] -= couplingUpdate2;
(*Base::A)[elemInIdxSLE][iElemIdx + i + 1] -= couplingUpdate3;
(*Base::A)[iElemIdx + i + 1][elemInIdxSLE] -= couplingUpdate3;
@@ -374,10 +374,10 @@ private:
//i = dim - 1
const Scalar couplingUpdate4 = qpGlobal[dim] * couplingUpdate1;
- (*Base::A)[elemInIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4;
- (*Base::A)[iElemIdx][elemInIdxSLE + dim + 1] -= couplingUpdate4;
- (*Base::A)[elemOutIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4;
- (*Base::A)[iElemIdx][elemOutIdxSLE + dim + 1] -= couplingUpdate4;
+ (*Base::A)[elemInIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
+ (*Base::A)[iElemIdxSLE][elemInIdxSLE + dim + 1] -= couplingUpdate4;
+ (*Base::A)[elemOutIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
+ (*Base::A)[iElemIdxSLE][elemOutIdxSLE + dim + 1] -= couplingUpdate4;
}
// === interface entries ===
--
GitLab