ParameterLib/Parameter.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "Parameter.h"
 
#include "BaseLib/ConfigTree.h"
#include "BaseLib/Error.h"
#include "ConstantParameter.h"
#include "CurveScaledParameter.h"
#include "FunctionParameter.h"
#include "GroupBasedParameter.h"
#include "MeshElementParameter.h"
#include "MeshNodeParameter.h"
#include "RandomFieldMeshElementParameter.h"
#include "RasterParameter.h"
#include "TimeDependentHeterogeneousParameter.h"
 
namespace ParameterLib
{
std::unique_ptr<ParameterBase> createParameter(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<MeshLib::Mesh>> const& meshes,
    std::vector<GeoLib::NamedRaster> const& named_rasters,
    std::map<std::string,
             std::unique_ptr<MathLib::PiecewiseLinearInterpolation>> const&
        curves)
{
    auto const name = config.getConfigParameter<std::string>("name");
    auto const type = config.peekConfigParameter<std::string>("type");
 
    // Either the mesh name is given, or the first mesh's name will be
    // taken.
    auto const mesh_name =
        config.getConfigParameter<std::string>("mesh", meshes[0]->getName());
 
    auto const& mesh = MeshLib::findMeshByName(meshes, mesh_name);
 
    // Create parameter based on the provided type.
    if (type == "Constant")
    {
        INFO("ConstantParameter: {:s}", name);
        return createConstantParameter(name, config);
    }
    if (type == "CurveScaled")
    {
        INFO("CurveScaledParameter: {:s}", name);
        return createCurveScaledParameter(name, config, curves);
    }
    if (type == "Function")
    {
        INFO("FunctionParameter: {:s}", name);
        return createFunctionParameter(name, config, curves);
    }
    if (type == "Group")
    {
        INFO("GroupBasedParameter: {:s}", name);
        return createGroupBasedParameter(name, config, mesh);
    }
    if (type == "MeshElement")
    {
        INFO("MeshElementParameter: {:s}", name);
        return createMeshElementParameter(name, config, mesh);
    }
    if (type == "MeshNode")
    {
        INFO("MeshNodeParameter: {:s}", name);
        return createMeshNodeParameter(name, config, mesh);
    }
    if (type == "RandomFieldMeshElement")
    {
        auto& mesh_var = MeshLib::findMeshByName(meshes, mesh_name);
        INFO("RandomFieldMeshElement: {:s}", name);
        return createRandomFieldMeshElementParameter(name, config, mesh_var);
    }
    if (type == "Raster")
    {
        INFO("RasterParameter: {:s}", name);
        return createRasterParameter(name, config, named_rasters);
    }
    if (type == "TimeDependentHeterogeneousParameter")
    {
        INFO("TimeDependentHeterogeneousParameter: {:s}", name);
        return createTimeDependentHeterogeneousParameter(name, config);
    }
 
    OGS_FATAL("Cannot construct a parameter of given type '{:s}'.", type);
}
 
std::optional<std::string> isDefinedOnSameMesh(ParameterBase const& parameter,
                                               MeshLib::Mesh const& mesh)
{
    // Arbitrary domain of definition.
    if (parameter.mesh() == nullptr)
    {
        return {};
    }
 
    // Equal meshes.
    if (*parameter.mesh() == mesh)
    {
        return {};
    }
 
    return "The parameter's domain of definition mesh '" +
           parameter.mesh()->getName() + "' differs from the used mesh '" +
           mesh.getName() +
           "'. The same mesh (the same name) has to be referenced in the "
           "project file. Possible reasons are:\n - the parameter used for the "
           "initial condition is not defined on the bulk mesh,\n - the "
           "parameter's domain of definition mesh differs from the boundary "
           "condition or source term domain of definition mesh.";
}
 
Parameter<double>& getNamedOrCreateInlineParameter(
    BaseLib::ConfigTree const& config,
    std::vector<std::unique_ptr<ParameterBase>>& parameters,
    std::string const& property_name,
    std::string const& tag_name,
    std::string const& inline_suffix)
{
    std::string const raw = config.getConfigParameter<std::string>(tag_name);
 
    // try to parse number(s) (if empty - no inline-values)
    std::size_t bad_idx = 0;
    if (auto values = BaseLib::tryParseVector<double>(raw, &bad_idx); values)
    {
        if (values->empty())
        {
            OGS_FATAL(
                "Empty inline value list in <{:s}> for property '{:s}'. "
                "Provide at least one numeric value (e.g., \"1.23 4.56\") or "
                "specify the name of an existing parameter. "
                "Raw input was: \"{:s}\".",
                tag_name, property_name, raw);
        }
 
        // collect all existing parameter names for checks against uniqueness
        std::vector<std::string> existing_names;
        existing_names.reserve(parameters.size());
        std::transform(parameters.begin(), parameters.end(),
                       std::back_inserter(existing_names),
                       [](auto const& p) { return p->name; });
 
        // assign inline values as constant parameter
        std::string pname = BaseLib::getUniqueName(
            existing_names, property_name + "_" + inline_suffix);
        using P = ConstantParameter<double>;
        parameters.push_back(std::make_unique<P>(pname, std::move(*values)));
        return static_cast<Parameter<double>&>(*parameters.back());
    }
 
    // No inline values -> 'raw' is just a parameter name, now find that
    // parameter.
    return findParameter<double>(raw, parameters, 0, nullptr);
}
}  // namespace ParameterLib