diff --git a/meshmode/mesh/__init__.py b/meshmode/mesh/__init__.py index cc32e962cd785693c9e0d064276d4b550577fedd..472e6c4f1d2f27bce7f8b9e541c480f216d27b29 100644 --- a/meshmode/mesh/__init__.py +++ b/meshmode/mesh/__init__.py @@ -872,7 +872,7 @@ class Mesh(Record): if elem < grp.nelements: return igrp elem -= grp.nelements - raise RuntimeError("Could not find group with element ", elem) + raise RuntimeError("Could not find group with element %d" % elem) # Design experience: Try not to add too many global data structures to the # mesh. Let the element groups be responsible for that at the mesh level. diff --git a/meshmode/mesh/processing.py b/meshmode/mesh/processing.py index 19d8e9f788fc84ce32dace56517c46ba03ca51e4..c028c63fbb8491de6c5dfda1fbb51745d8aae177 100644 --- a/meshmode/mesh/processing.py +++ b/meshmode/mesh/processing.py @@ -185,6 +185,8 @@ def partition_mesh(mesh, part_per_element, part_nr): n_elem = np.count_nonzero( part_per_element[:rank_neighbor] == n_part_num) + # I cannot compute the group because the other + # partitions have not been built yet. adj_grps[igrp].add_connection( elem, face,