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

flow_ebos FIP: create the element iterator object as part of the for statement

Browse Source 
IMO, the previous solution was easier to read and a few lines shorter,
but I don't really care about how this is done...
This commit is contained in:
Andreas Lauser 2017-02-21 13:53:08 +01:00
parent 1f7be0d324
commit 9dab0323ce
1 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
opm/autodiff/BlackoilModelEbos.hpp
View File

@ -1013,9 +1013,11 @@ namespace Opm {
}
ElementContext elemCtx(ebosSimulator_);
auto elemIt = elemCtx.gridView().template begin</*codim=*/0>();
const auto& elemEndIt = elemCtx.gridView().template end</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
for (auto elemIt = elemCtx.gridView().template begin</*codim=*/0>();
elemIt != elemEndIt;
++elemIt)
{
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity) {
continue;
@ -1085,8 +1087,10 @@ namespace Opm {
std::vector<double> hcpv(dims, 0.0);
std::vector<double> pres(dims, 0.0);
elemIt = elemCtx.gridView().template begin</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
for (auto elemIt = elemCtx.gridView().template begin</*codim=*/0>();
elemIt != elemEndIt;
++elemIt)
{
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity) {
continue;
@ -1121,8 +1125,10 @@ namespace Opm {
comm.sum(hcpv.data(), hcpv.size());
comm.sum(pres.data(), pres.size());
elemIt = elemCtx.gridView().template begin</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
for (auto elemIt = elemCtx.gridView().template begin</*codim=*/0>();
elemIt != elemEndIt;
++elemIt)
{
const auto& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity) {
continue;