#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 14 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 94 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)))
           throw MooseException(
               "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 32 of file NearestNodeThread.h.

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

◆ _mesh

const MooseMesh& NearestNodeThread::_mesh

protected

Definition at line 36 of file NearestNodeThread.h.

Referenced by operator()().

◆ _nearest_node_info

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

Definition at line 28 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