libMesh
sibling_coupling.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2017 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 // Local Includes
21 #include "libmesh/sibling_coupling.h"
22 #include "libmesh/elem.h"
23 #include "libmesh/remote_elem.h"
24 
25 namespace libMesh
26 {
27 
28 void SiblingCoupling::operator()
29  (const MeshBase::const_element_iterator & range_begin,
30  const MeshBase::const_element_iterator & range_end,
32  map_type & coupled_elements)
33 {
34  LOG_SCOPE("operator()", "SiblingCoupling");
35 
36  for (MeshBase::const_element_iterator elem_it = range_begin;
37  elem_it != range_end; ++elem_it)
38  {
39  std::vector<const Elem *> active_siblings;
40 
41  const Elem * const elem = *elem_it;
42 
43  const Elem * parent = elem->parent();
44  if (!parent)
45  continue;
46 
47 #ifdef LIBMESH_ENABLE_AMR
48  parent->active_family_tree(active_siblings);
49 #endif
50 
51  for (std::size_t a=0; a != active_siblings.size(); ++a)
52  {
53  const Elem * sibling = active_siblings[a];
54 
55  if (sibling->processor_id() != p)
56  coupled_elements.insert
57  (std::make_pair(sibling, _dof_coupling));
58  }
59  }
60 }
61 
62 
63 } // namespace libMesh
The definition of the const_element_iterator struct.
Definition: mesh_base.h:1494
const Elem * parent() const
Definition: elem.h:2346
This is the base class from which all geometric element types are derived.
Definition: elem.h:89
uint8_t processor_id_type
Definition: id_types.h:99
The libMesh namespace provides an interface to certain functionality in the library.
void active_family_tree(std::vector< const Elem * > &active_family, const bool reset=true) const
Same as the family_tree() member, but only adds the active children.
Definition: elem.C:1724
std::unordered_map< const Elem *, const CouplingMatrix * > map_type
What elements do we care about and what variables do we care about on each element?
processor_id_type processor_id() const
Definition: dof_object.h:694