libMesh::GhostPointNeighbors Class Reference

#include <ghost_point_neighbors.h>

Inheritance diagram for libMesh::GhostPointNeighbors:

Public Types
typedef std::unordered_map< const Elem *, const CouplingMatrix * >	map_type

Public Member Functions
	GhostPointNeighbors (const MeshBase &mesh)
virtual void	operator() (const MeshBase::const_element_iterator &range_begin, const MeshBase::const_element_iterator &range_end, processor_id_type p, map_type &coupled_elements)
virtual void	mesh_reinit ()
virtual void	dofmap_reinit ()
virtual void	redistribute ()
virtual void	delete_remote_elements ()

Static Public Member Functions
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 ()

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)

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

Private Attributes
const MeshBase &	_mesh

Detailed Description

This class implements the default geometry ghosting in libMesh: point neighbors and interior_parent elements are ghosted.

Author

Roy H. Stogner

Date

2016

Definition at line 36 of file ghost_point_neighbors.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.

map_type

typedef std::unordered_map<const Elem*, const CouplingMatrix*> libMesh::GhostingFunctor::map_type

inherited

What elements do we care about and what variables do we care about on each element?

Definition at line 154 of file ghosting_functor.h.

Constructor & Destructor Documentation

GhostPointNeighbors()

libMesh::GhostPointNeighbors::GhostPointNeighbors ( const MeshBase &  mesh )

inline

Constructor.

Definition at line 43 of file ghost_point_neighbors.h.

: _mesh(mesh) {}

Member Function Documentation

delete_remote_elements()

virtual void libMesh::GhostingFunctor::delete_remote_elements ( )

inlinevirtualinherited

GhostingFunctor subclasses with relatively long-lasting caches may want to delete the no-longer-relevant parts of those caches after a redistribution is complete.

Reimplemented in libMesh::PointNeighborCoupling, and libMesh::DefaultCoupling.

Definition at line 204 of file ghosting_functor.h.

{};

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().

{
  _enable_print_counter = false;
  return;
}

dofmap_reinit()

virtual void libMesh::GhostingFunctor::dofmap_reinit ( )

inlinevirtualinherited

For algebraic ghosting or coupling functors we also call dofmap_reinit() later, after dofs have been distributed on the new mesh but before the functors have been queried for send_list or sparsity pattern calculations.

Definition at line 185 of file ghosting_functor.h.

{};

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.

{
  _enable_print_counter = true;
  return;
}

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().

{
#if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)

  std::ostringstream oss;

  oss << '\n'
      << " ---------------------------------------------------------------------------- \n"
      << "| Reference count information                                                |\n"
      << " ---------------------------------------------------------------------------- \n";

  for (const auto & pr : _counts)
    {
      const std::string name(pr.first);
      const unsigned int creations    = pr.second.first;
      const unsigned int destructions = pr.second.second;

      oss << "| " << name << " reference count information:\n"
          << "|  Creations:    " << creations    << '\n'
          << "|  Destructions: " << destructions << '\n';
    }

  oss << " ---------------------------------------------------------------------------- \n";

  return oss.str();

#else

  return "";

#endif
}

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().

{
  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
  std::pair<unsigned int, unsigned int> & p = _counts[name];

  p.first++;
}

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().

{
  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
  std::pair<unsigned int, unsigned int> & p = _counts[name];

  p.second++;
}

mesh_reinit()

virtual void libMesh::GhostingFunctor::mesh_reinit ( )

inlinevirtualinherited

GhostingFunctor subclasses which cache data will need to initialize that cache. We call mesh_reinit() whenever the relevant Mesh has changed, but before remote elements on a distributed mesh are deleted.

Reimplemented in libMesh::PointNeighborCoupling, and libMesh::DefaultCoupling.

Definition at line 177 of file ghosting_functor.h.

{};

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.

  { return _n_objects; }

operator()()

void libMesh::GhostPointNeighbors::operator() ( const MeshBase::const_element_iterator &  range_begin, const MeshBase::const_element_iterator &  range_end, processor_id_type  p, GhostPointNeighbors::map_type &  coupled_elements  )

virtual

For the specified range of active elements, find their point neighbors and interior_parent elements, ignoring those on processor p.

Implements libMesh::GhostingFunctor.

Definition at line 33 of file ghost_point_neighbors.C.

References libMesh::as_range(), and libMesh::remote_elem.

{
  // Using the connected_nodes set rather than point_neighbors() gives
  // us correct results even in corner cases, such as where two
  // elements meet only at a corner.  ;-)

  std::unordered_set<const Node *> connected_nodes;

  // Links between boundary and interior elements on mixed
  // dimensional meshes also give us correct ghosting in this way.
  std::unordered_set<const Elem *> interior_parents;

  // We also preserve neighbors and their neighboring children for
  // active local elements - in most cases this is redundant with the
  // node check, but for non-conforming Tet4 meshes and
  // non-level-one-conforming 2D+3D meshes it is possible for an
  // element and its coarse neighbor to not share any vertices.
  //
  // We also preserve interior parents for active pid elements


  // This code is just for geometric coupling, so we use a null
  // CouplingMatrix pointer.  We'll declare that here so as to avoid
  // confusing the insert() calls later.
  CouplingMatrix * nullcm = nullptr;

  for (const auto & elem : as_range(range_begin, range_end))
    {
      if (elem->processor_id() != p)
        coupled_elements.insert (std::make_pair(elem,nullcm));

      for (auto neigh : elem->neighbor_ptr_range())
        {
          if (neigh && neigh != remote_elem)
            {
#ifdef LIBMESH_ENABLE_AMR
              if (!neigh->active())
                {
                  std::vector<const Elem*> family;
                  neigh->active_family_tree_by_neighbor(family, elem);

                  for (const Elem * f : family)
                    if (f->processor_id() != p)
                      coupled_elements.insert
                        (std::make_pair(f, nullcm));
                }
              else
#endif
                if (neigh->processor_id() != p)
                  coupled_elements.insert
                    (std::make_pair(neigh, nullcm));
            }
        }

      // It is possible that a refined boundary element will not
      // touch any nodes of its interior_parent, in TRI3/TET4 and in
      // non-level-one rule cases.  So we can't just rely on node
      // connections to preserve interior_parent().  However, trying
      // to preserve interior_parent() manually only works if it's on
      // the same Mesh, which is *not* guaranteed!  So we'll
      // double-check later to make sure our interior parents are in
      // the mesh before we connect them.
      if (elem->dim() < LIBMESH_DIM &&
          elem->interior_parent() &&
          elem->interior_parent()->processor_id() != p)
        interior_parents.insert (elem->interior_parent());

      // Add nodes connected to active local elements
      for (auto n : elem->node_index_range())
        connected_nodes.insert (elem->node_ptr(n));
    }

  // Connect any interior_parents who are really in our mesh
  for (const auto & elem : _mesh.element_ptr_range())
    {
      std::unordered_set<const Elem *>::iterator ip_it =
        interior_parents.find(elem);

      if (ip_it != interior_parents.end())
        {
          coupled_elements.insert
            (std::make_pair(elem, nullcm));

          // Shrink the set ASAP to speed up subsequent searches
          interior_parents.erase(ip_it);
        }
    }

  // Connect any active elements which are connected to our range's
  // elements' nodes by addin elements connected to nodes on active
  // local elements.
  for (const auto & elem : _mesh.active_element_ptr_range())
    if (elem->processor_id() != p)
      for (auto & n : elem->node_ref_range())
        if (connected_nodes.count(&n))
          coupled_elements.insert(std::make_pair(elem, nullcm));
}

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().

{
  if (_enable_print_counter)
    out_stream << ReferenceCounter::get_info();
}

redistribute()

virtual void libMesh::GhostingFunctor::redistribute ( )

inlinevirtualinherited

GhostingFunctor subclasses with relatively long-lasting caches may want to redistribute those caches whenever the relevant Mesh is redistributed; we will give them an opportunity when that happens. At the point in the code where this is called, element processor ids have been set to their new destinations, and those elements have been copied to their new destinations, but the elements have not yet been deleted by the processors which previously held them..

Reimplemented in libMesh::PointNeighborCoupling, and libMesh::DefaultCoupling.

Definition at line 197 of file ghosting_functor.h.

{};

Member Data Documentation

_counts

ReferenceCounter::Counts libMesh::ReferenceCounter::_counts

staticprotectedinherited

Actually holds the data.

Definition at line 122 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::get_info(), libMesh::ReferenceCounter::increment_constructor_count(), and libMesh::ReferenceCounter::increment_destructor_count().

_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().

_mesh

const MeshBase& libMesh::GhostPointNeighbors::_mesh

private

Definition at line 57 of file ghost_point_neighbors.h.

_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().

---

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

• ghost_point_neighbors.h
• ghost_point_neighbors.C