Commit 7523aa37 authored by Martin Schneider's avatar Martin Schneider Committed by Timo Koch
Browse files

[wmpfa] Always make the flip index set

parent 1b03b307
......@@ -239,21 +239,18 @@ public:
scvfIndicesOfScv_[eIdx] = scvfsIndexSet;
}
// Make the flip index set for network, surface, and periodic grids
if (this->isPeriodic())
// Make the flip index set
flipScvfIndices_.resize(scvfs_.size());
for (auto&& scvf : scvfs_)
{
flipScvfIndices_.resize(scvfs_.size());
for (auto&& scvf : scvfs_)
{
if (scvf.boundary())
continue;
flipScvfIndices_[scvf.index()].resize(scvf.numOutsideScvs());
const auto insideScvIdx = scvf.insideScvIdx();
// check which outside scvf has the insideScvIdx index in its outsideScvIndices
for (unsigned int i = 0; i < scvf.numOutsideScvs(); ++i)
flipScvfIndices_[scvf.index()][i] = findFlippedScvfIndex_(insideScvIdx, scvf.outsideScvIdx(i));
}
if (scvf.boundary())
continue;
flipScvfIndices_[scvf.index()].resize(scvf.numOutsideScvs());
const auto insideScvIdx = scvf.insideScvIdx();
// check which outside scvf has the insideScvIdx index in its outsideScvIndices
for (unsigned int i = 0; i < scvf.numOutsideScvs(); ++i)
flipScvfIndices_[scvf.index()][i] = findFlippedScvfIndex_(insideScvIdx, scvf.outsideScvIdx(i));
}
// build the connectivity map for an effecient assembly
......
Markdown is supported
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