MeshLib::anonymous_namespace{ConvertToLinearMesh.cpp} Namespace Reference

Functions
template<typename T_ELEMENT >
T_ELEMENT *	createLinearElement (MeshLib::Element const *e, std::vector< MeshLib::Node * > const &vec_new_nodes)

Function Documentation

createLinearElement()

template<typename T_ELEMENT >

T_ELEMENT* MeshLib::anonymous_namespace{ConvertToLinearMesh.cpp}::createLinearElement	(	MeshLib::Element const *	e,
		std::vector< MeshLib::Node * > const &	vec_new_nodes
	)

Definition at line 31 of file ConvertToLinearMesh.cpp.

{
    auto const n_base_nodes = T_ELEMENT::n_base_nodes;
    auto** nodes = new MeshLib::Node*[n_base_nodes];
    for (unsigned i = 0; i < e->getNumberOfBaseNodes(); i++)
    {
        auto const it = find_if(
            begin(vec_new_nodes), end(vec_new_nodes),
            [node_i = e->getNode(i)](Node* const new_node)
            {
                return *node_i ==
                       *new_node;  // coordinate comparison up to epsilon
            });
        if (it == end(vec_new_nodes))
        {
            OGS_FATAL(
                "A base node {:d} (with original global node id {:d}) not "
                "found in the list for element {:d}.",
                i, e->getNode(i)->getID(), e->getID());
        }
        nodes[i] = const_cast<MeshLib::Node*>(*it);
    }
    return new T_ELEMENT(nodes);
}

References MathLib::Point3dWithID::getID(), MeshLib::Element::getID(), MeshLib::Element::getNode(), MeshLib::Element::getNumberOfBaseNodes(), and OGS_FATAL.