IntegrationPointWriter.cpp File Reference

Detailed Description

Copyright

Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org) Distributed under a Modified BSD License. See accompanying file LICENSE.txt or http://www.opengeosys.org/project/license

Definition in file IntegrationPointWriter.cpp.

#include "IntegrationPointWriter.h"
#include <nlohmann/json.hpp>
#include "MeshLib/Mesh.h"

Include dependency graph for IntegrationPointWriter.cpp:

Go to the source code of this file.

Namespaces
	ProcessLib

Functions
static ProcessLib::IntegrationPointMetaData	addIntegrationPointData (MeshLib::Mesh &mesh, ProcessLib::IntegrationPointWriter const &writer)
static void	addIntegrationPointMetaData (MeshLib::Mesh &mesh, std::vector< ProcessLib::IntegrationPointMetaData > const &meta_data)
static ProcessLib::IntegrationPointMetaData	extractIntegrationPointMetaData (json const &meta_data, std::string const &name)
void	ProcessLib::addIntegrationPointWriter (MeshLib::Mesh &mesh, std::vector< std::unique_ptr< IntegrationPointWriter >> const &integration_point_writer)
IntegrationPointMetaData	ProcessLib::getIntegrationPointMetaData (MeshLib::Mesh const &mesh, std::string const &name)

Function Documentation

addIntegrationPointData()

static ProcessLib::IntegrationPointMetaData addIntegrationPointData ( MeshLib::Mesh &  mesh, ProcessLib::IntegrationPointWriter const &  writer  )

static

Adds the integration point data and creates meta data for it.

Returns meta data for the written integration point data.

Definition at line 22 of file IntegrationPointWriter.cpp.

{
    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.clear();
 
    for (const auto& element_ip_values : ip_values)
    {
        std::copy(element_ip_values.begin(), element_ip_values.end(),
                  std::back_inserter(field_data));
    }
 
    return {writer.name(), writer.numberOfComponents(),
            writer.integrationOrder()};
}

References MathLib::LinAlg::copy(), MeshLib::Mesh::getNumberOfElements(), ProcessLib::IntegrationPointWriter::integrationOrder(), MeshLib::IntegrationPoint, ProcessLib::IntegrationPointWriter::name(), ProcessLib::IntegrationPointWriter::numberOfComponents(), and ProcessLib::IntegrationPointWriter::values().

Referenced by ProcessLib::addIntegrationPointWriter().

addIntegrationPointMetaData()

static void addIntegrationPointMetaData ( MeshLib::Mesh &  mesh, std::vector< ProcessLib::IntegrationPointMetaData > const &  meta_data  )

static

Adds integration point meta data as char mesh property encoded in JSON format, which is then stored as VTK's field data.

Definition at line 47 of file IntegrationPointWriter.cpp.

{
    json json_meta_data;
    json_meta_data["integration_point_arrays"] = json::array();
 
    for (auto const& md : meta_data)
    {
        json_meta_data["integration_point_arrays"].push_back(
            {{"name", md.name},
             {"number_of_components", md.n_components},
             {"integration_order", md.integration_order}});
    }
 
    // Store the field data.
    std::string const json_string = json_meta_data.dump();
    auto& dictionary = *MeshLib::getOrCreateMeshProperty<char>(
        mesh, "IntegrationPointMetaData",
        MeshLib::MeshItemType::IntegrationPoint, 1);
    dictionary.clear();
    std::copy(json_string.begin(), json_string.end(),
              std::back_inserter(dictionary));
}

References MathLib::LinAlg::copy(), and MeshLib::IntegrationPoint.

Referenced by ProcessLib::addIntegrationPointWriter().

extractIntegrationPointMetaData()

static ProcessLib::IntegrationPointMetaData extractIntegrationPointMetaData ( json const &  meta_data, std::string const &  name  )

static

For the given json object and the name extract integration point meta data, or fail if no meta data was found for the given name.

Definition at line 74 of file IntegrationPointWriter.cpp.

{
    auto const& ip_meta_data = meta_data["integration_point_arrays"];
    if (auto const it =
            find_if(cbegin(ip_meta_data), cend(ip_meta_data),
                    [&name](auto const& md) { return md["name"] == name; });
        it != cend(ip_meta_data))
    {
        return {name, (*it)["number_of_components"],
                (*it)["integration_order"]};
    }
    OGS_FATAL("No integration point meta data with name '{:s}' found.", name);
}

References MaterialPropertyLib::name, and OGS_FATAL.

Referenced by ProcessLib::getIntegrationPointMetaData().