From f08ee2048d16809f5742b6acb586f9fca8cce9b9 Mon Sep 17 00:00:00 2001 From: Thomas Fischer <thomas.fischer@ufz.de> Date: Fri, 16 Sep 2022 13:05:25 +0200 Subject: [PATCH] [MPL] Unify function parameter names in decl/def --- .../MPL/Properties/AverageMolarMass.cpp | 15 +++++----- MaterialLib/MPL/Properties/AverageMolarMass.h | 2 +- ...pillaryPressureRegularizedVanGenuchten.cpp | 4 +-- .../CapillaryPressureVanGenuchten.cpp | 4 +-- .../SaturationBrooksCorey.cpp | 15 +++++----- .../SaturationExponential.cpp | 11 ++++--- .../SaturationExponential.h | 6 ++-- .../SaturationLiakopoulos.cpp | 15 +++++----- .../SaturationVanGenuchten.cpp | 19 ++++++------ .../MPL/Properties/ClausiusClapeyron.cpp | 13 ++++----- .../MPL/Properties/ClausiusClapeyron.h | 6 ++-- MaterialLib/MPL/Properties/Curve.cpp | 4 +-- MaterialLib/MPL/Properties/Curve.h | 2 +- .../Properties/Density/WaterVapourDensity.cpp | 6 ++-- .../Properties/Density/WaterVapourDensity.h | 2 +- .../EmbeddedFracturePermeability.cpp | 4 +-- .../Properties/EmbeddedFracturePermeability.h | 2 +- .../Enthalpy/LinearWaterVapourLatentHeat.cpp | 4 +-- .../Enthalpy/LinearWaterVapourLatentHeat.h | 2 +- ...apourLatentHeatWithCriticalTemperature.cpp | 4 +-- ...rVapourLatentHeatWithCriticalTemperature.h | 2 +- MaterialLib/MPL/Properties/Exponential.cpp | 4 +-- MaterialLib/MPL/Properties/Exponential.h | 2 +- MaterialLib/MPL/Properties/Function.cpp | 8 ++--- MaterialLib/MPL/Properties/Function.h | 2 +- MaterialLib/MPL/Properties/IdealGasLaw.cpp | 29 +++++++++---------- .../Properties/IdealGasLawBinaryMixture.cpp | 7 ++--- MaterialLib/MPL/Properties/Linear.cpp | 11 ++++--- MaterialLib/MPL/Properties/Linear.h | 2 +- ...rthotropicEmbeddedFracturePermeability.cpp | 4 +-- .../OrthotropicEmbeddedFracturePermeability.h | 2 +- MaterialLib/MPL/Properties/Parameter.cpp | 2 +- MaterialLib/MPL/Properties/Parameter.h | 2 +- .../PermeabilityOrthotropicPowerLaw.cpp | 4 +-- .../Properties/PorosityFromMassBalance.cpp | 3 +- .../RelPermUdellNonwettingPhase.cpp | 4 +-- .../RelPermUdellNonwettingPhase.h | 2 +- .../RelPermVanGenuchten.cpp | 4 +-- .../SaturationDependentSwelling.cpp | 4 +-- .../TransportPorosityFromMassBalance.cpp | 3 +- 40 files changed, 115 insertions(+), 126 deletions(-) diff --git a/MaterialLib/MPL/Properties/AverageMolarMass.cpp b/MaterialLib/MPL/Properties/AverageMolarMass.cpp index 54923693c1f..d5f51a7362e 100644 --- a/MaterialLib/MPL/Properties/AverageMolarMass.cpp +++ b/MaterialLib/MPL/Properties/AverageMolarMass.cpp @@ -83,12 +83,12 @@ PropertyDataType AverageMolarMass::value( } PropertyDataType AverageMolarMass::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const { - if ((primary_variable != Variable::phase_pressure) && - (primary_variable != Variable::temperature)) + if ((variable != Variable::phase_pressure) && + (variable != Variable::temperature)) { OGS_FATAL( "AverageMolarMass::dValue is implemented for derivatives with " @@ -114,7 +114,7 @@ PropertyDataType AverageMolarMass::dValue( // for binary phases only. auto const dxnC = phase->property(PropertyType::mole_fraction) .template dValue<Eigen::Vector2d>( - variable_array, primary_variable, pos, t, dt)[0]; + variable_array, variable, pos, t, dt)[0]; auto const M_0 = phase->component(0) .property(PropertyType::molar_mass) @@ -128,10 +128,9 @@ PropertyDataType AverageMolarMass::dValue( } // namespace MaterialPropertyLib PropertyDataType AverageMolarMass::d2Value( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable1*/, Variable const /*primary_variable2*/, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& /*variable_array*/, Variable const /*variable1*/, + Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { OGS_FATAL("AverageMolarMass::d2Value is not yet implemented."); diff --git a/MaterialLib/MPL/Properties/AverageMolarMass.h b/MaterialLib/MPL/Properties/AverageMolarMass.h index 72b2864c83f..331d011174e 100644 --- a/MaterialLib/MPL/Properties/AverageMolarMass.h +++ b/MaterialLib/MPL/Properties/AverageMolarMass.h @@ -38,7 +38,7 @@ public: double const /*t*/, double const /*dt*/) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureRegularizedVanGenuchten.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureRegularizedVanGenuchten.cpp index 02069f2a1b1..88b1b66cddc 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureRegularizedVanGenuchten.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureRegularizedVanGenuchten.cpp @@ -68,11 +68,11 @@ PropertyDataType CapillaryPressureRegularizedVanGenuchten::value( } PropertyDataType CapillaryPressureRegularizedVanGenuchten::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::liquid_saturation) + if (variable != Variable::liquid_saturation) { OGS_FATAL( "CapillaryPressureRegularizedVanGenuchten::dValue is implemented " diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureVanGenuchten.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureVanGenuchten.cpp index c44eb6e4d0c..44487d00674 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureVanGenuchten.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/CapillaryPressureVanGenuchten.cpp @@ -107,11 +107,11 @@ PropertyDataType CapillaryPressureVanGenuchten::value( } PropertyDataType CapillaryPressureVanGenuchten::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::liquid_saturation) + if (variable != Variable::liquid_saturation) { OGS_FATAL( "CapillaryPressureVanGenuchten::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationBrooksCorey.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationBrooksCorey.cpp index 36836d20b39..9780668df52 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationBrooksCorey.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationBrooksCorey.cpp @@ -56,11 +56,11 @@ PropertyDataType SaturationBrooksCorey::value( } PropertyDataType SaturationBrooksCorey::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::capillary_pressure) + if (variable != Variable::capillary_pressure) { OGS_FATAL( "SaturationBrooksCorey::dValue is implemented for derivatives with " @@ -85,13 +85,12 @@ PropertyDataType SaturationBrooksCorey::dValue( } PropertyDataType SaturationBrooksCorey::d2Value( - VariableArray const& variable_array, Variable const primary_variable1, - Variable const primary_variable2, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& variable_array, Variable const variable1, + Variable const variable2, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { - if ((primary_variable1 != Variable::capillary_pressure) && - (primary_variable2 != Variable::capillary_pressure)) + if ((variable1 != Variable::capillary_pressure) && + (variable2 != Variable::capillary_pressure)) { OGS_FATAL( "SaturationBrooksCorey::d2Value is implemented for derivatives " diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.cpp index af2a0869b29..de7c3aa50b5 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.cpp @@ -58,11 +58,11 @@ PropertyDataType SaturationExponential::value( } PropertyDataType SaturationExponential::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::capillary_pressure) + if (variable != Variable::capillary_pressure) { OGS_FATAL( "SaturationExponential::dValue is implemented for derivatives with " @@ -83,10 +83,9 @@ PropertyDataType SaturationExponential::dValue( } PropertyDataType SaturationExponential::d2Value( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable1*/, Variable const /*primary_variable2*/, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& /*variable_array*/, Variable const /*variable1*/, + Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { OGS_FATAL("SaturationExponential::d2Value() is not implemented."); } diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.h b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.h index f60d1624b91..9cdde6acff5 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.h +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationExponential.h @@ -68,13 +68,13 @@ public: double const /*t*/, double const /*dt*/) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; PropertyDataType d2Value(VariableArray const& variable_array, - Variable const /*primary_variable1*/, - Variable const /*primary_variable2*/, + Variable const /*variable1*/, + Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationLiakopoulos.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationLiakopoulos.cpp index ae2017caa62..7536285d87d 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationLiakopoulos.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationLiakopoulos.cpp @@ -37,11 +37,11 @@ PropertyDataType SaturationLiakopoulos::value( } PropertyDataType SaturationLiakopoulos::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::capillary_pressure) + if (variable != Variable::capillary_pressure) { OGS_FATAL( "SaturationLiakopoulos::dValue is implemented for derivatives with " @@ -60,13 +60,12 @@ PropertyDataType SaturationLiakopoulos::dValue( } PropertyDataType SaturationLiakopoulos::d2Value( - VariableArray const& variable_array, Variable const primary_variable1, - Variable const primary_variable2, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& variable_array, Variable const variable1, + Variable const variable2, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { - if ((primary_variable1 != Variable::capillary_pressure) && - (primary_variable2 != Variable::capillary_pressure)) + if ((variable1 != Variable::capillary_pressure) && + (variable2 != Variable::capillary_pressure)) { OGS_FATAL( "SaturationLiakopoulos::d2Value is implemented for derivatives " diff --git a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationVanGenuchten.cpp b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationVanGenuchten.cpp index 472d3be853c..2f53c7f119b 100644 --- a/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationVanGenuchten.cpp +++ b/MaterialLib/MPL/Properties/CapillaryPressureSaturation/SaturationVanGenuchten.cpp @@ -61,11 +61,11 @@ PropertyDataType SaturationVanGenuchten::value( } PropertyDataType SaturationVanGenuchten::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::capillary_pressure) + if (variable != Variable::capillary_pressure) { OGS_FATAL( "SaturationVanGenuchten::dValue is implemented for derivatives " @@ -98,15 +98,14 @@ PropertyDataType SaturationVanGenuchten::dValue( } PropertyDataType SaturationVanGenuchten::d2Value( - VariableArray const& variable_array, Variable const primary_variable1, - Variable const primary_variable2, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& variable_array, Variable const variable1, + Variable const variable2, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { - (void)primary_variable1; - (void)primary_variable2; - assert((primary_variable1 == Variable::capillary_pressure) && - (primary_variable2 == Variable::capillary_pressure) && + (void)variable1; + (void)variable2; + assert((variable1 == Variable::capillary_pressure) && + (variable2 == Variable::capillary_pressure) && "SaturationVanGenuchten::d2Value is implemented for derivatives " "with respect to capillary pressure only."); diff --git a/MaterialLib/MPL/Properties/ClausiusClapeyron.cpp b/MaterialLib/MPL/Properties/ClausiusClapeyron.cpp index c844d61ca68..10172c0a627 100644 --- a/MaterialLib/MPL/Properties/ClausiusClapeyron.cpp +++ b/MaterialLib/MPL/Properties/ClausiusClapeyron.cpp @@ -88,7 +88,7 @@ PropertyDataType ClausiusClapeyron::value( } PropertyDataType ClausiusClapeyron::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const { @@ -103,7 +103,7 @@ PropertyDataType ClausiusClapeyron::dValue( { return 0.; } - if (primary_variable == Variable::phase_pressure) + if (variable == Variable::phase_pressure) { return 0.; } @@ -112,7 +112,7 @@ PropertyDataType ClausiusClapeyron::dValue( const double dh = variable_array.enthalpy_of_evaporation; const double p_vap = std::get<double>(value(variable_array, pos, t, dt)); - if (primary_variable == Variable::temperature) + if (variable == Variable::temperature) { return p_vap * M * dh / (R * T * T); } @@ -122,10 +122,9 @@ PropertyDataType ClausiusClapeyron::dValue( } PropertyDataType ClausiusClapeyron::d2Value( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable1*/, Variable const /*primary_variable2*/, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& /*variable_array*/, Variable const /*variable1*/, + Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { OGS_FATAL("ClausiusClapeyron::d2Value is not implemented."); } diff --git a/MaterialLib/MPL/Properties/ClausiusClapeyron.h b/MaterialLib/MPL/Properties/ClausiusClapeyron.h index 455fcc6e50c..4c02b9a94d1 100644 --- a/MaterialLib/MPL/Properties/ClausiusClapeyron.h +++ b/MaterialLib/MPL/Properties/ClausiusClapeyron.h @@ -44,13 +44,13 @@ public: double const t, double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; PropertyDataType d2Value(VariableArray const& variable_array, - Variable const primary_variable1, - Variable const primary_variable2, + Variable const variable1, + Variable const variable2, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; diff --git a/MaterialLib/MPL/Properties/Curve.cpp b/MaterialLib/MPL/Properties/Curve.cpp index eb7f75a0f80..bc291a2ce0f 100644 --- a/MaterialLib/MPL/Properties/Curve.cpp +++ b/MaterialLib/MPL/Properties/Curve.cpp @@ -29,11 +29,11 @@ PropertyDataType Curve::value(VariableArray const& variable_array, } PropertyDataType Curve::dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != independent_variable_) + if (variable != independent_variable_) { return 0.0; } diff --git a/MaterialLib/MPL/Properties/Curve.h b/MaterialLib/MPL/Properties/Curve.h index e29f452c65e..f212e5e2e4d 100644 --- a/MaterialLib/MPL/Properties/Curve.h +++ b/MaterialLib/MPL/Properties/Curve.h @@ -37,7 +37,7 @@ public: /// This method will compute the derivative of a property with respect to /// the given primary variable. PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/Density/WaterVapourDensity.cpp b/MaterialLib/MPL/Properties/Density/WaterVapourDensity.cpp index a294d8a6051..f68b6142834 100644 --- a/MaterialLib/MPL/Properties/Density/WaterVapourDensity.cpp +++ b/MaterialLib/MPL/Properties/Density/WaterVapourDensity.cpp @@ -52,7 +52,7 @@ PropertyDataType WaterVapourDensity::value( } PropertyDataType WaterVapourDensity::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { @@ -60,7 +60,7 @@ PropertyDataType WaterVapourDensity::dValue( double const T = variable_array.temperature; double const water_density = variable_array.density; - if (primary_variable == Variable::temperature) + if (variable == Variable::temperature) { double const h = humidity(T, p, water_density); double const rho_v = h * saturatedVaporDensity(T); @@ -72,7 +72,7 @@ PropertyDataType WaterVapourDensity::dValue( SpecificGasConstant::WaterVapour); } - if (primary_variable == Variable::phase_pressure) + if (variable == Variable::phase_pressure) { double const factor = 1.0 / diff --git a/MaterialLib/MPL/Properties/Density/WaterVapourDensity.h b/MaterialLib/MPL/Properties/Density/WaterVapourDensity.h index 098b0fc510a..967342f481c 100644 --- a/MaterialLib/MPL/Properties/Density/WaterVapourDensity.h +++ b/MaterialLib/MPL/Properties/Density/WaterVapourDensity.h @@ -71,7 +71,7 @@ public: /// \return The derivative of water vapour density with respect to /// temperature or phase (water) pressure. PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp index 6998b4e8e7a..45e2f537b20 100644 --- a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp +++ b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp @@ -92,11 +92,11 @@ PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::value( template <int DisplacementDim> PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const /*dt*/) const { - if (primary_variable != Variable::mechanical_strain) + if (variable != Variable::mechanical_strain) { OGS_FATAL( "EmbeddedFracturePermeability::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.h b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.h index 8f9a10ded24..6c41e86f2ec 100644 --- a/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.h +++ b/MaterialLib/MPL/Properties/EmbeddedFracturePermeability.h @@ -84,7 +84,7 @@ public: ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.cpp b/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.cpp index 5afd5d2ee55..146217ca77b 100644 --- a/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.cpp +++ b/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.cpp @@ -31,11 +31,11 @@ PropertyDataType LinearWaterVapourLatentHeat::value( } PropertyDataType LinearWaterVapourLatentHeat::dValue( - VariableArray const& /*variable_array*/, Variable const primary_variable, + VariableArray const& /*variable_array*/, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable == Variable::temperature) + if (variable == Variable::temperature) { return -2369.2; } diff --git a/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.h b/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.h index d685fe1c611..975cd0591e6 100644 --- a/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.h +++ b/MaterialLib/MPL/Properties/Enthalpy/LinearWaterVapourLatentHeat.h @@ -57,7 +57,7 @@ public: double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.cpp b/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.cpp index 81a1b654640..449796bf91b 100644 --- a/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.cpp +++ b/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.cpp @@ -71,11 +71,11 @@ PropertyDataType WaterVapourLatentHeatWithCriticalTemperature::value( } PropertyDataType WaterVapourLatentHeatWithCriticalTemperature::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::temperature) + if (variable != Variable::temperature) { OGS_FATAL( "WaterVapourLatentHeatWithCriticalTemperature::dValue is " diff --git a/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.h b/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.h index 5c0d02bb6e9..87a8f6c2ce3 100644 --- a/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.h +++ b/MaterialLib/MPL/Properties/Enthalpy/WaterVapourLatentHeatWithCriticalTemperature.h @@ -82,7 +82,7 @@ public: double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/Exponential.cpp b/MaterialLib/MPL/Properties/Exponential.cpp index 8e924c8ae5a..40ddf9e4d89 100644 --- a/MaterialLib/MPL/Properties/Exponential.cpp +++ b/MaterialLib/MPL/Properties/Exponential.cpp @@ -39,11 +39,11 @@ PropertyDataType Exponential::value( } PropertyDataType Exponential::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (exponent_data_.type != primary_variable) + if (exponent_data_.type != variable) { return 0.; } diff --git a/MaterialLib/MPL/Properties/Exponential.h b/MaterialLib/MPL/Properties/Exponential.h index 252bef62e55..65757129e29 100644 --- a/MaterialLib/MPL/Properties/Exponential.h +++ b/MaterialLib/MPL/Properties/Exponential.h @@ -55,7 +55,7 @@ public: /// This method will compute the derivative of a property with respect /// to the given primary variable. PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/Function.cpp b/MaterialLib/MPL/Properties/Function.cpp index f4faaa95dcb..97e0ae21a85 100644 --- a/MaterialLib/MPL/Properties/Function.cpp +++ b/MaterialLib/MPL/Properties/Function.cpp @@ -191,21 +191,21 @@ PropertyDataType Function::value(VariableArray const& variable_array, } PropertyDataType Function::dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { auto const it = std::find_if(begin(dvalue_expressions_), end(dvalue_expressions_), - [&primary_variable](auto const& v) - { return v.first == primary_variable; }); + [&variable](auto const& v) + { return v.first == variable; }); if (it == end(dvalue_expressions_)) { OGS_FATAL( "Requested derivative with respect to the variable {:s} not " "provided for Function-type property {:s}.", - variable_enum_to_string[static_cast<int>(primary_variable)], name_); + variable_enum_to_string[static_cast<int>(variable)], name_); } return evaluateExpressions(symbol_values_, variable_array, it->second); diff --git a/MaterialLib/MPL/Properties/Function.h b/MaterialLib/MPL/Properties/Function.h index 4c531baeb7d..bcd0f0d5c1b 100644 --- a/MaterialLib/MPL/Properties/Function.h +++ b/MaterialLib/MPL/Properties/Function.h @@ -38,7 +38,7 @@ public: double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; diff --git a/MaterialLib/MPL/Properties/IdealGasLaw.cpp b/MaterialLib/MPL/Properties/IdealGasLaw.cpp index bf824a08070..d29340d27cd 100644 --- a/MaterialLib/MPL/Properties/IdealGasLaw.cpp +++ b/MaterialLib/MPL/Properties/IdealGasLaw.cpp @@ -41,7 +41,7 @@ PropertyDataType IdealGasLaw::value( } PropertyDataType IdealGasLaw::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { @@ -51,14 +51,14 @@ PropertyDataType IdealGasLaw::dValue( const double molar_mass = variable_array.molar_mass; // todo: add molar mass derivatives - if (primary_variable == Variable::temperature) + if (variable == Variable::temperature) { // extend to take temperature-dependent molar mass into account return -pressure * molar_mass / gas_constant / temperature / temperature; } - if (primary_variable == Variable::phase_pressure) + if (variable == Variable::phase_pressure) { // extend to take pressure-dependent molar mass into account return molar_mass / gas_constant / temperature; @@ -72,10 +72,9 @@ PropertyDataType IdealGasLaw::dValue( } PropertyDataType IdealGasLaw::d2Value( - VariableArray const& variable_array, Variable const primary_variable1, - Variable const primary_variable2, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& variable_array, Variable const variable1, + Variable const variable2, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { const double gas_constant = MaterialLib::PhysicalConstant::IdealGasConstant; const double pressure = variable_array.phase_pressure; @@ -83,25 +82,25 @@ PropertyDataType IdealGasLaw::d2Value( const double molar_mass = variable_array.molar_mass; // todo: add molar mass derivatives - if ((primary_variable1 == Variable::phase_pressure) && - (primary_variable2 == Variable::phase_pressure)) + if ((variable1 == Variable::phase_pressure) && + (variable2 == Variable::phase_pressure)) { // d2rho_dp2 // extend to take pressure-dependent molar mass into account return 0.; } - if ((primary_variable1 == Variable::temperature) && - (primary_variable2 == Variable::temperature)) + if ((variable1 == Variable::temperature) && + (variable2 == Variable::temperature)) { // d2rho_dT2 // extend to take temperature-dependent molar mass into account return 2. * molar_mass * pressure / gas_constant / temperature / temperature / temperature; } - if (((primary_variable1 == Variable::phase_pressure) && - (primary_variable2 == Variable::temperature)) || - ((primary_variable1 == Variable::temperature) && - (primary_variable2 == Variable::phase_pressure))) + if (((variable1 == Variable::phase_pressure) && + (variable2 == Variable::temperature)) || + ((variable1 == Variable::temperature) && + (variable2 == Variable::phase_pressure))) { // d2rho_dpdT or d2rho_dTdp // extend to take pressure-temperature-dependent molar mass into account diff --git a/MaterialLib/MPL/Properties/IdealGasLawBinaryMixture.cpp b/MaterialLib/MPL/Properties/IdealGasLawBinaryMixture.cpp index 2894aa48312..81d2ef38d64 100644 --- a/MaterialLib/MPL/Properties/IdealGasLawBinaryMixture.cpp +++ b/MaterialLib/MPL/Properties/IdealGasLawBinaryMixture.cpp @@ -67,10 +67,9 @@ PropertyDataType IdealGasLawBinaryMixture::dValue( } PropertyDataType IdealGasLawBinaryMixture::d2Value( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable1*/, Variable const /*primary_variable2*/, - ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, - double const /*dt*/) const + VariableArray const& /*variable_array*/, Variable const /*variable1*/, + Variable const /*variable2*/, ParameterLib::SpatialPosition const& /*pos*/, + double const /*t*/, double const /*dt*/) const { OGS_FATAL("IdealGasLawBinaryMixture::d2Value is not implemented."); diff --git a/MaterialLib/MPL/Properties/Linear.cpp b/MaterialLib/MPL/Properties/Linear.cpp index df3d87d09d5..6a39545a37c 100644 --- a/MaterialLib/MPL/Properties/Linear.cpp +++ b/MaterialLib/MPL/Properties/Linear.cpp @@ -46,15 +46,14 @@ PropertyDataType Linear::value(VariableArray const& variable_array, } PropertyDataType Linear::dValue(VariableArray const& /*variable_array*/, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - auto const independent_variable = - std::find_if(independent_variables_.begin(), - independent_variables_.end(), - [&primary_variable](auto const& iv) -> bool - { return iv.type == primary_variable; }); + auto const independent_variable = std::find_if( + independent_variables_.begin(), + independent_variables_.end(), + [&variable](auto const& iv) -> bool { return iv.type == variable; }); return independent_variable != independent_variables_.end() ? std::get<double>(value_) * diff --git a/MaterialLib/MPL/Properties/Linear.h b/MaterialLib/MPL/Properties/Linear.h index 484a9d1d45f..44765b88a68 100644 --- a/MaterialLib/MPL/Properties/Linear.h +++ b/MaterialLib/MPL/Properties/Linear.h @@ -44,7 +44,7 @@ public: /// This method will compute the derivative of a property with respect to /// the given primary variable. PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp index 29c46798ff5..80310580364 100644 --- a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp +++ b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.cpp @@ -87,11 +87,11 @@ OrthotropicEmbeddedFracturePermeability<DisplacementDim>::value( template <int DisplacementDim> PropertyDataType OrthotropicEmbeddedFracturePermeability<DisplacementDim>::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const /*dt*/) const { - if (primary_variable != Variable::mechanical_strain) + if (variable != Variable::mechanical_strain) { OGS_FATAL( "OrthotropicEmbeddedFracturePermeability::dValue is implemented " diff --git a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.h b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.h index 34ffadedf8c..16fae69ddc8 100644 --- a/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.h +++ b/MaterialLib/MPL/Properties/OrthotropicEmbeddedFracturePermeability.h @@ -87,7 +87,7 @@ public: ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/Parameter.cpp b/MaterialLib/MPL/Properties/Parameter.cpp index c097d9d197a..85e65e178ca 100644 --- a/MaterialLib/MPL/Properties/Parameter.cpp +++ b/MaterialLib/MPL/Properties/Parameter.cpp @@ -35,7 +35,7 @@ PropertyDataType Parameter::value(VariableArray const& /*variable_array*/, } PropertyDataType Parameter::dValue(VariableArray const& /*variable_array*/, - Variable const /*primary_variable*/, + Variable const /*variable*/, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const diff --git a/MaterialLib/MPL/Properties/Parameter.h b/MaterialLib/MPL/Properties/Parameter.h index 077282ff590..f1adaf5be7b 100644 --- a/MaterialLib/MPL/Properties/Parameter.h +++ b/MaterialLib/MPL/Properties/Parameter.h @@ -39,7 +39,7 @@ public: /// This method will compute the derivative of a property with respect to /// the given primary variable. PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const override; diff --git a/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp b/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp index 5bfc924f87c..9093af7c66a 100644 --- a/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp +++ b/MaterialLib/MPL/Properties/PermeabilityOrthotropicPowerLaw.cpp @@ -86,11 +86,11 @@ PropertyDataType PermeabilityOrthotropicPowerLaw<DisplacementDim>::value( } template <int DisplacementDim> PropertyDataType PermeabilityOrthotropicPowerLaw<DisplacementDim>::dValue( - VariableArray const& /*variable_array*/, Variable const primary_variable, + VariableArray const& /*variable_array*/, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::total_strain) + if (variable != Variable::total_strain) { OGS_FATAL( "PermeabilityOrthotropicPowerLaw::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp b/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp index 3c0196874d5..97a139d53ce 100644 --- a/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp +++ b/MaterialLib/MPL/Properties/PorosityFromMassBalance.cpp @@ -63,8 +63,7 @@ PropertyDataType PorosityFromMassBalance::value( } PropertyDataType PorosityFromMassBalance::dValue( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable*/, + VariableArray const& /*variable_array*/, Variable const /*variable*/, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.cpp b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.cpp index 4a5f231b4e6..760b8d5aa68 100644 --- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.cpp +++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.cpp @@ -67,11 +67,11 @@ PropertyDataType RelPermUdellNonwettingPhase::value( return std::max(min_relative_permeability_, S_e_g * S_e_g * S_e_g); } PropertyDataType RelPermUdellNonwettingPhase::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::liquid_saturation) + if (variable != Variable::liquid_saturation) { OGS_FATAL( "RelPermUdellNonwettingPhase::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.h b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.h index f7df53cef14..f49ec4c2e6a 100644 --- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.h +++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermUdellNonwettingPhase.h @@ -66,7 +66,7 @@ public: ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; PropertyDataType dValue(VariableArray const& variable_array, - Variable const primary_variable, + Variable const variable, ParameterLib::SpatialPosition const& pos, double const t, double const dt) const override; }; diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermVanGenuchten.cpp b/MaterialLib/MPL/Properties/RelativePermeability/RelPermVanGenuchten.cpp index eaaf1e248eb..beee8a916fe 100644 --- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermVanGenuchten.cpp +++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermVanGenuchten.cpp @@ -54,11 +54,11 @@ PropertyDataType RelPermVanGenuchten::value( } PropertyDataType RelPermVanGenuchten::dValue( - VariableArray const& variable_array, Variable const primary_variable, + VariableArray const& variable_array, Variable const variable, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::liquid_saturation) + if (variable != Variable::liquid_saturation) { OGS_FATAL( "RelativePermeabilityVanGenuchten::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp b/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp index fbeb1d60a42..68d8db3d8ff 100644 --- a/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp +++ b/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp @@ -103,11 +103,11 @@ PropertyDataType SaturationDependentSwelling::value( PropertyDataType SaturationDependentSwelling::dValue( VariableArray const& variable_array, - VariableArray const& variable_array_prev, Variable const primary_variable, + VariableArray const& variable_array_prev, Variable const variable, ParameterLib::SpatialPosition const& pos, double const /*t*/, double const /*dt*/) const { - if (primary_variable != Variable::liquid_saturation) + if (variable != Variable::liquid_saturation) { OGS_FATAL( "SaturationDependentSwelling::dValue is implemented for " diff --git a/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp b/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp index 93938219847..8a257335386 100644 --- a/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp +++ b/MaterialLib/MPL/Properties/TransportPorosityFromMassBalance.cpp @@ -65,8 +65,7 @@ PropertyDataType TransportPorosityFromMassBalance::value( } PropertyDataType TransportPorosityFromMassBalance::dValue( - VariableArray const& /*variable_array*/, - Variable const /*primary_variable*/, + VariableArray const& /*variable_array*/, Variable const /*variable*/, ParameterLib::SpatialPosition const& /*pos*/, double const /*t*/, double const /*dt*/) const { -- GitLab