Changeset 8497

Timestamp:

Aug 7, 2012, 12:22:48 PM (13 years ago)

Author:

steve

Message:

Move over to opl machine

Location:

trunk/anuga_core/source/anuga_parallel

Files:

• parallel_generic_communications.py (modified) (4 diffs)
• run_parallel_sw_merimbula.py (modified) (1 diff)
• run_parallel_sw_rectangular_cross_profile.py (modified) (4 diffs)

trunk/anuga_core/source/anuga_parallel/parallel_generic_communications.py

@@ -41 +41 @@
     t0 = time.time()
 
-    allreduce = False
+    allreduce = True
     
     if allreduce:
-        #pypar.allreduce(domain.local_timestep, pypar.MIN,
-        #              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
     else:
@@ -62 +62 @@
         domain.communication_broadcast_time += time.time()-t0
 
-    old_fux_timestep = domain.flux_timestep
+    #old_flux_timestep = domain.flux_timestep
     domain.flux_timestep = domain.global_timestep[0]
     
@@ -93 +93 @@
                         Xout[:,i] = num.take(Q_cv, Idf)
 
-                    pypar.send(Xout, int(send_proc), use_buffer=True)#, bypass=True)
+                    pypar.send(Xout, int(send_proc), use_buffer=True, bypass=True)
 
 
@@ -103 +103 @@
                 X   = domain.ghost_recv_dict[iproc][2]
 
-                X = pypar.receive(int(iproc), buffer=X)#, bypass=True)
+                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.py

@@ -109 +109 @@
 #--------------------------------------------------------------------------
 
-domain.set_flow_algorithm('2_0')   
-
+domain.set_flow_algorithm('2_0')
 
 #domain.smooth = False

trunk/anuga_core/source/anuga_parallel/run_parallel_sw_rectangular_cross_profile.py

@@ -46 +46 @@
 #mesh_filename = "test-100.tsh" ; x0 = 0.25 ; x1 = 0.5
 
-finaltime = 0.5
-yieldstep = 0.05
+
 verbose = True
 
@@ -69 +68 @@
 #--------------------------------------------------------------------------
 if myid == 0:
-    length = 2.0*numprocs
+    length = 2.0
     width = 2.0
-    dx = dy = 0.04
+    dx = dy = 0.02
     domain = rectangular_cross_domain(int(length/dx), int(width/dy),
                                               len1=length, len2=width)
@@ -115 +114 @@
 Br = Reflective_boundary(domain)      # Solid reflective wall
 
-domain.set_boundary({'top' :Br, 'bottom' :Br, 'left' :Br, 'right' :Br })
+domain.set_boundary({'top' :Br, 'bottom' :Br, 'left' :Br, 'right' :Br, 'ghost' : None })
 
 """
@@ -133 +132 @@
 
 t0 = time.time()
-
+finaltime = 2.0
+yieldstep = 0.5
 
 s = """