Changeset 8856 for trunk/anuga_core/source/anuga_parallel/distribute_mesh.py

Timestamp:

Apr 30, 2013, 8:42:00 PM (12 years ago)

Author:

steve

Message:

Got an error when distributing a small mesh over 17 processors

File:

• trunk/anuga_core/source/anuga_parallel/distribute_mesh.py (modified) (4 diffs)

trunk/anuga_core/source/anuga_parallel/distribute_mesh.py

@@ -132 +132 @@
     # Wrapper for old pmesh_divide_metis which does not return tri_index or r_tri_index
     nodes, ttriangles, boundary, triangles_per_proc, quantities, tri_index, r_tri_index = pmesh_divide_metis_helper(domain, n_procs)
+
+    print triangles_per_proc
+    assert num.all(triangles_per_proc>0), 'Metis created a partition where one mesh has 0 triangles'
+
     return nodes, ttriangles, boundary, triangles_per_proc, quantities
 
 def pmesh_divide_metis_with_map(domain, n_procs):
+
     return pmesh_divide_metis_helper(domain, n_procs)
 
@@ -208 +213 @@
 
         triangles_per_proc = num.bincount(epart)
+
+        msg = """Metis created a partition where at least one submesh has no triangles.
+        Try using a smaller number of mpi processors"""
+        assert num.all(triangles_per_proc>0), msg
+
         proc_sum = num.zeros(n_procs+1,num.int)
         proc_sum[1:] = num.cumsum(triangles_per_proc)
@@ -400 +410 @@
 
 
+    print triangles_per_proc
+    
     nodes = mesh.nodes
     triangles = mesh.triangles
@@ -450 +462 @@
 
         ids = num.unique(subtriangles.flat)
+
+            
         lnodes = nodes[ids]
-#        print nodes.shape
-#        print ids.shape
-#        print lnodes.shape
+#       print nodes.shape
+#       print ids.shape
+#       print lnodes.shape
         x = num.concatenate((num.reshape(ids, (-1,1)),lnodes ), 1)
-#        print x
-#        print y
+#       print x
+#       print y
         node_list.append(x)
+
+            
         # Move to the next processor