MeshLib::ProjectPointOnMesh Namespace Reference

Functions
Element const *	getProjectedElement (std::vector< const Element * > const &elements, Node const &node)
double	getElevation (Element const &element, Node const &node)

Function Documentation

getElevation()

double MeshLib::ProjectPointOnMesh::getElevation	(	Element const &	element,
		Node const &	node
	)

Returns the z-coordinate of a point projected onto the plane defined by a mesh element.

Definition at line 58 of file ProjectPointOnMesh.cpp.

{
    Eigen::Vector3d const v =
        Eigen::Map<Eigen::Vector3d const>(node.getCoords()) -
        Eigen::Map<Eigen::Vector3d const>(element.getNode(0)->getCoords());
    auto const n = FaceRule::getSurfaceNormal(element).normalized();
    return node[2] - n.dot(v) * n[2];
}

References MathLib::TemplatePoint< T, DIM >::getCoords(), MeshLib::Element::getNode(), and MeshLib::FaceRule::getSurfaceNormal().

Referenced by main(), and voteMatId().

getProjectedElement()

Element const * MeshLib::ProjectPointOnMesh::getProjectedElement	(	std::vector< const Element * > const &	elements,
		Node const &	node
	)

Returns the element in which the given node is located when projected onto a mesh, or nullptr if no such element was found.

Definition at line 21 of file ProjectPointOnMesh.cpp.

{
    auto is_right_of = [&node](Node const& a, Node const& b) {
        return GeoLib::getOrientationFast(node, a, b) ==
               GeoLib::Orientation::CW;
    };
 
    for (auto const* e : elements)
    {
        auto const* nodes = e->getNodes();
        if (e->getGeomType() == MeshElemType::TRIANGLE)
        {
            auto const& a = *nodes[0];
            auto const& b = *nodes[1];
            auto const& c = *nodes[2];
            if (!is_right_of(a, b) && !is_right_of(b, c) && !is_right_of(c, a))
            {
                return e;
            }
        }
        else if (e->getGeomType() == MeshElemType::QUAD)
        {
            auto const& a = *nodes[0];
            auto const& b = *nodes[1];
            auto const& c = *nodes[2];
            auto const& d = *nodes[3];
            if (!is_right_of(a, b) && !is_right_of(b, c) &&
                !is_right_of(c, d) && !is_right_of(d, a))
            {
                return e;
            }
        }
    }
    return nullptr;
}

References MaterialPropertyLib::c, GeoLib::CW, GeoLib::getOrientationFast(), MeshLib::QUAD, and MeshLib::TRIANGLE.

Referenced by getProjectedElement(), and main().