CreateFunction.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include <vector>
 
#include "BaseLib/ConfigTree.h"
#include "Function.h"
#include "MathLib/InterpolationAlgorithms/PiecewiseLinearInterpolation.h"
 
namespace MaterialPropertyLib
{
std::unique_ptr<Function> createFunction(
    BaseLib::ConfigTree const& config,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves)
{
    config.checkConfigParameter("type", "Function");
 
    // Second access for storage.
    auto property_name = config.peekConfigParameter<std::string>("name");
 
    DBUG("Create Function property {:s}.", property_name);
 
    std::vector<std::string> value_expressions;
    auto const& value_config = config.getConfigSubtree("value");

    for (auto const& p : value_config.getConfigSubtreeList("expression"))
    {
        value_expressions.emplace_back(p.getValue<std::string>());
    }
 
    // For each derivative a name of the variable and the list of expressions.
    std::vector<std::pair<std::string, std::vector<std::string>>>
        dvalue_expressions;
    for (auto const& dvalue_config : config.getConfigSubtreeList("dvalue"))
    {
        auto variable_name =
            dvalue_config.getConfigParameter<std::string>("variable_name");
 
        std::vector<std::string> expressions;
        auto const& expression_configs =
            dvalue_config.getConfigSubtreeList("expression");
 
        expressions.reserve(expression_configs.size());
        std::transform(std::begin(expression_configs),
                       std::end(expression_configs),
                       std::back_inserter(expressions),
                       [](BaseLib::ConfigTree const& p)
                       { return p.getValue<std::string>(); });
 
        dvalue_expressions.emplace_back(std::move(variable_name),
                                        std::move(expressions));
    }
 
    return std::make_unique<MaterialPropertyLib::Function>(
        std::move(property_name),
        value_expressions,
        dvalue_expressions,
        curves);
}
}  // namespace MaterialPropertyLib