anonymous_namespace{AssemblyMixin.cpp} Namespace Reference

Functions
void	checkResiduumNamesVsProcessVariables (std::vector< std::vector< std::string > > const &per_process_residuum_names, std::vector< std::vector< std::reference_wrapper< ProcessLib::ProcessVariable > > > const &per_process_pvs)
std::vector< std::vector< std::reference_wrapper< MeshLib::PropertyVector< double > > > >	createResiduumVectors (MeshLib::Mesh &mesh, std::vector< std::vector< std::string > > const &per_process_residuum_names, std::vector< std::vector< std::reference_wrapper< ProcessLib::ProcessVariable > > > const &per_process_pvs)

Function Documentation

checkResiduumNamesVsProcessVariables()

void anonymous_namespace{AssemblyMixin.cpp}::checkResiduumNamesVsProcessVariables	(	std::vector< std::vector< std::string > > const &	per_process_residuum_names,
		std::vector< std::vector< std::reference_wrapper< ProcessLib::ProcessVariable > > > const &	per_process_pvs )

Definition at line 26 of file AssemblyMixin.cpp.

{
    if (per_process_pvs.size() != per_process_residuum_names.size())
    {
        OGS_FATAL(
            "The number of passed residuum names ({}) does not match the "
            "number of processes ({}).",
            per_process_residuum_names.size(), per_process_pvs.size());
    }
 
    auto check_sizes = [](std::size_t const process_id, auto const& rns_pvs)
    {
        auto const& [rns, pvs] = rns_pvs;
 
        if (rns.size() != pvs.size())
        {
            OGS_FATAL(
                "The number of passed residuum names ({}) does not match the "
                "number of process variables ({}) for process {}.",
                rns.size(), pvs.size(), process_id);
        }
    };
    for (auto const& [process_id, rns_pvs] :
         ranges::views::zip(per_process_residuum_names, per_process_pvs) |
             ranges::views::enumerate)
    {
        check_sizes(process_id, rns_pvs);
    }
}

References OGS_FATAL.

Referenced by createResiduumVectors().

createResiduumVectors()

std::vector< std::vector< std::reference_wrapper< MeshLib::PropertyVector< double > > > > anonymous_namespace{AssemblyMixin.cpp}::createResiduumVectors	(	MeshLib::Mesh &	mesh,
		std::vector< std::vector< std::string > > const &	per_process_residuum_names,
		std::vector< std::vector< std::reference_wrapper< ProcessLib::ProcessVariable > > > const &	per_process_pvs )

Definition at line 62 of file AssemblyMixin.cpp.

{
    checkResiduumNamesVsProcessVariables(per_process_residuum_names,
                                         per_process_pvs);
 
    auto create_mesh_property_for_residuum =
        [&mesh](std::pair<std::string const&,
                          ProcessLib::ProcessVariable const&> const&
                    residuum_name_process_variable)
        -> std::reference_wrapper<MeshLib::PropertyVector<double>>
    {
        auto const& [rn, pv] = residuum_name_process_variable;
        return *MeshLib::getOrCreateMeshProperty<double>(
            mesh, rn, MeshLib::MeshItemType::Node,
            pv.getNumberOfGlobalComponents());
    };
 
    auto create_mesh_properties =
        [&create_mesh_property_for_residuum](auto const& rns_pvs)
    {
        auto const& [rns, pvs] = rns_pvs;
        return ranges::views::zip(rns, pvs) |
               ranges::views::transform(create_mesh_property_for_residuum) |
               ranges::to<std::vector<
                   std::reference_wrapper<MeshLib::PropertyVector<double>>>>;
    };
 
    return ranges::views::zip(per_process_residuum_names, per_process_pvs) |
           ranges::views::transform(create_mesh_properties) |
           ranges::to<std::vector<std::vector<
               std::reference_wrapper<MeshLib::PropertyVector<double>>>>>;
}

References checkResiduumNamesVsProcessVariables(), MeshLib::getOrCreateMeshProperty(), and MeshLib::Node.