[Libmesh-users] Elem::active_family_tree_by_side

[Manav Bhatia]

Hi, 

      I am curious if SerialMesh with AMR uses the RemoteElem.  In the 
following function (from elem.C), I have marked the lines of interest with 
“****”. 

      This method gets called by line 2180 of nemesis_io_helper.C, which wants 
to find active elements with a common side, but also on the local processor. 
But, the method in elem.C does not seem to make a distinction about elements 
that may have a different pid assigned, but may not be remote elems. 
Elem::is_remote() is false by default. 

      Eventually, error on line 2211 of nemesis_io_helper.C is thrown. 
  
      Maybe it would make sense to add a check for (processor_id != 
this->processor_id) in Elem::is_remote()? Any thoughts? 

Thanks,
Manav


void Elem::active_family_tree_by_side (std::vector<const Elem*>& family,
                                       const unsigned int s,
                                       const bool reset) const
{
  // The "family tree" doesn't include subactive elements
  libmesh_assert(!this->subactive());

  // Clear the vector if the flag reset tells us to.
  if (reset)
    family.clear();

  libmesh_assert_less (s, this->n_sides());

  // Add an active element to the family tree.
  if (this->active())
    family.push_back(this);

  // Or recurse into an ancestor element's children.
  // Do not clear the vector any more.
  else
    for (unsigned int c=0; c<this->n_children(); c++)
  **** if (!this->child(c)->is_remote() && this->is_child_on_side(c, s))  *****
        this->child(c)->active_family_tree_by_side (family, s, false);
}

------------------------------------------------------------------------------
Is your legacy SCM system holding you back? Join Perforce May 7 to find out:
&#149; 3 signs your SCM is hindering your productivity
&#149; Requirements for releasing software faster
&#149; Expert tips and advice for migrating your SCM now
http://p.sf.net/sfu/perforce
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users