MeshNodesOnPoint.cpp

Go to the documentation of this file.

 
#include "MeshNodesOnPoint.h"
 
#include "MeshLib/Mesh.h"
 
namespace MeshGeoToolsLib
{
MeshNodesOnPoint::MeshNodesOnPoint(MeshLib::Mesh const& mesh,
                                   GeoLib::Grid<MeshLib::Node> const& mesh_grid,
                                   GeoLib::Point const& pnt,
                                   double epsilon_radius,
                                   SearchAllNodes search_all_nodes)
    : _mesh(mesh), _pnt(pnt)
{
    std::vector<std::size_t> vec_ids(
        mesh_grid.getPointsInEpsilonEnvironment(pnt, epsilon_radius));
    if (search_all_nodes == SearchAllNodes::Yes)
    {
        _msh_node_ids = vec_ids;
    }
    else
    {
        for (auto id : vec_ids)
        {
            if (MeshLib::isBaseNode(*_mesh.getNode(id),
                                    _mesh.getElementsConnectedToNode(id)))
            {
                _msh_node_ids.push_back(id);
            }
        }
    }
}
 
}  // end namespace MeshGeoToolsLib