SelectSolidConstitutiveRelation.h

Go to the documentation of this file.

 
#pragma once
 
#include <spdlog/fmt/ranges.h>
 
#include <map>
#include <range/v3/view/map.hpp>
 
#include "MeshLib/PropertyVector.h"
 
namespace MaterialLib
{
namespace Solids
{
template <typename SolidMaterialsMap>
auto& selectSolidConstitutiveRelation(
    SolidMaterialsMap const& constitutive_relations,
    MeshLib::PropertyVector<int> const* const material_ids,
    std::size_t const element_id)
{
    // Multiple constitutive relations and no material ids should not be valid.
    if (constitutive_relations.size() > 1 && material_ids == nullptr)
    {
        OGS_FATAL(
            "There are {} constitutive relations provided in the project file "
            "but no MaterialIDs could be found in the mesh.",
            constitutive_relations.size());
    }
 
    int const material_id = ((constitutive_relations.size() == 1 &&
                              constitutive_relations.begin()->first == 0) ||
                             material_ids == nullptr)
                                ? 0
                                : (*material_ids)[element_id];
 
    auto const constitutive_relation = constitutive_relations.find(material_id);
    if (constitutive_relation == end(constitutive_relations))
    {
        OGS_FATAL(
            "No constitutive relation found for material id {:d} and element "
            "{:d}. There are {:d} constitutive relations available, "
            "corresponding to the ids: {}",
            material_id, element_id, constitutive_relations.size(),
            fmt::join(constitutive_relations | ranges::views::keys, " "));
    }
 
    if (constitutive_relation->second == nullptr)
    {
        OGS_FATAL(
            "The constitutive relation found for material id {:d} and element "
            "{:d} is a nullptr, which is impossible.",
            material_id, element_id);
    }
 
    return *constitutive_relation->second;
}
}  // namespace Solids
}  // namespace MaterialLib