fixed: do not process aquifers attached to overlap cells

This commit is contained in:
Arne Morten Kvarving 2020-12-15 08:47:29 +01:00
parent a74e363c3c
commit 91336f30e7
2 changed files with 4 additions and 0 deletions

View File

@ -118,6 +118,8 @@ protected:
const auto& elemEndIt = gridView.template end</*codim=*/ 0>();
for (; elemIt != elemEndIt; ++elemIt) {
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity)
continue;
unsigned cell_index = elemMapper.index(elem);
int idx = this->cellToConnectionIdx_[cell_index];

View File

@ -120,6 +120,8 @@ protected:
const auto& elemEndIt = gridView.template end</*codim=*/ 0>();
for (; elemIt != elemEndIt; ++elemIt) {
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity)
continue;
unsigned cell_index = elemMapper.index(elem);
int idx = this->cellToConnectionIdx_[cell_index];