Skip to content
Snippets Groups Projects

Fix/cc internal dirichlet

Merged Kilian Weishaupt requested to merge fix/cc-internal-dirichlet into master
2 files
+ 7
7
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -284,7 +284,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
@@ -454,7 +454,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
@@ -579,7 +579,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
@@ -660,7 +660,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;
}
}
Loading