CubicLawPermeability.h

Go to the documentation of this file.

 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/VariableType.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialPropertyLib
{
struct CubicLawPermeability final : public Property
{
    explicit CubicLawPermeability(std::string name,
                                  ParameterLib::Parameter<double> const& b)
        : _b(b)
    {
        name_ = std::move(name);
    }
 
    void checkScale() const override
    {
        if (!std::holds_alternative<Medium*>(scale_))
        {
            OGS_FATAL(
                "The property 'CubicLawPermeability' is implemented on the "
                "'media' scaleonly.");
        }
    }
 
    PropertyDataType value(VariableArray const& variable_array,
                           ParameterLib::SpatialPosition const& pos,
                           double const t, double const dt) const override;
    PropertyDataType dValue(VariableArray const& variable_array,
                            Variable const variable,
                            ParameterLib::SpatialPosition const& pos,
                            double const t, double const dt) const override;
 
private:
    ParameterLib::Parameter<double> const& _b;
};
}  // namespace MaterialPropertyLib