eclproblem: use elements range generator

This commit is contained in:
Arne Morten Kvarving 2022-10-12 14:27:20 +02:00
parent be84bc16be
commit 3599e0e5de

View File

@ -2467,13 +2467,7 @@ private:
auto& sol = this->model().solution(/*timeIdx=*/0); auto& sol = this->model().solution(/*timeIdx=*/0);
const auto& gridView = this->gridView(); const auto& gridView = this->gridView();
ElementContext elemCtx(simulator); ElementContext elemCtx(simulator);
auto elemIt = gridView.template begin</*codim=*/0>(); for (const auto& elem : elements(gridView, Dune::Partitions::interior)) {
const auto& elemEndIt = gridView.template end</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity)
continue;
elemCtx.updatePrimaryStencil(elem); elemCtx.updatePrimaryStencil(elem);
int elemIdx = elemCtx.globalSpaceIndex(/*spaceIdx=*/0, /*timeIdx=*/0); int elemIdx = elemCtx.globalSpaceIndex(/*spaceIdx=*/0, /*timeIdx=*/0);
initial(sol[elemIdx], elemCtx, /*spaceIdx=*/0, /*timeIdx=*/0); initial(sol[elemIdx], elemCtx, /*spaceIdx=*/0, /*timeIdx=*/0);