Linear.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include "MaterialLib/MPL/Property.h"
#include "MaterialLib/MPL/VariableType.h"
 
namespace MaterialPropertyLib
{
using StringOrVariable = std::variant<std::string, Variable>;
struct IndependentVariable
{
    StringOrVariable type;
    VariableType reference_condition;  // scalar or vector
    VariableType slope;                // scalar or matrix
    std::optional<double> min;
    std::optional<double> max;
 
    double valueClamped(VariableArray const& variable_array,
                        ParameterLib::SpatialPosition const& pos,
                        double const t) const;
 
    double valueUnclamped(VariableArray const& variable_array,
                          ParameterLib::SpatialPosition const& pos,
                          double const t) const;
};

class Linear final : public Property
{
public:
    Linear(std::string name,
           PropertyDataType const& property_reference_value,
           std::vector<IndependentVariable> const& vs);

    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;

    PropertyDataType d2Value(VariableArray const& variable_array,
                             Variable const pv1, Variable const pv2,
                             ParameterLib::SpatialPosition const& /*pos*/,
                             double const /*t*/,
                             double const /*dt*/) const override;
 
private:
    std::vector<IndependentVariable> const independent_variables_;
};
}  // namespace MaterialPropertyLib