OGS
List of all members | Public Member Functions | Private Attributes
MeshGeoToolsLib::BoundaryElementsAtPoint Class Referencefinal

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:
[legend]

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
mesha mesh object
mshNodeSearchera MeshNodeSearcher object which is internally used to search mesh nodes
pointa point object where edges are searched
multiple_nodes_allowedAllows 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 23 of file BoundaryElementsAtPoint.cpp.

26 : _point(point)
27{
28 auto const node_ids = mshNodeSearcher.getMeshNodeIDs(_point);
29
30#ifdef USE_PETSC
31 std::size_t const number_of_found_nodes_at_rank = node_ids.size();
32 std::size_t const number_of_total_found_nodes = BaseLib::MPI::allreduce(
33 number_of_found_nodes_at_rank, MPI_SUM, BaseLib::MPI::Mpi{});
34
35 if (number_of_total_found_nodes == 0)
36 {
37 OGS_FATAL(
38 "BoundaryElementsAtPoint: the mesh node searcher was unable to "
39 "locate the point ({:f}, {:f}, {:f}) in the mesh.",
40 _point[0], _point[1], _point[2]);
41 }
42
43 if (number_of_found_nodes_at_rank == 0)
44 {
45 return;
46 }
47#else
48 if (node_ids.empty())
49 {
50 OGS_FATAL(
51 "BoundaryElementsAtPoint: the mesh node searcher was unable to "
52 "locate the point ({:f}, {:f}, {:f}) in the mesh.",
53 _point[0], _point[1], _point[2]);
54 }
55#endif
56
57 if (node_ids.size() == 1)
58 {
59 std::array<MeshLib::Node*, 1> const nodes = {
60 {const_cast<MeshLib::Node*>(mesh.getNode(node_ids[0]))}};
61
62 _boundary_elements.push_back(new MeshLib::Point{nodes, node_ids[0]});
63 return;
64 }
65
66 auto& mesh_nodes =
67 const_cast<std::vector<MeshLib::Node*>&>(mesh.getNodes());
68 std::size_t const nearest_node_id =
69 *std::min_element(node_ids.begin(), node_ids.end(),
70 [&mesh_nodes, &point](auto id0, auto id1)
71 {
72 return MathLib::sqrDist(*mesh_nodes[id0], point) <
73 MathLib::sqrDist(*mesh_nodes[id1], point);
74 });
75
76 if (!multiple_nodes_allowed)
77 {
78 OGS_FATAL(
79 "BoundaryElementsAtPoint: the mesh node searcher found {:d} points "
80 "near the requested point ({:f}, {:f}, {:f}) in the mesh, while "
81 "exactly one is expected. Node (id={:d}) ({:f}, {:f}, {:f}) has "
82 "distance {:f}.",
83 node_ids.size(), _point[0], _point[1], _point[2],
84 mesh_nodes[nearest_node_id]->getID(),
85 (*mesh_nodes[nearest_node_id])[0],
86 (*mesh_nodes[nearest_node_id])[1],
87 (*mesh_nodes[nearest_node_id])[2],
88 MathLib::sqrDist(*mesh_nodes[nearest_node_id], point));
89 }
90 WARN(
91 "BoundaryElementsAtPoint: the mesh node searcher found {:d} points "
92 "near the requested point ({:f}, {:f}, {:f}) in the mesh, while "
93 "exactly one is expected. Node (id={:d}) ({:f}, {:f}, {:f}) has "
94 "distance {:f}.",
95 node_ids.size(), _point[0], _point[1], _point[2],
96 mesh_nodes[nearest_node_id]->getID(), (*mesh_nodes[nearest_node_id])[0],
97 (*mesh_nodes[nearest_node_id])[1], (*mesh_nodes[nearest_node_id])[2],
98 MathLib::sqrDist(*mesh_nodes[nearest_node_id], point));
99
100 std::array<MeshLib::Node*, 1> const nodes = {
101 {const_cast<MeshLib::Node*>(mesh.getNode(nearest_node_id))}};
102
103 _boundary_elements.push_back(new MeshLib::Point{nodes, nearest_node_id});
104}
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26
WARN
void WARN(fmt::format_string< Args... > fmt, Args &&... args)
Definition Logging.h:40
MathLib::Point3dWithID::getID
std::size_t getID() const
Definition Point3dWithID.h:62
MeshGeoToolsLib::BoundaryElementsAtPoint::_boundary_elements
std::vector< MeshLib::Element * > _boundary_elements
Definition BoundaryElementsAtPoint.h:62
BaseLib::MPI::allreduce
static T allreduce(T const &value, MPI_Op const &mpi_op, Mpi const &mpi)
Definition MPI.h:128
MathLib::sqrDist
double sqrDist(MathLib::Point3d const &p0, MathLib::Point3d const &p1)
Definition Point3d.cpp:26

References _boundary_elements, _point, BaseLib::MPI::allreduce(), 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 106 of file BoundaryElementsAtPoint.cpp.

107{
108 for (auto p : _boundary_elements)
109 {
110 delete p;
111 }
112}

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.

56 {
57 return _boundary_elements;
58 }

References _boundary_elements.

◆ getPoint()

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

Definition at line 49 of file BoundaryElementsAtPoint.h.

50 {
51 return _point;
52 }

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: