Commit 438cc76a authored by Dennis Gläser's avatar Dennis Gläser Committed by Timo Koch
Browse files

[mpfa][localassembler] detect zero transmissibilities also in outside scvs

parent b52068c7
......@@ -250,11 +250,13 @@ private:
{
if (abs(wijk[faceIdx][0][localDir]) < wijZeroThresh)
{
insideZeroWij = !curGlobalScvf.boundary();
if (!curIsDirichlet)
{
insideZeroWij = true;
faceMarkers.emplace_back( std::make_pair(curLocalDofIdx, faceIdx) );
}
}
}
// if we are not on a Dirichlet face, add entries associated with unknown face pressures
// i.e. in matrix C and maybe A (if current face is not a Dirichlet face)
......@@ -308,7 +310,8 @@ private:
const auto& otherLocalScvf = iv.localScvf(otherLocalScvfIdx);
const auto otherLocalDofIdx = otherLocalScvf.localDofIndex();
if (otherLocalDofIdx == curLocalDofIdx && insideZeroWij)
// check for zero transmissibilities (skip if inside has been zero already)
if (otherLocalDofIdx == curLocalDofIdx && !insideZeroWij)
if (abs(wijk[faceIdx][idxOnScvf][localDir]) < wijZeroThresh)
faceMarkers.emplace_back( std::make_pair(curLocalDofIdx, faceIdx) );
......
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