CreateConstitutiveRelationsGeneric.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 "BaseLib/ConfigTree.h"
#include "MaterialLib/Utils/MediaCreation.h"
#include "MeshLib/PropertyVector.h"
#include "ParameterLib/Parameter.h"
 
namespace MaterialLib::Solids
{
template <typename SolidConstitutiveRelation>
std::map<int, std::shared_ptr<SolidConstitutiveRelation>>
createConstitutiveRelationsGeneric(
    std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const& parameters,
    std::optional<ParameterLib::CoordinateSystem> const&
        local_coordinate_system,
    MeshLib::PropertyVector<int> const* const material_ids,
    BaseLib::ConfigTree const& config,
    std::unique_ptr<SolidConstitutiveRelation> (*create_constitutive_relation)(
        std::vector<std::unique_ptr<ParameterLib::ParameterBase>> const&,
        std::optional<ParameterLib::CoordinateSystem> const&,
        BaseLib::ConfigTree const&))
{
    auto const constitutive_relation_configs =
        config.getConfigSubtreeList("constitutive_relation");
 
    std::map<int, std::shared_ptr<SolidConstitutiveRelation>>
        constitutive_relations;
 
    for (auto const& constitutive_relation_config :
         constitutive_relation_configs)
    {
        auto create = [&create_constitutive_relation,
                       &parameters,
                       &local_coordinate_system,
                       &constitutive_relation_config](int const /*id*/)
        {
            return create_constitutive_relation(parameters,
                                                local_coordinate_system,
                                                constitutive_relation_config);
        };
 
        auto const material_id_string =
            constitutive_relation_config.getConfigAttribute<std::string>("id",
                                                                         "0");
 
        auto const material_ids_of_this_constitutive_relation =
            MaterialLib::parseMaterialIdString(material_id_string,
                                               material_ids);
 
        for (auto const& material_id :
             material_ids_of_this_constitutive_relation)
        {
            MaterialLib::createMediumForId(
                material_id,
                constitutive_relations,
                material_ids_of_this_constitutive_relation,
                create);
        }
    }
 
    DBUG("Found {:d} constitutive relations.", constitutive_relations.size());
 
    return constitutive_relations;
}
 
}  // namespace MaterialLib::Solids