Special adaptor that works with subdomains of the Mesh. More...

#include <NanoflannMeshAdaptor.h>

Public Types
typedef Real	coord_t
	libMesh `Point` coordinate type More...

Public Member Functions
	NanoflannMeshSubdomainAdaptor (const MeshBase &mesh, subdomain_id_type s)

size_t	kdtree_get_point_count () const
	Must return the number of data points. More...

coord_t	kdtree_distance (const coord_t *p1, const size_t idx_p2, size_t size) const
	Returns the distance between the vector "p1[0:size-1]" and the data point with index "idx_p2" stored in _mesh. More...

coord_t	kdtree_get_pt (const size_t idx, int dim) const
	Returns the dim'th component of the idx'th point in the class. More...

template<class BBOX >
bool	kdtree_get_bbox (BBOX &) const
	Optional bounding-box computation: return false to default to a standard bbox computation loop. More...

Private Attributes
const MeshBase &	_mesh

subdomain_id_type	_sid

std::set< dof_id_type >	_legal_point_indices

Detailed Description

template<unsigned int Dim>
class NanoflannMeshSubdomainAdaptor< Dim >

Special adaptor that works with subdomains of the Mesh.

When nanoflann asks for the distance between points p1 and p2, and p2 is not a Point of an element in the required subdomain, a "large" distance is returned. Likewise, when the coordinates of a point not in the subdomain are requested, a point at "infinity" is returned.

Definition at line 118 of file NanoflannMeshAdaptor.h.

Member Typedef Documentation

◆ coord_t

template<unsigned int Dim>

typedef Real NanoflannMeshSubdomainAdaptor< Dim >::coord_t

libMesh Point coordinate type

Definition at line 144 of file NanoflannMeshAdaptor.h.

Constructor & Destructor Documentation

◆ NanoflannMeshSubdomainAdaptor()

template<unsigned int Dim>

NanoflannMeshSubdomainAdaptor< Dim >::NanoflannMeshSubdomainAdaptor	(	const MeshBase &	mesh,
		subdomain_id_type	s
	)

inline

Definition at line 131 of file NanoflannMeshAdaptor.h.

                                                                             : _mesh(mesh), _sid(s)
   {
     // Loop over the elements of the Mesh, for those in the requested
     // subdomain, add its node ids to the _legal_point_indices set.
     for (const auto & elem : _mesh.active_element_ptr_range())
       if (elem->subdomain_id() == _sid)
         for (MooseIndex(elem->n_vertices()) n = 0; n < elem->n_vertices(); ++n)
           _legal_point_indices.insert(elem->node_id(n));
   }

Member Function Documentation

◆ kdtree_distance()

template<unsigned int Dim>

coord_t NanoflannMeshSubdomainAdaptor< Dim >::kdtree_distance	(	const coord_t *	p1,
		const size_t	idx_p2,
		size_t	size
	)		const

inline

Returns the distance between the vector "p1[0:size-1]" and the data point with index "idx_p2" stored in _mesh.

Definition at line 155 of file NanoflannMeshAdaptor.h.

   {
     libmesh_assert_equal_to(size, Dim);
 
     // If this is not a valid point, then return a "large" distance.
     if (!_legal_point_indices.count(static_cast<dof_id_type>(idx_p2)))
       return std::numeric_limits<coord_t>::max();
 
     // Construct a libmesh Point object from the input coord_t.  This
     // assumes LIBMESH_DIM==3.
     Point point1(p1[0], size > 1 ? p1[1] : 0., size > 2 ? p1[2] : 0.);
 
     // Get the referred-to point from the Mesh
     const Point & point2 = _mesh.point(idx_p2);
 
     // Compute Euclidean distance, squared
     return (point1 - point2).norm_sq();
   }

◆ kdtree_get_bbox()

template<unsigned int Dim>

template<class BBOX >

bool NanoflannMeshSubdomainAdaptor< Dim >::kdtree_get_bbox ( BBOX & ) const

inline

Optional bounding-box computation: return false to default to a standard bbox computation loop.

Definition at line 195 of file NanoflannMeshAdaptor.h.

   {
     return false;
   }

◆ kdtree_get_point_count()

template<unsigned int Dim>

size_t NanoflannMeshSubdomainAdaptor< Dim >::kdtree_get_point_count ( ) const

inline

Must return the number of data points.

Definition at line 149 of file NanoflannMeshAdaptor.h.

149 { return _mesh.n_nodes(); }

NanoflannMeshSubdomainAdaptor::_mesh

const MeshBase & _mesh

Definition: NanoflannMeshAdaptor.h:122

libMesh::MeshBase::n_nodes

virtual dof_id_type n_nodes() const=0

◆ kdtree_get_pt()

template<unsigned int Dim>

coord_t NanoflannMeshSubdomainAdaptor< Dim >::kdtree_get_pt	(	const size_t	idx,
		int	dim
	)		const

inline

Returns the dim'th component of the idx'th point in the class.

Definition at line 177 of file NanoflannMeshAdaptor.h.

   {
     libmesh_assert_less(dim, (int)Dim);
     libmesh_assert_less(idx, _mesh.n_nodes());
     libmesh_assert_less(dim, 3);
 
     // If this is not a valid point, then return a "large" distance.
     if (!_legal_point_indices.count(static_cast<dof_id_type>(idx)))
       return std::numeric_limits<coord_t>::max();
 
     return _mesh.point(idx)(dim);
   }