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

changed: consistently use std::array

Browse Source
This commit is contained in:
Arne Morten Kvarving 2022-08-05 10:35:52 +02:00
parent 874e4f818c
commit 6ed8a38ba5
3 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
ebos/ecloutputblackoilmodule.hh
View File

@ -804,7 +804,7 @@ public:
if (!this->oilPressure_.empty()) { if (!this->oilPressure_.empty()) {
// this assumes that capillary pressures only depend on the phase saturations // this assumes that capillary pressures only depend on the phase saturations
// and possibly on temperature. (this is always the case for ECL problems.) // and possibly on temperature. (this is always the case for ECL problems.)
Dune::FieldVector< Scalar, numPhases > pc(0); std::array<Scalar, numPhases> pc = {0};
const MaterialLawParams& matParams = simulator_.problem().materialLawParams(elemIdx); const MaterialLawParams& matParams = simulator_.problem().materialLawParams(elemIdx);
MaterialLaw::capillaryPressures(pc, matParams, fs); MaterialLaw::capillaryPressures(pc, matParams, fs);
Valgrind::CheckDefined(this->oilPressure_[elemIdx]); Valgrind::CheckDefined(this->oilPressure_[elemIdx]);

2
ebos/eclproblem.hh
View File

@ -2625,7 +2625,7 @@ private:
// this assumes that capillary pressures only depend on the phase saturations // this assumes that capillary pressures only depend on the phase saturations
// and possibly on temperature. (this is always the case for ECL problems.) // and possibly on temperature. (this is always the case for ECL problems.)
Dune::FieldVector<Scalar, numPhases> pc(0.0); std::array<Scalar, numPhases> pc = {0};
const auto& matParams = materialLawParams(dofIdx); const auto& matParams = materialLawParams(dofIdx);
MaterialLaw::capillaryPressures(pc, matParams, dofFluidState); MaterialLaw::capillaryPressures(pc, matParams, dofFluidState);
Valgrind::CheckDefined(pressure); Valgrind::CheckDefined(pressure);

6
ebos/equil/equilibrationhelpers.cc
View File

@ -352,8 +352,7 @@ operator()(double s) const
fluidState.setSaturation(FluidSystem::gasPhaseIdx, 0.0); fluidState.setSaturation(FluidSystem::gasPhaseIdx, 0.0);
fluidState.setSaturation(phase_, s); fluidState.setSaturation(phase_, s);
double pc[FluidSystem::numPhases]; std::array<double, FluidSystem::numPhases> pc{0.0};
std::fill(pc, pc + FluidSystem::numPhases, 0.0);
using MaterialLaw = typename MaterialLawManager::MaterialLaw; using MaterialLaw = typename MaterialLawManager::MaterialLaw;
MaterialLaw::capillaryPressures(pc, matParams, fluidState); MaterialLaw::capillaryPressures(pc, matParams, fluidState);
double sign = (phase_ == FluidSystem::waterPhaseIdx)? -1.0 : 1.0; double sign = (phase_ == FluidSystem::waterPhaseIdx)? -1.0 : 1.0;
@ -388,8 +387,7 @@ operator()(double s) const
fluidState.setSaturation(phase1_, s); fluidState.setSaturation(phase1_, s);
fluidState.setSaturation(phase2_, 1.0 - s); fluidState.setSaturation(phase2_, 1.0 - s);
double pc[FluidSystem::numPhases]; std::array<double, FluidSystem::numPhases> pc {0.0};
std::fill(pc, pc + FluidSystem::numPhases, 0.0);
using MaterialLaw = typename MaterialLawManager::MaterialLaw; using MaterialLaw = typename MaterialLawManager::MaterialLaw;
MaterialLaw::capillaryPressures(pc, matParams, fluidState); MaterialLaw::capillaryPressures(pc, matParams, fluidState);