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: • 3 signs your SCM is hindering your productivity • Requirements for releasing software faster • 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