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

[md][fcsubdomainlocalassember] Revise incorporation of Dirichlet BCs

(cherry picked from commit 3b56a4b8)
parent 71ce45bc
No related branches found
No related tags found
1 merge request!2982[md][fcsubdomainlocalassember] Revise incorporation of Dirichlet BCs
Pipeline #12065 failed
+1
......@@ -388,10 +388,13 @@ public:
// For other dofs, add the contribution of the partial derivative.
for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
{
if (bcTypes.isCouplingDirichlet(eqIdx))
A[scv.dofIndex()][globalJ][eqIdx][pvIdx] = partialDerivs[scv.localDofIndex()][eqIdx];
else if (bcTypes.isDirichlet(eqIdx))
A[scv.dofIndex()][globalJ][eqIdx][pvIdx] = 0.0;
for (int pvIdx = 0; pvIdx < GridView::dimension; ++pvIdx)
{
if (bcTypes.isCouplingDirichlet(pvIdx) && pvIdx == scv.dofAxis()) // TODO?
A[scv.dofIndex()][globalJ][eqIdx][eqIdx] = partialDerivs[scv.localDofIndex()][eqIdx];
else if (bcTypes.isDirichlet(pvIdx) && pvIdx == scv.dofAxis())
A[scv.dofIndex()][globalJ][eqIdx][eqIdx] = 0.0;
}
}
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment