From 2c2b5249f02fc08f15576bcd0c02198c807193bd Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Wed, 9 Sep 2020 08:06:02 +0200
Subject: [PATCH] [md][cclocalassembler] Use undeflected privars for setting
internal Dirichlet
---
dumux/multidomain/subdomaincclocalassembler.hh | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dumux/multidomain/subdomaincclocalassembler.hh b/dumux/multidomain/subdomaincclocalassembler.hh
index a4657ac225..12bd85cde5 100644
--- a/dumux/multidomain/subdomaincclocalassembler.hh
+++ b/dumux/multidomain/subdomaincclocalassembler.hh
@@ -415,7 +415,7 @@ public:
{
if (internalDirichletConstraintsOwnElement[eqIdx])
{
- origResiduals[0][eqIdx] = this->curElemVolVars()[scv].priVars()[eqIdx] - dirichletValues[eqIdx];
+ origResiduals[0][eqIdx] = origVolVars.priVars()[eqIdx] - dirichletValues[eqIdx];
A[globalI][globalI][eqIdx][pvIdx] = (eqIdx == pvIdx) ? 1.0 : 0.0;
}
else
@@ -717,7 +717,7 @@ public:
{
if (internalDirichletConstraints[eqIdx])
{
- residual[eqIdx] = this->curElemVolVars()[scv].priVars()[eqIdx] - dirichletValues[eqIdx];
+ residual[eqIdx] = origVolVars.priVars()[eqIdx] - dirichletValues[eqIdx];
A[globalI][globalI][eqIdx][pvIdx] = (eqIdx == pvIdx) ? 1.0 : 0.0;
}
else
@@ -846,7 +846,7 @@ public:
{
if (internalDirichletConstraints[eqIdx])
{
- residual[eqIdx] = this->curElemVolVars()[scv].priVars()[eqIdx] - dirichletValues[eqIdx];
+ residual[eqIdx] = volVars.priVars()[eqIdx] - dirichletValues[eqIdx];
A[globalI][globalI][eqIdx][pvIdx] = (eqIdx == pvIdx) ? 1.0 : 0.0;
// inner faces
--
GitLab