Also allow remote_elems to neighbor active local 1D elements

libMesh · May 3, 2024 · efdaf92 · efdaf92
1 parent 0144e6d
commit efdaf92
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/src/mesh/unstructured_mesh.C b/src/mesh/unstructured_mesh.C
@@ -1068,11 +1068,13 @@ void UnstructuredMesh::find_neighbors (const bool reset_remote_elements,
                       libmesh_assert(neigh_has_remote_children ||
                                      neigh->dim() == 1);
 
-                      // And let's double-check that we don't have
-                      // a remote_elem neighboring an active local element
+                      // And let's double-check that we don't have a
+                      // remote_elem neighboring an active local element
+                      // unless this is a 1D Elem, for the reasons
+                      // discussed above.
                       if (current_elem->active())
-                        libmesh_assert_not_equal_to (current_elem->processor_id(),
-                                                     this->processor_id());
+                        libmesh_assert (current_elem->dim() == 1 ||
+                                        current_elem->processor_id() != this->processor_id());
 #endif // DEBUG
                       neigh = const_cast<RemoteElem *>(remote_elem);
                     }