Skip to content
Snippets Groups Projects
Commit 2c2b5249 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[md][cclocalassembler] Use undeflected privars for setting internal Dirichlet

parent f1d6ce6c
No related branches found
No related tags found
1 merge request!2240Fix/cc internal dirichlet
Hide whitespace changes
Inline Side-by-side
dumux/multidomain/subdomaincclocalassembler.hh
+ 3
− 3
View file @ 2c2b5249
......@@ -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
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment