IntegrationPointWriter.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "IntegrationPointWriter.h"
 
#include <range/v3/view/join.hpp>
 
#include "MeshLib/Mesh.h"
#include "MeshLib/Utils/getOrCreateMeshProperty.h"

static MeshLib::IntegrationPointMetaDataSingleField addIntegrationPointData(
    MeshLib::Mesh& mesh, MeshLib::IntegrationPointWriter const& writer)
{
    auto const& ip_values = writer.values(/*t, x, dof_table*/);
    assert(ip_values.size() == mesh.getNumberOfElements());
 
    // create field data and fill it with nodal values, and an offsets cell
    // array indicating where the cell's integration point data starts.
    auto& field_data = *MeshLib::getOrCreateMeshProperty<double>(
        mesh, writer.name(), MeshLib::MeshItemType::IntegrationPoint,
        writer.numberOfComponents());
 
    field_data.assign(ip_values | ranges::views::join);
 
    return {writer.name(), writer.numberOfComponents(),
            writer.integrationOrder()};
}

static void addIntegrationPointMetaDataSingleField(
    MeshLib::Mesh& mesh, MeshLib::IntegrationPointMetaData const& ip_meta_data)
{
    // Store the field data.
    auto& dictionary = *MeshLib::getOrCreateMeshProperty<char>(
        mesh, "IntegrationPointMetaData",
        MeshLib::MeshItemType::IntegrationPoint, 1);
    dictionary.assign(ip_meta_data.toJsonString());
}
 
namespace MeshLib
{
void addIntegrationPointDataToMesh(
    MeshLib::Mesh& mesh,
    std::vector<std::unique_ptr<IntegrationPointWriter>> const&
        integration_point_writer)
{
    auto meta_data = IntegrationPointMetaData{
        integration_point_writer |
        ranges::views::transform(
            [&](auto const& ip_writer)
            { return addIntegrationPointData(mesh, *ip_writer); })};
 
    if (!meta_data.empty())
    {
        addIntegrationPointMetaDataSingleField(mesh, meta_data);
    }
}
 
std::optional<IntegrationPointMetaData> getIntegrationPointMetaData(
    MeshLib::Properties const& properties)
{
    if (!properties.existsPropertyVector<char>("IntegrationPointMetaData"))
    {
        return std::nullopt;
    }
    auto const& mesh_property_ip_meta_data =
        *properties.template getPropertyVector<char>(
            "IntegrationPointMetaData");
 
    if (mesh_property_ip_meta_data.getMeshItemType() !=
        MeshLib::MeshItemType::IntegrationPoint)
    {
        OGS_FATAL("IntegrationPointMetaData array must be field data.");
    }
 
    // Find the current integration point data entry and extract the
    // meta data.
    return IntegrationPointMetaData{std::string_view{
        mesh_property_ip_meta_data.data(), mesh_property_ip_meta_data.size()}};
}
}  // namespace MeshLib