Changeset 8856

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

Location:

trunk/anuga_core/source/anuga_parallel

Files:

• distribute_mesh.py (modified) (4 diffs)
• run_parallel_sw_merimbula.py (modified) (2 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
 

trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

@@ -46 +46 @@
 #mesh_filename = "merimbula_17156.tsh"   ; x0 = 756000.0 ; x1 = 756500.0; yieldstep = 50; finaltime = 500
 #mesh_filename = "merimbula_43200_1.tsh"   ; x0 = 756000.0 ; x1 = 756500.0; yieldstep = 50; finaltime = 500
-#mesh_filename = "test-100.tsh" ; x0 = 200.0 ; x1 = 300.0; yieldstep = 1; finaltime = 50
+mesh_filename = "test-100.tsh" ; x0 = 200.0 ; x1 = 300.0; yieldstep = 1; finaltime = 50
 #mesh_filename = "test-20.tsh" ; x0 = 250.0 ; x1 = 350.0; yieldstep = 1; finaltime = 50
 
@@ -172 +172 @@
 t0 = time.time()
 for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
-    if myid == 0:
-        domain.write_time()
+    print 'P:'+ str(myid) + ' '+ domain.timestepping_statistics()
+    #if myid == 0:
+    #    domain.write_time()