From 81ac95418cdb084de471e57833994fc4e45621a0 Mon Sep 17 00:00:00 2001 From: Arne Morten Kvarving Date: Fri, 16 May 2025 09:10:52 +0200 Subject: [PATCH] richardsintensivequantities: cosmetics --- opm/models/richards/richardsintensivequantities.hh | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/opm/models/richards/richardsintensivequantities.hh b/opm/models/richards/richardsintensivequantities.hh index 9c054a7382..252011516c 100644 --- a/opm/models/richards/richardsintensivequantities.hh +++ b/opm/models/richards/richardsintensivequantities.hh @@ -72,14 +72,13 @@ class RichardsIntensiveQuantities using DimMatrix = Dune::FieldMatrix; using ScalarPhaseVector = Dune::FieldVector; using PhaseVector = Dune::FieldVector; - using Toolbox = Opm::MathToolbox; + using Toolbox = MathToolbox; public: //! The type returned by the fluidState() method - using FluidState = Opm::ImmiscibleFluidState; + using FluidState = ImmiscibleFluidState; - RichardsIntensiveQuantities() - {} + RichardsIntensiveQuantities() = default; RichardsIntensiveQuantities(const RichardsIntensiveQuantities& other) = default; @@ -115,7 +114,7 @@ public: // reference pressure if the medium is fully // saturated by the wetting phase const Evaluation& pW = priVars.makeEvaluation(pressureWIdx, timeIdx); - Evaluation pN = + const Evaluation pN = Toolbox::max(elemCtx.problem().referencePressure(elemCtx, dofIdx, /*timeIdx=*/0), pW + (pC[gasPhaseIdx] - pC[liquidPhaseIdx])); @@ -147,8 +146,9 @@ public: MaterialLaw::relativePermeabilities(relativePermeability_, materialParams, fluidState_); // mobilities - for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) - mobility_[phaseIdx] = relativePermeability_[phaseIdx]/fluidState_.viscosity(phaseIdx); + for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) { + mobility_[phaseIdx] = relativePermeability_[phaseIdx] / fluidState_.viscosity(phaseIdx); + } // porosity porosity_ = problem.porosity(elemCtx, dofIdx, timeIdx);