use compressed field properties setting up MaterialLawManager

This commit is contained in:
Arne Morten Kvarving 2020-01-28 11:16:31 +01:00
parent 36761e9347
commit 3c050f4dda
2 changed files with 2 additions and 7 deletions

View File

@ -2326,11 +2326,6 @@ private:
////////////////////////////////
// fluid-matrix interactions (saturation functions; relperm/capillary pressure)
size_t numDof = this->model().numGridDof();
std::vector<int> compressedToCartesianElemIdx(numDof);
for (unsigned elemIdx = 0; elemIdx < numDof; ++elemIdx)
compressedToCartesianElemIdx[elemIdx] = vanguard.cartesianIndex(elemIdx);
materialLawManager_ = std::make_shared<EclMaterialLawManager>();
if (comm.rank() == 0)
materialLawManager_->initFromDeck(deck, eclState);
@ -2338,7 +2333,7 @@ private:
EclMpiSerializer ser(comm);
ser.broadcast(*materialLawManager_);
materialLawManager_->initParamsForElements(eclState, compressedToCartesianElemIdx);
materialLawManager_->initParamsForElements(eclState, this->model().numGridDof());
////////////////////////////////
}

View File

@ -56,7 +56,7 @@ namespace Opm {
const auto dims = Opm::UgGridHelpers::cartDims(grid);
const auto& compressedToCartesianIdx = Opm::compressedToCartesian(nc, global_cell);
scaledEpsInfo_.resize(nc);
EclEpsGridProperties epsGridProperties(eclState, false, compressedToCartesianIdx);
EclEpsGridProperties epsGridProperties(eclState, false);
const std::string tag = "Scaled endpoints";
for (int c = 0; c < nc; ++c) {
const std::string satnumIdx = std::to_string(epsGridProperties.satRegion(c));