Detailed Description

Definition at line 18 of file CellAverageData.h.

#include <CellAverageData.h>

Collaboration diagram for ProcessLib::CellAverageData:

Public Member Functions
	CellAverageData (MeshLib::Mesh &mesh)

MeshLib::PropertyVector< double > &	getOrCreatePropertyVector (std::string const &name, unsigned const num_comp)

Private Attributes
MeshLib::Mesh const &	mesh_

std::map< std::string, MeshLib::PropertyVector< double > * >	cell_averages_

Constructor & Destructor Documentation

◆ CellAverageData()

ProcessLib::CellAverageData::CellAverageData ( MeshLib::Mesh & mesh )

inlineexplicit

Definition at line 20 of file CellAverageData.h.

20: mesh_{mesh} {}

ProcessLib::CellAverageData::mesh_

MeshLib::Mesh const & mesh_

Definition CellAverageData.h:26

Member Function Documentation

◆ getOrCreatePropertyVector()

MeshLib::PropertyVector< double > & ProcessLib::CellAverageData::getOrCreatePropertyVector	(	std::string const &	name,
		unsigned const	num_comp )

Definition at line 17 of file CellAverageData.cpp.

{
    if (auto const it = cell_averages_.find(name); it != cell_averages_.end())
    {
        auto& prop_vec = *it->second;
        auto const num_comp_mesh = prop_vec.getNumberOfGlobalComponents();
        if (num_comp_mesh == static_cast<int>(num_comp))
        {
            return prop_vec;
        }
 
        OGS_FATAL(
            "The requested property '{}' has {} components, but the one "
            "present in the mesh has {} components.",
            name, num_comp, num_comp_mesh);
    }
 
    auto const name_in_mesh = name + "_avg";
    auto [it, emplaced] = cell_averages_.emplace(
        name, MeshLib::getOrCreateMeshProperty<double>(
                  const_cast<MeshLib::Mesh&>(mesh_), name_in_mesh,
                  MeshLib::MeshItemType::Cell, num_comp));
 
    if (!it->second)
    {
        OGS_FATAL("The cell property '{}' could not be added to the mesh.",
                  name_in_mesh);
    }
 
    if (!emplaced)
    {
        OGS_FATAL(
            "Internal logic error. Something very bad happened. The cell "
            "property '{}' was not added to the list of cell averages to "
            "compute. There is some very strange inconsistency in the "
            "code. Trouble ahead!",
            name_in_mesh);
    }
 
    return *it->second;
}