Commit 995ed673 authored by Kilian Weishaupt's avatar Kilian Weishaupt
Browse files

[md][localassembler] Fix Dirichlet conditions

parent 0e3c4c8c
Pipeline #12603 waiting for manual action with stages
......@@ -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