Changeset 8489

Timestamp:

Aug 5, 2012, 12:36:26 PM (12 years ago)

Author:

steve

Message:

Added some output to run_parallel_sw_merimbula_profile

Location:

trunk/anuga_core/source/anuga_parallel

Files:

• parallel_generic_communications.py (modified) (4 diffs)
• run_parallel_sw_merimbula_profile.py (modified) (6 diffs)

trunk/anuga_core/source/anuga_parallel/parallel_generic_communications.py

@@ -43 +43 @@
     #Determine minimum timestep across all processors
     pypar.reduce(domain.local_timestep, pypar.MIN, 0,
-                      buffer=domain.global_timestep)#,
-                     #bypass=True)
+                      buffer=domain.global_timestep,
+                      bypass=True)
+
+    #pypar.allreduce(domain.local_timestep, pypar.MIN,
+    #                  buffer=domain.global_timestep,
+    #                  bypass=True)
 
     domain.communication_reduce_time += time.time()-t0
@@ -51 +55 @@
     #Broadcast minimal timestep to all processors
     t0 = time.time()
-    pypar.broadcast(domain.global_timestep, 0)#,
-                        #bypass=True)
+    pypar.broadcast(domain.global_timestep, 0, bypass=True)
 
     domain.communication_broadcast_time += time.time()-t0
 
-    old_timestep = domain.flux_timestep
+    old_fux_timestep = domain.flux_timestep
     domain.flux_timestep = domain.global_timestep[0]
     
     
-    #Compute minimal timestep on local process
-    #Domain.update_timestep(domain, yieldstep, finaltime)
-
 
 def communicate_ghosts(domain):
@@ -90 +90 @@
                         Xout[:,i] = num.take(Q_cv, Idf)
 
-                    pypar.send(Xout, int(send_proc), use_buffer=True)
+                    pypar.send(Xout, int(send_proc), use_buffer=True, bypass=True)
 
 
@@ -100 +100 @@
                 X   = domain.ghost_recv_dict[iproc][2]
 
-                X = pypar.receive(int(iproc), buffer=X)
+                X = pypar.receive(int(iproc), buffer=X, bypass=True)
 
                 for i, q in enumerate(domain.conserved_quantities):

trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula_profile.py

@@ -8 +8 @@
    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.
+   Will produce sww files with names domain_Pn_m.sww where n is number of processors and
+   m in [0, n-1] refers to specific processor that owned this part of the partitioned mesh.
 """
 
@@ -46 +46 @@
 #mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5
 
-finaltime = 50
-yieldstep = finaltime
+finaltime = 500
+yieldstep = 50
 verbose = True
 
@@ -82 +82 @@
 
 #domain.smooth = False
+for id in range(numprocs):
+    if myid == id:
+        print 'Process ID %g' %myid
+        print 'Number of triangles %g ' %domain.get_number_of_triangles()
+
+    barrier()
 
 
@@ -128 +134 @@
 
 
+barrier()
+
 if myid == 0:
     import pstats
     p = pstats.Stats(prof_file)
-    p.sort_stats('cumulative').print_stats(25)
+    #p.sort_stats('cumulative').print_stats(25)
 
 
     p.sort_stats('time').print_stats(25)
 
+barrier()
 
+
+if myid == 1:
+    import pstats
+    p = pstats.Stats(prof_file)
+    #p.sort_stats('cumulative').print_stats(25)
+
+
+    p.sort_stats('time').print_stats(25)
 
     #p.print_stats()
@@ -146 +163 @@
 #profiler.close()
 
-
+barrier()
 
 if myid == 0:
+    print 'Process ID %g' %myid
+    print 'Number of processors %g ' %numprocs
+    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
+
+barrier()
+
+if myid == 1:
+    print 'Process ID %g' %myid
     print 'Number of processors %g ' %numprocs
     print 'That took %.2f seconds' %(time.time()-t0)
@@ -156 +184 @@
 
 
+
 finalize()