OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

fix incrementer in boundarycontext

Browse Source
This commit is contained in:
Tor Harald Sandve 2021-01-13 13:01:41 +01:00
parent 732931d952
commit 1b4f0cf2c5
2 changed files with 12 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
opm/models/discretization/common/fvbaseboundarycontext.hh
View File

@ -77,8 +77,15 @@ public:
void increment() void increment()
{ {
const auto& iend = gridView().iend(element()); const auto& iend = gridView().iend(element());
while (intersectionIt_ != iend && !intersectionIt_->boundary())
++ intersectionIt_; if(intersectionIt_ == iend)
return;
++intersectionIt_;
// iterate to the next boundary intersection
while (intersectionIt_ != iend && !intersectionIt_->boundary()) {
++intersectionIt_;
}
} }
/*! /*!

3
opm/models/discretization/common/fvbaselocalresidual.hh
View File

@ -425,6 +425,9 @@ protected:
return; return;
BoundaryContext boundaryCtx(elemCtx); BoundaryContext boundaryCtx(elemCtx);
// move the iterator to the first boundary
if(!boundaryCtx.intersection(0).boundary())
boundaryCtx.increment();
// evaluate the boundary for all boundary faces of the current context // evaluate the boundary for all boundary faces of the current context
size_t numBoundaryFaces = boundaryCtx.numBoundaryFaces(/*timeIdx=*/0); size_t numBoundaryFaces = boundaryCtx.numBoundaryFaces(/*timeIdx=*/0);