Commit 12498bb0 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

Merge branch 'cherry-pick-995ed673' into 'master'

[md][fclocalassembler] Fix Dirichlet conditions

See merge request !2991
parents e44273bf 2a82f953
Pipeline #12697 passed with stages
in 0 seconds
......@@ -388,12 +388,12 @@ public:
// For other dofs, add the contribution of the partial derivative.
for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
{
for (int pvIdx = 0; pvIdx < GridView::dimension; ++pvIdx)
for (int pvIdx = 0; pvIdx < JacobianBlock::block_type::cols; ++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;
if (bcTypes.isCouplingDirichlet(eqIdx))
A[scv.dofIndex()][globalJ][eqIdx][pvIdx] = partialDerivs[scv.localDofIndex()][pvIdx];
else if (bcTypes.isDirichlet(eqIdx))
A[scv.dofIndex()][globalJ][eqIdx][pvIdx] = 0.0;
}
}
}
......
Supports Markdown
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