#include <NearestNodeThread.h>

Public Member Functions
	NearestNodeThread (const MooseMesh &mesh, std::map< dof_id_type, std::vector< dof_id_type >> &neighbor_nodes)

	NearestNodeThread (NearestNodeThread &x, Threads::split split)

void	operator() (const NodeIdRange &range)
	Save a patch of nodes that are close to each of the slave nodes to speed the search algorithm TODO: This needs to be updated at some point in time. More...

void	join (const NearestNodeThread &other)

Public Attributes
std::map< dof_id_type, NearestNodeLocator::NearestNodeInfo >	_nearest_node_info

Real	_max_patch_percentage

Protected Attributes
const MooseMesh &	_mesh

std::map< dof_id_type, std::vector< dof_id_type > > &	_neighbor_nodes

Detailed Description

Definition at line 15 of file NearestNodeThread.h.

Constructor & Destructor Documentation

◆ NearestNodeThread() [1/2]

NearestNodeThread::NearestNodeThread	(	const MooseMesh &	mesh,
		std::map< dof_id_type, std::vector< dof_id_type >> &	neighbor_nodes
	)

Definition at line 18 of file NearestNodeThread.C.

   : _max_patch_percentage(0.0), _mesh(mesh), _neighbor_nodes(neighbor_nodes)
 {
 }

◆ NearestNodeThread() [2/2]

NearestNodeThread::NearestNodeThread	(	NearestNodeThread &	x,
		Threads::split	split
	)

Definition at line 25 of file NearestNodeThread.C.

   : _max_patch_percentage(x._max_patch_percentage),
     _mesh(x._mesh),
     _neighbor_nodes(x._neighbor_nodes)
 {
 }

Member Function Documentation

◆ join()

void NearestNodeThread::join ( const NearestNodeThread & other )

Definition at line 93 of file NearestNodeThread.C.

 {
   // Did the other one go further through the patch than this one?
   if (other._max_patch_percentage > _max_patch_percentage)
     _max_patch_percentage = other._max_patch_percentage;
 
   _nearest_node_info.insert(other._nearest_node_info.begin(), other._nearest_node_info.end());
 }

◆ operator()()

void NearestNodeThread::operator() ( const NodeIdRange & range )

Save a patch of nodes that are close to each of the slave nodes to speed the search algorithm TODO: This needs to be updated at some point in time.

If the hits into this data structure approach "the end" then it may be time to update

Definition at line 39 of file NearestNodeThread.C.

 {
   for (const auto & node_id : range)
   {
     const Node & node = _mesh.nodeRef(node_id);
 
     const Node * closest_node = NULL;
     Real closest_distance = std::numeric_limits<Real>::max();
 
     const std::vector<dof_id_type> & neighbor_nodes = _neighbor_nodes[node_id];
 
     unsigned int n_neighbor_nodes = neighbor_nodes.size();
 
     for (unsigned int k = 0; k < n_neighbor_nodes; k++)
     {
       const Node * cur_node = &_mesh.nodeRef(neighbor_nodes[k]);
       Real distance = ((*cur_node) - node).norm();
 
       if (distance < closest_distance)
       {
         Real patch_percentage = static_cast<Real>(k) / static_cast<Real>(n_neighbor_nodes);
 
         // Save off the maximum we had to go through the patch to find the closes node
         if (patch_percentage > _max_patch_percentage)
           _max_patch_percentage = patch_percentage;
 
         closest_distance = distance;
         closest_node = cur_node;
       }
     }
 
     if (closest_distance == std::numeric_limits<Real>::max())
     {
       for (unsigned int k = 0; k < n_neighbor_nodes; k++)
       {
         const Node * cur_node = &_mesh.nodeRef(neighbor_nodes[k]);
         if (std::isnan((*cur_node)(0)) || std::isinf((*cur_node)(0)) ||
             std::isnan((*cur_node)(1)) || std::isinf((*cur_node)(1)) ||
             std::isnan((*cur_node)(2)) || std::isinf((*cur_node)(2)))
           mooseError("Failure in NearestNodeThread because solution contains inf or not-a-number "
                      "entries.  This is likely due to a failed factorization of the Jacobian "
                      "matrix.");
       }
       mooseError("Unable to find nearest node!");
     }
 
     NearestNodeLocator::NearestNodeInfo & info = _nearest_node_info[node.id()];
 
     info._nearest_node = closest_node;
     info._distance = closest_distance;
   }
 }

Member Data Documentation

◆ _max_patch_percentage

Real NearestNodeThread::_max_patch_percentage

Definition at line 33 of file NearestNodeThread.h.

Referenced by NearestNodeLocator::findNodes(), join(), and operator()().

◆ _mesh

const MooseMesh& NearestNodeThread::_mesh

protected

Definition at line 37 of file NearestNodeThread.h.

Referenced by operator()().

◆ _nearest_node_info

std::map<dof_id_type, NearestNodeLocator::NearestNodeInfo> NearestNodeThread::_nearest_node_info

Definition at line 29 of file NearestNodeThread.h.

Referenced by NearestNodeLocator::findNodes(), join(), and operator()().

◆ _neighbor_nodes

std::map<dof_id_type, std::vector<dof_id_type> >& NearestNodeThread::_neighbor_nodes