libMesh::PointLocatorTree Class Reference

#include <point_locator_tree.h>

Inheritance diagram for libMesh::PointLocatorTree:

Public Member Functions

 PointLocatorTree (const MeshBase &mesh, const PointLocatorBase *master=nullptr)
 
 PointLocatorTree (const MeshBase &mesh, const Trees::BuildType build_type, const PointLocatorBase *master=nullptr)
 
 ~PointLocatorTree ()
 
virtual void clear () override
 
void init (Trees::BuildType build_type)
 
virtual void init () override
 
virtual const Elemoperator() (const Point &p, const std::set< subdomain_id_type > *allowed_subdomains=nullptr) const override
 
virtual void operator() (const Point &p, std::set< const Elem *> &candidate_elements, const std::set< subdomain_id_type > *allowed_subdomains=nullptr) const override
 
const Elemperform_linear_search (const Point &p, const std::set< subdomain_id_type > *allowed_subdomains, bool use_close_to_point, Real close_to_point_tolerance=TOLERANCE) const
 
std::set< const Elem * > perform_fuzzy_linear_search (const Point &p, const std::set< subdomain_id_type > *allowed_subdomains, Real close_to_point_tolerance=TOLERANCE) const
 
virtual void enable_out_of_mesh_mode () override
 
virtual void disable_out_of_mesh_mode () override
 
void set_target_bin_size (unsigned int target)
 
unsigned int get_target_bin_size () const
 
virtual const Nodelocate_node (const Point &p, const std::set< subdomain_id_type > *allowed_subdomains=nullptr, Real tol=TOLERANCE) const
 
bool initialized () const
 
virtual void set_close_to_point_tol (Real close_to_point_tol)
 
virtual void unset_close_to_point_tol ()
 
const MeshBaseget_mesh () const
 

Static Public Member Functions

static std::unique_ptr< PointLocatorBasebuild (PointLocatorType t, const MeshBase &mesh, const PointLocatorBase *master=nullptr)
 
static std::string get_info ()
 
static void print_info (std::ostream &out=libMesh::out)
 
static unsigned int n_objects ()
 
static void enable_print_counter_info ()
 
static void disable_print_counter_info ()
 

Public Attributes

bool _verbose
 

Protected Types

typedef std::map< std::string, std::pair< unsigned int, unsigned int > > Counts
 

Protected Member Functions

void increment_constructor_count (const std::string &name)
 
void increment_destructor_count (const std::string &name)
 

Protected Attributes

TreeBase_tree
 
const Elem_element
 
bool _out_of_mesh_mode
 
unsigned int _target_bin_size
 
Trees::BuildType _build_type
 
const PointLocatorBase_master
 
const MeshBase_mesh
 
bool _initialized
 
bool _use_close_to_point_tol
 
Real _close_to_point_tol
 

Static Protected Attributes

static Counts _counts
 
static Threads::atomic< unsigned int > _n_objects
 
static Threads::spin_mutex _mutex
 
static bool _enable_print_counter = true
 

Detailed Description

This is a point locator. It locates points in space using a tree: given a mesh they return the element and local coordinates for a given point in global coordinates. Use PointLocatorBase::build() to create objects of this type at run time.

Author
Daniel Dreyer
Date
2003

Definition at line 48 of file point_locator_tree.h.

Member Typedef Documentation

◆ Counts

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts
protectedinherited

Data structure to log the information. The log is identified by the class name.

Definition at line 117 of file reference_counter.h.

Constructor & Destructor Documentation

◆ PointLocatorTree() [1/2]

libMesh::PointLocatorTree::PointLocatorTree ( const MeshBase mesh,
const PointLocatorBase master = nullptr 
)

Constructor. Needs the mesh in which the points should be located. Optionally takes a master interpolator. This master helps in saving memory by reducing the number of trees in use. Only the master locator holds a tree, the others simply use the master's tree.

Definition at line 37 of file point_locator_tree.C.

References _build_type, and init().

38  :
39  PointLocatorBase (mesh,master),
40  _tree (nullptr),
41  _element (nullptr),
42  _out_of_mesh_mode(false),
43  _target_bin_size (200),
45 {
46  this->init(_build_type);
47 }
virtual void init() override
MeshBase & mesh
PointLocatorBase(const MeshBase &mesh, const PointLocatorBase *master)

◆ PointLocatorTree() [2/2]

libMesh::PointLocatorTree::PointLocatorTree ( const MeshBase mesh,
const Trees::BuildType  build_type,
const PointLocatorBase master = nullptr 
)

Constructor. Needs the mesh in which the points should be located. Allows the user to specify the method to use when building the tree. Optionally takes a master interpolator. This master helps in saving memory by reducing the number of trees in use. Only the master locator holds a tree, the others simply use the master's tree. Allows the user to specify the build type.

Definition at line 51 of file point_locator_tree.C.

References _build_type, and init().

53  :
54  PointLocatorBase (mesh,master),
55  _tree (nullptr),
56  _element (nullptr),
57  _out_of_mesh_mode(false),
58  _target_bin_size (200),
59  _build_type(build_type)
60 {
61  this->init(_build_type);
62 }
virtual void init() override
MeshBase & mesh
PointLocatorBase(const MeshBase &mesh, const PointLocatorBase *master)

◆ ~PointLocatorTree()

libMesh::PointLocatorTree::~PointLocatorTree ( )

Destructor.

Definition at line 66 of file point_locator_tree.C.

References clear().

67 {
68  this->clear ();
69 }
virtual void clear() override

Member Function Documentation

◆ build()

std::unique_ptr< PointLocatorBase > libMesh::PointLocatorBase::build ( PointLocatorType  t,
const MeshBase mesh,
const PointLocatorBase master = nullptr 
)
staticinherited

Builds an PointLocator for the mesh mesh. Optionally takes a master PointLocator to save memory. An std::unique_ptr<PointLocatorBase> is returned to prevent memory leak. This way the user need not remember to delete the object.

Definition at line 62 of file point_locator_base.C.

References libMesh::Trees::ELEMENTS, libMesh::Trees::LOCAL_ELEMENTS, mesh, libMesh::TREE, libMesh::TREE_ELEMENTS, and libMesh::TREE_LOCAL_ELEMENTS.

Referenced by libMesh::MeshBase::point_locator(), and libMesh::MeshBase::sub_point_locator().

65 {
66  switch (t)
67  {
68  case TREE:
69  return libmesh_make_unique<PointLocatorTree>(mesh, /*Trees::NODES,*/ master);
70 
71  case TREE_ELEMENTS:
72  return libmesh_make_unique<PointLocatorTree>(mesh, Trees::ELEMENTS, master);
73 
75  return libmesh_make_unique<PointLocatorTree>(mesh, Trees::LOCAL_ELEMENTS, master);
76 
77  default:
78  libmesh_error_msg("ERROR: Bad PointLocatorType = " << t);
79  }
80 }
MeshBase & mesh

◆ clear()

void libMesh::PointLocatorTree::clear ( )
overridevirtual

Clears the locator. This function frees dynamic memory with "delete".

Implements libMesh::PointLocatorBase.

Definition at line 73 of file point_locator_tree.C.

References libMesh::PointLocatorBase::_initialized, libMesh::PointLocatorBase::_master, and _tree.

Referenced by ~PointLocatorTree().

74 {
75  // only delete the tree when we are the master
76  if (this->_tree != nullptr)
77  {
78  if (this->_master == nullptr)
79  // we own the tree
80  delete this->_tree;
81  else
82  // someone else owns and therefore deletes the tree
83  this->_tree = nullptr;
84 
85  // make sure operator () throws an assertion
86  this->_initialized = false;
87  }
88 }
const PointLocatorBase * _master

◆ disable_out_of_mesh_mode()

void libMesh::PointLocatorTree::disable_out_of_mesh_mode ( void  )
overridevirtual

Disables out-of-mesh mode (default). If asked to find a point that is contained in no mesh at all, the point locator will now crash.

Implements libMesh::PointLocatorBase.

Definition at line 339 of file point_locator_tree.C.

References _out_of_mesh_mode.

340 {
341  _out_of_mesh_mode = false;
342 }

◆ disable_print_counter_info()

void libMesh::ReferenceCounter::disable_print_counter_info ( )
staticinherited

Definition at line 106 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit().

107 {
108  _enable_print_counter = false;
109  return;
110 }

◆ enable_out_of_mesh_mode()

void libMesh::PointLocatorTree::enable_out_of_mesh_mode ( )
overridevirtual

Enables out-of-mesh mode. In this mode, if asked to find a point that is contained in no mesh at all, the point locator will return nullptr instead of crashing. Per default, this mode is off.

Implements libMesh::PointLocatorBase.

Definition at line 331 of file point_locator_tree.C.

References _out_of_mesh_mode.

332 {
333  // Out-of-mesh mode should now work properly even on meshes with
334  // non-affine elements.
335  _out_of_mesh_mode = true;
336 }

◆ enable_print_counter_info()

void libMesh::ReferenceCounter::enable_print_counter_info ( )
staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 100 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

101 {
102  _enable_print_counter = true;
103  return;
104 }

◆ get_info()

std::string libMesh::ReferenceCounter::get_info ( )
staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

48 {
49 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
50 
51  std::ostringstream oss;
52 
53  oss << '\n'
54  << " ---------------------------------------------------------------------------- \n"
55  << "| Reference count information |\n"
56  << " ---------------------------------------------------------------------------- \n";
57 
58  for (const auto & pr : _counts)
59  {
60  const std::string name(pr.first);
61  const unsigned int creations = pr.second.first;
62  const unsigned int destructions = pr.second.second;
63 
64  oss << "| " << name << " reference count information:\n"
65  << "| Creations: " << creations << '\n'
66  << "| Destructions: " << destructions << '\n';
67  }
68 
69  oss << " ---------------------------------------------------------------------------- \n";
70 
71  return oss.str();
72 
73 #else
74 
75  return "";
76 
77 #endif
78 }
std::string name(const ElemQuality q)
Definition: elem_quality.C:42

◆ get_mesh()

const MeshBase & libMesh::PointLocatorBase::get_mesh ( ) const
inherited

Get a const reference to this PointLocator's mesh.

Definition at line 96 of file point_locator_base.C.

References libMesh::PointLocatorBase::_mesh.

Referenced by libMesh::PeriodicBoundaries::neighbor().

97 {
98  return _mesh;
99 }

◆ get_target_bin_size()

unsigned int libMesh::PointLocatorTree::get_target_bin_size ( ) const

Get the target bin size.

Definition at line 351 of file point_locator_tree.C.

References _target_bin_size.

Referenced by init().

352 {
353  return _target_bin_size;
354 }

◆ increment_constructor_count()

void libMesh::ReferenceCounter::increment_constructor_count ( const std::string &  name)
inlineprotectedinherited

Increments the construction counter. Should be called in the constructor of any derived class that will be reference counted.

Definition at line 181 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

182 {
183  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
184  std::pair<unsigned int, unsigned int> & p = _counts[name];
185 
186  p.first++;
187 }
std::string name(const ElemQuality q)
Definition: elem_quality.C:42
spin_mutex spin_mtx
Definition: threads.C:29

◆ increment_destructor_count()

void libMesh::ReferenceCounter::increment_destructor_count ( const std::string &  name)
inlineprotectedinherited

Increments the destruction counter. Should be called in the destructor of any derived class that will be reference counted.

Definition at line 194 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

195 {
196  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
197  std::pair<unsigned int, unsigned int> & p = _counts[name];
198 
199  p.second++;
200 }
std::string name(const ElemQuality q)
Definition: elem_quality.C:42
spin_mutex spin_mtx
Definition: threads.C:29

◆ init() [1/2]

void libMesh::PointLocatorTree::init ( Trees::BuildType  build_type)

Initializes the locator, so that the operator() methods can be used. This function allocates dynamic memory with "new".

Definition at line 99 of file point_locator_tree.C.

References _build_type, _element, libMesh::PointLocatorBase::_initialized, libMesh::PointLocatorBase::_master, libMesh::PointLocatorBase::_mesh, _tree, std::abs(), libMesh::MeshTools::create_bounding_box(), libMesh::err, get_target_bin_size(), libMesh::PointLocatorBase::initialized(), libMesh::MeshBase::mesh_dimension(), and libMesh::Real.

100 {
101  libmesh_assert (!this->_tree);
102 
103  if (this->_initialized)
104  {
105  // Warn that we are already initialized
106  libMesh::err << "Warning: PointLocatorTree already initialized! Will ignore this call..." << std::endl;
107 
108  // Further warn if we try to init() again with a different build_type
109  if (_build_type != build_type)
110  {
111  libMesh::err << "Warning: PointLocatorTree is using build_type = " << _build_type << ".\n"
112  << "Your requested build_type, " << build_type << " will not be used!" << std::endl;
113  }
114  }
115 
116  else
117  {
118  // Let the requested build_type override the _build_type we were
119  // constructed with. This is no big deal since we have not been
120  // initialized before.
121  _build_type = build_type;
122 
123  if (this->_master == nullptr)
124  {
125  LOG_SCOPE("init(no master)", "PointLocatorTree");
126 
127  if (this->_mesh.mesh_dimension() == 3)
129  else
130  {
131  // A 1D/2D mesh in 3D space needs special consideration.
132  // If the mesh is planar XY, we want to build a QuadTree
133  // to search efficiently. If the mesh is truly a manifold,
134  // then we need an octree
135 #if LIBMESH_DIM > 2
136  bool is_planar_xy = false;
137 
138  // Build the bounding box for the mesh. If the delta-z bound is
139  // negligibly small then we can use a quadtree.
140  {
141  BoundingBox bbox = MeshTools::create_bounding_box(this->_mesh);
142 
143  const Real
144  Dx = bbox.second(0) - bbox.first(0),
145  Dz = bbox.second(2) - bbox.first(2);
146 
147  if (std::abs(Dz/(Dx + 1.e-20)) < 1e-10)
148  is_planar_xy = true;
149  }
150 
151  if (!is_planar_xy)
153  else
154 #endif
155 #if LIBMESH_DIM > 1
157 #else
159 #endif
160  }
161  }
162 
163  else
164  {
165  // We are _not_ the master. Let our Tree point to
166  // the master's tree. But for this we first transform
167  // the master in a state for which we are friends.
168  // And make sure the master has a tree!
169  const PointLocatorTree * my_master =
170  cast_ptr<const PointLocatorTree *>(this->_master);
171 
172  if (my_master->initialized())
173  this->_tree = my_master->_tree;
174  else
175  libmesh_error_msg("ERROR: Initialize master first, then servants!");
176  }
177 
178  // Not all PointLocators may own a tree, but all of them
179  // use their own element pointer. Let the element pointer
180  // be unique for every interpolator.
181  // Suppose the interpolators are used concurrently
182  // at different locations in the mesh, then it makes quite
183  // sense to have unique start elements.
184  this->_element = nullptr;
185  }
186 
187  // ready for take-off
188  this->_initialized = true;
189 }
double abs(double a)
libMesh::BoundingBox create_bounding_box(const MeshBase &mesh)
Definition: mesh_tools.C:386
Tree< 8 > OctTree
Definition: tree.h:134
unsigned int get_target_bin_size() const
Tree< 4 > QuadTree
Definition: tree.h:128
PointLocatorTree(const MeshBase &mesh, const PointLocatorBase *master=nullptr)
OStreamProxy err(std::cerr)
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
unsigned int mesh_dimension() const
Definition: mesh_base.C:126
const PointLocatorBase * _master
Tree< 2 > BinaryTree
Definition: tree.h:122

◆ init() [2/2]

void libMesh::PointLocatorTree::init ( )
overridevirtual

Initializes the locator, so that the operator() methods can be used. This function allocates dynamic memory with "new".

Implements libMesh::PointLocatorBase.

Definition at line 92 of file point_locator_tree.C.

References _build_type.

Referenced by PointLocatorTree().

93 {
94  this->init(_build_type);
95 }
virtual void init() override

◆ initialized()

bool libMesh::PointLocatorBase::initialized ( ) const
inherited
Returns
true when this object is properly initialized and ready for use, false otherwise.

Definition at line 55 of file point_locator_base.C.

References libMesh::PointLocatorBase::_initialized.

Referenced by init().

56 {
57  return this->_initialized;
58 }

◆ locate_node()

const Node * libMesh::PointLocatorBase::locate_node ( const Point p,
const std::set< subdomain_id_type > *  allowed_subdomains = nullptr,
Real  tol = TOLERANCE 
) const
virtualinherited
Returns
A pointer to a Node with global coordinates p or nullptr if no such Node can be found.

Virtual subclasses can override for efficiency, but the base class has a default implementation that works based on element lookup.

Optionally allows the user to restrict the subdomains searched; with such a restriction, only a Node belonging to an element on one or more of those subdomains will be returned.

Will only return a Node whose distance from p is less than tol multiplied by the size of a semilocal element which contains p.

Definition at line 104 of file point_locator_base.C.

References libMesh::TensorTools::norm_sq(), libMesh::PointLocatorBase::operator()(), and libMesh::Real.

107 {
108  std::set<const Elem *> candidate_elements;
109  this->operator()(p, candidate_elements, allowed_subdomains);
110 
111  for (const auto & elem : candidate_elements)
112  {
113  const int elem_n_nodes = elem->n_nodes();
114  const Real hmax = elem->hmax();
115  const Real dist_tol_sq = (tol * hmax) * (tol * hmax);
116 
117  for (int n=0; n != elem_n_nodes; ++n)
118  if ((elem->point(n) - p).norm_sq() < dist_tol_sq)
119  return elem->node_ptr(n);
120  }
121 
122  return nullptr;
123 }
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual const Elem * operator()(const Point &p, const std::set< subdomain_id_type > *allowed_subdomains=nullptr) const =0

◆ n_objects()

static unsigned int libMesh::ReferenceCounter::n_objects ( )
inlinestaticinherited

Prints the number of outstanding (created, but not yet destroyed) objects.

Definition at line 83 of file reference_counter.h.

References libMesh::ReferenceCounter::_n_objects.

84  { return _n_objects; }
static Threads::atomic< unsigned int > _n_objects

◆ operator()() [1/2]

const Elem * libMesh::PointLocatorTree::operator() ( const Point p,
const std::set< subdomain_id_type > *  allowed_subdomains = nullptr 
) const
overridevirtual

Locates the element in which the point with global coordinates p is located, optionally restricted to a set of allowed subdomains. The mutable _element member is used to cache the result and allow it to be used during the next call to operator().

Implements libMesh::PointLocatorBase.

Definition at line 193 of file point_locator_tree.C.

References libMesh::PointLocatorBase::_close_to_point_tol, _element, libMesh::PointLocatorBase::_initialized, _out_of_mesh_mode, _tree, libMesh::PointLocatorBase::_use_close_to_point_tol, libMesh::PointLocatorBase::_verbose, libMesh::Elem::active(), libMesh::Elem::contains_point(), libMesh::TreeBase::find_element(), libMesh::out, and perform_linear_search().

195 {
196  libmesh_assert (this->_initialized);
197 
198  LOG_SCOPE("operator()", "PointLocatorTree");
199 
200  // If we're provided with an allowed_subdomains list and have a cached element, make sure it complies
201  if (allowed_subdomains && this->_element && !allowed_subdomains->count(this->_element->subdomain_id())) this->_element = nullptr;
202 
203  // First check the element from last time before asking the tree
204  if (this->_element==nullptr || !(this->_element->contains_point(p)))
205  {
206  // ask the tree
207  this->_element = this->_tree->find_element (p,allowed_subdomains);
208 
209  if (this->_element == nullptr)
210  {
211  // If we haven't found the element, we may want to do a linear
212  // search using a tolerance.
214  {
215  if (_verbose)
216  {
217  libMesh::out << "Performing linear search using close-to-point tolerance "
219  << std::endl;
220  }
221 
222  this->_element =
223  this->perform_linear_search(p,
224  allowed_subdomains,
225  /*use_close_to_point*/ true,
227 
228  return this->_element;
229  }
230 
231  // No element seems to contain this point. In theory, our
232  // tree now correctly handles curved elements. In
233  // out-of-mesh mode this is sometimes expected, and we can
234  // just return nullptr without searching further. Out of
235  // out-of-mesh mode, something must have gone wrong.
236  libmesh_assert_equal_to (_out_of_mesh_mode, true);
237 
238  return this->_element;
239  }
240  }
241 
242  // If we found an element, it should be active
243  libmesh_assert (!this->_element || this->_element->active());
244 
245  // If we found an element and have a restriction list, they better match
246  libmesh_assert (!this->_element || !allowed_subdomains || allowed_subdomains->count(this->_element->subdomain_id()));
247 
248  // return the element
249  return this->_element;
250 }
virtual bool contains_point(const Point &p, Real tol=TOLERANCE) const
Definition: elem.C:2301
virtual const Elem * find_element(const Point &p, const std::set< subdomain_id_type > *allowed_subdomains=nullptr, Real relative_tol=TOLERANCE) const =0
const Elem * perform_linear_search(const Point &p, const std::set< subdomain_id_type > *allowed_subdomains, bool use_close_to_point, Real close_to_point_tolerance=TOLERANCE) const
OStreamProxy out(std::cout)
bool active() const
Definition: elem.h:2390

◆ operator()() [2/2]

void libMesh::PointLocatorTree::operator() ( const Point p,
std::set< const Elem *> &  candidate_elements,
const std::set< subdomain_id_type > *  allowed_subdomains = nullptr 
) const
overridevirtual

Locates a set of elements in proximity to the point with global coordinates p Pure virtual. Optionally allows the user to restrict the subdomains searched.

Implements libMesh::PointLocatorBase.

Definition at line 253 of file point_locator_tree.C.

References libMesh::PointLocatorBase::_close_to_point_tol, libMesh::PointLocatorBase::_initialized, and perform_fuzzy_linear_search().

256 {
257  libmesh_assert (this->_initialized);
258 
259  LOG_SCOPE("operator() - Version 2", "PointLocatorTree");
260 
261  // forward call to perform_linear_search
262  candidate_elements = this->perform_fuzzy_linear_search(p, allowed_subdomains, _close_to_point_tol);
263 }
std::set< const Elem * > perform_fuzzy_linear_search(const Point &p, const std::set< subdomain_id_type > *allowed_subdomains, Real close_to_point_tolerance=TOLERANCE) const

◆ perform_fuzzy_linear_search()

std::set< const Elem * > libMesh::PointLocatorTree::perform_fuzzy_linear_search ( const Point p,
const std::set< subdomain_id_type > *  allowed_subdomains,
Real  close_to_point_tolerance = TOLERANCE 
) const

A method to check if "fat" point p is in multiple elements. This would happen if p is close to a face or node. This is important for evaluating MeshFunction on faces when discontinuous shape functions are used.

Definition at line 305 of file point_locator_tree.C.

References _build_type, libMesh::PointLocatorBase::_mesh, libMesh::MeshBase::active_element_ptr_range(), libMesh::MeshBase::active_local_element_ptr_range(), and libMesh::Trees::LOCAL_ELEMENTS.

Referenced by operator()().

308 {
309  LOG_SCOPE("perform_fuzzy_linear_search", "PointLocatorTree");
310 
311  std::set<const Elem *> candidate_elements;
312 
313  // The type of iterator depends on the Trees::BuildType
314  // used for this PointLocator. If it's
315  // TREE_LOCAL_ELEMENTS, we only want to double check
316  // local elements during this linear search.
317  SimpleRange<MeshBase::const_element_iterator> r =
321 
322  for (const auto & elem : r)
323  if ((!allowed_subdomains || allowed_subdomains->count(elem->subdomain_id())) && elem->close_to_point(p, close_to_point_tolerance))
324  candidate_elements.insert(elem);
325 
326  return candidate_elements;
327 }
virtual SimpleRange< element_iterator > active_element_ptr_range()=0
virtual SimpleRange< element_iterator > active_local_element_ptr_range()=0

◆ perform_linear_search()

const Elem * libMesh::PointLocatorTree::perform_linear_search ( const Point p,
const std::set< subdomain_id_type > *  allowed_subdomains,
bool  use_close_to_point,
Real  close_to_point_tolerance = TOLERANCE 
) const

As a fallback option, it's helpful to be able to do a linear search over the entire mesh. This can be used if operator() fails to find an element that contains p, for example. Optionally specify a "close to point" tolerance to use in the linear search. Return nullptr if no element is found.

Definition at line 267 of file point_locator_tree.C.

References _build_type, libMesh::PointLocatorBase::_mesh, libMesh::MeshBase::active_element_ptr_range(), libMesh::MeshBase::active_local_element_ptr_range(), and libMesh::Trees::LOCAL_ELEMENTS.

Referenced by operator()().

271 {
272  LOG_SCOPE("perform_linear_search", "PointLocatorTree");
273 
274  // The type of iterator depends on the Trees::BuildType
275  // used for this PointLocator. If it's
276  // TREE_LOCAL_ELEMENTS, we only want to double check
277  // local elements during this linear search.
278  SimpleRange<MeshBase::const_element_iterator> r =
282 
283  for (const auto & elem : r)
284  {
285  if (!allowed_subdomains ||
286  allowed_subdomains->count(elem->subdomain_id()))
287  {
288  if (!use_close_to_point)
289  {
290  if (elem->contains_point(p))
291  return elem;
292  }
293  else
294  {
295  if (elem->close_to_point(p, close_to_point_tolerance))
296  return elem;
297  }
298  }
299  }
300 
301  return nullptr;
302 }
virtual SimpleRange< element_iterator > active_element_ptr_range()=0
virtual SimpleRange< element_iterator > active_local_element_ptr_range()=0

◆ print_info()

void libMesh::ReferenceCounter::print_info ( std::ostream &  out = libMesh::out)
staticinherited

Prints the reference information, by default to libMesh::out.

Definition at line 87 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter, and libMesh::ReferenceCounter::get_info().

88 {
90  out_stream << ReferenceCounter::get_info();
91 }
static std::string get_info()

◆ set_close_to_point_tol()

void libMesh::PointLocatorBase::set_close_to_point_tol ( Real  close_to_point_tol)
virtualinherited

Set a tolerance to use when determining if a point is contained within the mesh.

Definition at line 82 of file point_locator_base.C.

References libMesh::PointLocatorBase::_close_to_point_tol, and libMesh::PointLocatorBase::_use_close_to_point_tol.

Referenced by libMesh::MeshFunction::set_point_locator_tolerance().

83 {
85  _close_to_point_tol = close_to_point_tol;
86 }

◆ set_target_bin_size()

void libMesh::PointLocatorTree::set_target_bin_size ( unsigned int  target)

Set the target bin size.

Definition at line 345 of file point_locator_tree.C.

References _target_bin_size.

346 {
347  _target_bin_size = target_bin_size;
348 }

◆ unset_close_to_point_tol()

void libMesh::PointLocatorBase::unset_close_to_point_tol ( )
virtualinherited

Specify that we do not want to use a user-specified tolerance to determine if a point is contained within the mesh.

Definition at line 89 of file point_locator_base.C.

References libMesh::PointLocatorBase::_close_to_point_tol, libMesh::PointLocatorBase::_use_close_to_point_tol, and libMesh::TOLERANCE.

Referenced by libMesh::MeshFunction::unset_point_locator_tolerance().

Member Data Documentation

◆ _build_type

Trees::BuildType libMesh::PointLocatorTree::_build_type
protected

How the underlying tree is built.

Definition at line 194 of file point_locator_tree.h.

Referenced by init(), perform_fuzzy_linear_search(), perform_linear_search(), and PointLocatorTree().

◆ _close_to_point_tol

Real libMesh::PointLocatorBase::_close_to_point_tol
protectedinherited

◆ _counts

ReferenceCounter::Counts libMesh::ReferenceCounter::_counts
staticprotectedinherited

◆ _element

const Elem* libMesh::PointLocatorTree::_element
mutableprotected

Pointer to the last element that was found by the tree. Chances are that this may be close to the next call to operator()...

Definition at line 178 of file point_locator_tree.h.

Referenced by init(), and operator()().

◆ _enable_print_counter

bool libMesh::ReferenceCounter::_enable_print_counter = true
staticprotectedinherited

Flag to control whether reference count information is printed when print_info is called.

Definition at line 141 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), and libMesh::ReferenceCounter::print_info().

◆ _initialized

bool libMesh::PointLocatorBase::_initialized
protectedinherited

true when properly initialized, false otherwise.

Definition at line 196 of file point_locator_base.h.

Referenced by clear(), init(), libMesh::PointLocatorBase::initialized(), and operator()().

◆ _master

const PointLocatorBase* libMesh::PointLocatorBase::_master
protectedinherited

Const pointer to our master, initialized to nullptr if none given. When using multiple PointLocators, one can be assigned master and be in charge of something that all can have access to.

Definition at line 186 of file point_locator_base.h.

Referenced by clear(), and init().

◆ _mesh

const MeshBase& libMesh::PointLocatorBase::_mesh
protectedinherited

constant reference to the mesh in which the point is looked for.

Definition at line 191 of file point_locator_base.h.

Referenced by libMesh::PointLocatorBase::get_mesh(), init(), perform_fuzzy_linear_search(), and perform_linear_search().

◆ _mutex

Threads::spin_mutex libMesh::ReferenceCounter::_mutex
staticprotectedinherited

Mutual exclusion object to enable thread-safe reference counting.

Definition at line 135 of file reference_counter.h.

◆ _n_objects

Threads::atomic< unsigned int > libMesh::ReferenceCounter::_n_objects
staticprotectedinherited

The number of objects. Print the reference count information when the number returns to 0.

Definition at line 130 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::n_objects(), libMesh::ReferenceCounter::ReferenceCounter(), and libMesh::ReferenceCounter::~ReferenceCounter().

◆ _out_of_mesh_mode

bool libMesh::PointLocatorTree::_out_of_mesh_mode
protected

true if out-of-mesh mode is enabled. See enable_out_of_mesh_mode() for details.

Definition at line 184 of file point_locator_tree.h.

Referenced by disable_out_of_mesh_mode(), enable_out_of_mesh_mode(), and operator()().

◆ _target_bin_size

unsigned int libMesh::PointLocatorTree::_target_bin_size
protected

Target bin size, which gets passed to the constructor of _tree.

Definition at line 189 of file point_locator_tree.h.

Referenced by get_target_bin_size(), and set_target_bin_size().

◆ _tree

TreeBase* libMesh::PointLocatorTree::_tree
protected

Pointer to our tree. The tree is built at run-time through init(). For servant PointLocators (not master), this simply points to the tree of the master.

Definition at line 171 of file point_locator_tree.h.

Referenced by clear(), init(), and operator()().

◆ _use_close_to_point_tol

bool libMesh::PointLocatorBase::_use_close_to_point_tol
protectedinherited

true if we will use a user-specified tolerance for locating the element.

Definition at line 202 of file point_locator_base.h.

Referenced by operator()(), libMesh::PointLocatorBase::set_close_to_point_tol(), and libMesh::PointLocatorBase::unset_close_to_point_tol().

◆ _verbose

bool libMesh::PointLocatorBase::_verbose
inherited

Boolean flag to indicate whether to print out extra info.

Definition at line 178 of file point_locator_base.h.

Referenced by operator()().


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