Changeset 7676

Timestamp:

Apr 6, 2010, 8:46:17 PM (15 years ago)

Author:

steve

Message:

 

File:

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

anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

@@ -1 +1 @@
 
 """Run parallel shallow water domain.
+
+   run using command like:
+
+   mpirun -np m python run_parallel_sw_merimbula.py
+
+   where m is the number of processors to be used.
+   
+   Will produce sww files with names domain_Pn_m.sww where m is number of processors and
+   n in [0, m-1] refers to specific processor that owned this part of the partitioned mesh.
 """
 
@@ -10 +19 @@
 import os
 import sys
+import time
 import pypar
 
@@ -38 +48 @@
 #--------------------------------------------------------------------------
 
-mesh_filename = "merimbula_10785_1.tsh"
+#mesh_filename = "merimbula_10785_1.tsh"
+mesh_filename = "merimbula_43200.tsh"
 #mesh_filename = "test-100.tsh"
-yieldstep = 1
-finaltime = 20
+yieldstep = 50
+finaltime = 200
 quantity = 'stage'
 verbose = True
@@ -88 +99 @@
 #------------------------------------------------------------------------------
 if myid == 0 and verbose: print 'EVOLVE'
-        
+
+t0 = time.time()
+
 for t in domain.evolve(yieldstep = yieldstep, finaltime = finaltime):
     if myid == 0:
         domain.write_time()
+
+
+if myid == 0:
+    print 'That took %.2f seconds' %(time.time()-t0)
+    print 'Communication time %.2f seconds'%domain.communication_time
+    print 'Reduction Communication time %.2f seconds'%domain.communication_reduce_time
+    print 'Broadcast time %.2f seconds'%domain.communication_broadcast_time
+
 
 finalize()