Skip to content
Snippets Groups Projects
Commit 8afe3668 authored by Hörl, Maximilian's avatar Hörl, Maximilian
Browse files

[bugfix] fix coupling index error

parent 304a0bf3
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 12 additions and 12 deletions
dune/mmdg/mmdg.hh
+ 12
− 12
View file @ 8afe3668
...@@ -350,10 +350,10 @@ private: ...@@ -350,10 +350,10 @@ private:
const Scalar couplingUpdate1 = quadratureFactor * betaEvaluation2; const Scalar couplingUpdate1 = quadratureFactor * betaEvaluation2;
(*Base::A)[elemInIdxSLE][iElemIdx] -= couplingUpdate1; (*Base::A)[elemInIdxSLE][iElemIdxSLE] -= couplingUpdate1;
(*Base::A)[iElemIdx][elemInIdxSLE] -= couplingUpdate1; (*Base::A)[iElemIdxSLE][elemInIdxSLE] -= couplingUpdate1;
(*Base::A)[elemOutIdxSLE][iElemIdx] -= couplingUpdate1; (*Base::A)[elemOutIdxSLE][iElemIdxSLE] -= couplingUpdate1;
(*Base::A)[iElemIdx][elemOutIdxSLE] -= couplingUpdate1; (*Base::A)[iElemIdxSLE][elemOutIdxSLE] -= couplingUpdate1;
for (int i = 0; i < dim - 1; i++) for (int i = 0; i < dim - 1; i++)
{ {
...@@ -361,10 +361,10 @@ private: ...@@ -361,10 +361,10 @@ private:
const Scalar couplingUpdate3 = const Scalar couplingUpdate3 =
(iFrame[i] * qpGlobal) * couplingUpdate1; (iFrame[i] * qpGlobal) * couplingUpdate1;
(*Base::A)[elemInIdxSLE + i + 1][iElemIdx] -= couplingUpdate2; (*Base::A)[elemInIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
(*Base::A)[iElemIdx][elemInIdxSLE + i + 1] -= couplingUpdate2; (*Base::A)[iElemIdxSLE][elemInIdxSLE + i + 1] -= couplingUpdate2;
(*Base::A)[elemOutIdxSLE + i + 1][iElemIdx] -= couplingUpdate2; (*Base::A)[elemOutIdxSLE + i + 1][iElemIdxSLE] -= couplingUpdate2;
(*Base::A)[iElemIdx][elemOutIdxSLE + i + 1] -= couplingUpdate2; (*Base::A)[iElemIdxSLE][elemOutIdxSLE + i + 1] -= couplingUpdate2;
(*Base::A)[elemInIdxSLE][iElemIdx + i + 1] -= couplingUpdate3; (*Base::A)[elemInIdxSLE][iElemIdx + i + 1] -= couplingUpdate3;
(*Base::A)[iElemIdx + i + 1][elemInIdxSLE] -= couplingUpdate3; (*Base::A)[iElemIdx + i + 1][elemInIdxSLE] -= couplingUpdate3;
...@@ -374,10 +374,10 @@ private: ...@@ -374,10 +374,10 @@ private:
//i = dim - 1 //i = dim - 1
const Scalar couplingUpdate4 = qpGlobal[dim] * couplingUpdate1; const Scalar couplingUpdate4 = qpGlobal[dim] * couplingUpdate1;
(*Base::A)[elemInIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4; (*Base::A)[elemInIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
(*Base::A)[iElemIdx][elemInIdxSLE + dim + 1] -= couplingUpdate4; (*Base::A)[iElemIdxSLE][elemInIdxSLE + dim + 1] -= couplingUpdate4;
(*Base::A)[elemOutIdxSLE + dim + 1][iElemIdx] -= couplingUpdate4; (*Base::A)[elemOutIdxSLE + dim + 1][iElemIdxSLE] -= couplingUpdate4;
(*Base::A)[iElemIdx][elemOutIdxSLE + dim + 1] -= couplingUpdate4; (*Base::A)[iElemIdxSLE][elemOutIdxSLE + dim + 1] -= couplingUpdate4;
} }
// === interface entries === // === interface entries ===
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment