MeshGeoToolsLib::BoundaryElementsAtPoint Class Reference

Detailed Description

This class collects point elements located at a given point elements.

Definition at line 30 of file BoundaryElementsAtPoint.h.

#include <BoundaryElementsAtPoint.h>

Collaboration diagram for MeshGeoToolsLib::BoundaryElementsAtPoint:

Public Member Functions
	BoundaryElementsAtPoint (MeshLib::Mesh const &mesh, MeshNodeSearcher const &mshNodeSearcher, GeoLib::Point const &point, const bool multiple_nodes_allowed)
	~BoundaryElementsAtPoint ()
GeoLib::Point const &	getPoint () const
std::vector< MeshLib::Element * > const &	getBoundaryElements () const
	Return the vector of boundary elements (i.e. points).

Private Attributes
GeoLib::Point const &	_point
std::vector< MeshLib::Element * >	_boundary_elements

Constructor & Destructor Documentation

BoundaryElementsAtPoint()

MeshGeoToolsLib::BoundaryElementsAtPoint::BoundaryElementsAtPoint	(	MeshLib::Mesh const &	mesh,
		MeshNodeSearcher const &	mshNodeSearcher,
		GeoLib::Point const &	point,
		const bool	multiple_nodes_allowed )

Constructor

Parameters

mesh	a mesh object
mshNodeSearcher	a MeshNodeSearcher object which is internally used to search mesh nodes
point	a point object where edges are searched
multiple_nodes_allowed	Allows to find multiple nodes within the search radius, the nearest node is returned (as point element). This enables to specify larger search radius to find possible other geometries that don't match exactly to the mesh.

Definition at line 26 of file BoundaryElementsAtPoint.cpp.

    : _point(point)
{
    auto const node_ids = mshNodeSearcher.getMeshNodeIDs(_point);
 
#ifdef USE_PETSC
    std::size_t const number_of_found_nodes_at_rank = node_ids.size();
    std::size_t number_of_total_found_nodes = 0;
 
    MPI_Allreduce(&number_of_found_nodes_at_rank, &number_of_total_found_nodes,
                  1, MPI_UNSIGNED_LONG, MPI_SUM, MPI_COMM_WORLD);
 
    if (number_of_total_found_nodes == 0)
    {
        OGS_FATAL(
            "BoundaryElementsAtPoint: the mesh node searcher was unable to "
            "locate the point ({:f}, {:f}, {:f}) in the mesh.",
            _point[0], _point[1], _point[2]);
    }
 
    if (number_of_found_nodes_at_rank == 0)
    {
        return;
    }
#else
    if (node_ids.empty())
    {
        OGS_FATAL(
            "BoundaryElementsAtPoint: the mesh node searcher was unable to "
            "locate the point ({:f}, {:f}, {:f}) in the mesh.",
            _point[0], _point[1], _point[2]);
    }
#endif
 
    if (node_ids.size() == 1)
    {
        std::array<MeshLib::Node*, 1> const nodes = {
            {const_cast<MeshLib::Node*>(mesh.getNode(node_ids[0]))}};
 
        _boundary_elements.push_back(new MeshLib::Point{nodes, node_ids[0]});
        return;
    }
 
    auto& mesh_nodes =
        const_cast<std::vector<MeshLib::Node*>&>(mesh.getNodes());
    std::size_t const nearest_node_id =
        *std::min_element(node_ids.begin(), node_ids.end(),
                          [&mesh_nodes, &point](auto id0, auto id1)
                          {
                              return MathLib::sqrDist(*mesh_nodes[id0], point) <
                                     MathLib::sqrDist(*mesh_nodes[id1], point);
                          });
 
    if (!multiple_nodes_allowed)
    {
        OGS_FATAL(
            "BoundaryElementsAtPoint: the mesh node searcher found {:d} points "
            "near the requested point ({:f}, {:f}, {:f}) in the mesh, while "
            "exactly one is expected. Node  (id={:d}) ({:f}, {:f}, {:f}) has "
            "distance {:f}.",
            node_ids.size(), _point[0], _point[1], _point[2],
            mesh_nodes[nearest_node_id]->getID(),
            (*mesh_nodes[nearest_node_id])[0],
            (*mesh_nodes[nearest_node_id])[1],
            (*mesh_nodes[nearest_node_id])[2],
            MathLib::sqrDist(*mesh_nodes[nearest_node_id], point));
    }
    WARN(
        "BoundaryElementsAtPoint: the mesh node searcher found {:d} points "
        "near the requested point ({:f}, {:f}, {:f}) in the mesh, while "
        "exactly one is expected. Node  (id={:d}) ({:f}, {:f}, {:f}) has "
        "distance {:f}.",
        node_ids.size(), _point[0], _point[1], _point[2],
        mesh_nodes[nearest_node_id]->getID(), (*mesh_nodes[nearest_node_id])[0],
        (*mesh_nodes[nearest_node_id])[1], (*mesh_nodes[nearest_node_id])[2],
        MathLib::sqrDist(*mesh_nodes[nearest_node_id], point));
 
    std::array<MeshLib::Node*, 1> const nodes = {
        {const_cast<MeshLib::Node*>(mesh.getNode(nearest_node_id))}};
 
    _boundary_elements.push_back(new MeshLib::Point{nodes, nearest_node_id});
}

References _boundary_elements, _point, MathLib::Point3dWithID::getID(), MeshGeoToolsLib::MeshNodeSearcher::getMeshNodeIDs(), MeshLib::Mesh::getNode(), MeshLib::Mesh::getNodes(), OGS_FATAL, MathLib::sqrDist(), and WARN().

~BoundaryElementsAtPoint()

MeshGeoToolsLib::BoundaryElementsAtPoint::~BoundaryElementsAtPoint

(

)

Definition at line 111 of file BoundaryElementsAtPoint.cpp.

{
    for (auto p : _boundary_elements)
    {
        delete p;
    }
}

References _boundary_elements.

Member Function Documentation

getBoundaryElements()

std::vector< MeshLib::Element * > const & MeshGeoToolsLib::BoundaryElementsAtPoint::getBoundaryElements ( ) const

inline

Return the vector of boundary elements (i.e. points).

Definition at line 55 of file BoundaryElementsAtPoint.h.

    {
        return _boundary_elements;
    }

References _boundary_elements.

getPoint()

GeoLib::Point const & MeshGeoToolsLib::BoundaryElementsAtPoint::getPoint ( ) const

inline

Definition at line 49 of file BoundaryElementsAtPoint.h.

    {
        return _point;
    }

References _point.

Referenced by MeshGeoToolsLib::BoundaryElementsSearcher::getBoundaryElements().

Member Data Documentation

_boundary_elements

std::vector<MeshLib::Element*> MeshGeoToolsLib::BoundaryElementsAtPoint::_boundary_elements

private

Definition at line 62 of file BoundaryElementsAtPoint.h.

Referenced by BoundaryElementsAtPoint(), ~BoundaryElementsAtPoint(), and getBoundaryElements().

_point

GeoLib::Point const& MeshGeoToolsLib::BoundaryElementsAtPoint::_point

private

Definition at line 61 of file BoundaryElementsAtPoint.h.

Referenced by BoundaryElementsAtPoint(), and getPoint().

---

The documentation for this class was generated from the following files:

• MeshGeoToolsLib/BoundaryElementsAtPoint.h
• MeshGeoToolsLib/BoundaryElementsAtPoint.cpp