anonymous_namespace{CollectAndInterpolateNodalDof.cpp} Namespace Reference

Functions
void	collectDofsToMatrixSingleComponentForSomeNodes (MeshLib::Element const &element, std::size_t const mesh_id, NumLib::LocalToGlobalIndexMap const &dof_table, GlobalVector const &x, int const variable, int const component, unsigned const num_nodes, Eigen::Ref< Eigen::VectorXd > all_nodal_dof_for_this_component)

Function Documentation

collectDofsToMatrixSingleComponentForSomeNodes()

void anonymous_namespace{CollectAndInterpolateNodalDof.cpp}::collectDofsToMatrixSingleComponentForSomeNodes	(	MeshLib::Element const &	element,
		std::size_t const	mesh_id,
		NumLib::LocalToGlobalIndexMap const &	dof_table,
		GlobalVector const &	x,
		int const	variable,
		int const	component,
		unsigned const	num_nodes,
		Eigen::Ref< Eigen::VectorXd >	all_nodal_dof_for_this_component )

Collects the degrees of freedom of the passed element from the passed global vector into a vector.

Note

num_nodes can be set "arbitrarily", e.g., to only collect d.o.f. on the base nodes of the passed mesh element.

Definition at line 22 of file CollectAndInterpolateNodalDof.cpp.

{
    bool dof_not_found = false;
 
    for (unsigned element_node_id = 0; element_node_id < num_nodes;
         ++element_node_id)
    {
        auto const& node = *element.getNode(element_node_id);
        auto const node_id = node.getID();
        MeshLib::Location const loc{mesh_id, MeshLib::MeshItemType::Node,
                                    node_id};
        auto const dof_idx = dof_table.getGlobalIndex(loc, variable, component);
 
        if (dof_idx == NumLib::MeshComponentMap::nop)
        {
            // We just skip this d.o.f. Actually we will also skip
            // all other nodes of this mesh element for this (var,
            // comp), because we assume that all linear nodes have a
            // lower node id than any higher order node.
            dof_not_found = true;
        }
        else
        {
            if (dof_not_found)
            {
                // We expect that for all mesh elements all linear
                // nodes have a lower node id than any higher order
                // node. I.e., there are no "holes" in the
                // primary_variables_mat. We rely on there being no
                // "holes" later on when interpolating the nodal
                // d.o.f. to the integration points.
                OGS_FATAL(
                    "This d.o.f. has been found in the d.o.f. "
                    "table, but before some d.o.f. has not been "
                    "found. Something has gone terribly wrong. "
                    "Some assumption in the implementation is "
                    "wrong.");
            }
            all_nodal_dof_for_this_component[element_node_id] = x[dof_idx];
        }
    }
}

References NumLib::LocalToGlobalIndexMap::getGlobalIndex(), MathLib::Point3dWithID::getID(), MeshLib::Element::getNode(), MeshLib::Node, NumLib::MeshComponentMap::nop, and OGS_FATAL.