14 #include "libmesh/point_locator_base.h" 15 #include "libmesh/elem.h" 16 #include "libmesh/enum_point_locator_type.h" 17 #include "libmesh/point.h" 33 std::pair<std::vector<Point>, std::vector<unsigned int>> & multi_point_list =
_points[elem];
35 const unsigned int npoint = multi_point_list.first.size();
36 mooseAssert(npoint == multi_point_list.second.size(),
37 "Different sizes for location and multiplicity data");
39 for (
unsigned int i = 0; i < npoint; ++i)
43 multi_point_list.second[i]++;
48 multi_point_list.first.push_back(p);
49 multi_point_list.second.push_back(1);
62 std::vector<Point> & point_list =
_points[elem].first;
64 for (
const auto & pt : point_list)
85 unsigned pl_needs_rebuild =
_elements.size();
115 const std::set<SubdomainID> &
blocks)
129 mooseError(
"Error, PointLocator is not initialized!");
136 const Elem * elem = (*_point_locator)(p, &
blocks);
142 dof_id_type elem_id = elem ? elem->
id() : DofObject::invalid_id;
149 return min_elem_id == elem_id ? elem : NULL;
bool hasPoint(const Elem *elem, const Point &p)
Return true if we have Point 'p' in Element 'elem'.
const Elem * findPoint(const Point &p, const MooseMesh &mesh, const std::set< SubdomainID > &blocks)
Used by client DiracKernel classes to determine the Elem in which the Point p resides.
std::set< const Elem * > _elements
The list of elements that need distributions.
void mooseError(Args &&... args)
Emit an error message with the given stringified, concatenated args and terminate the application...
static constexpr Real TOLERANCE
const Parallel::Communicator & comm() const
MultiPointMap _points
The list of physical xyz Points that need to be evaluated in each element.
The following methods are specializations for using the libMesh::Parallel::packed_range_* routines fo...
void addPoint(const Elem *elem, const Point &p)
Adds a point source.
std::unique_ptr< libMesh::PointLocatorBase > _point_locator
The DiracKernelInfo object manages a PointLocator object which is used by all DiracKernels to find Po...
void min(const T &r, T &o, Request &req) const
MooseMesh wraps a libMesh::Mesh object and enhances its capabilities by caching additional data and s...
void clearPoints()
Remove all of the current points and elements.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
void max(const T &r, T &o, Request &req) const
const Real _point_equal_distance_sq
threshold distance squared below which two points are considered identical
void updatePointLocator(const MooseMesh &mesh)
Called during FEProblemBase::meshChanged() to update the PointLocator object used by the DiracKernels...
virtual ~DiracKernelInfo()
bool pointsFuzzyEqual(const Point &, const Point &)
Check if two points are equal with respect to a tolerance.