Changeset 1582

Timestamp:

Jul 5, 2005, 4:48:45 PM (19 years ago)

Author:

linda

Message:

fixed bug in build_commun for receiving quantities

Location:

inundation/ga/storm_surge/parallel

Files:

• build_commun.py (modified) (2 diffs)
• run_parallel_sw_merimbula.py (modified) (2 diffs)

inundation/ga/storm_surge/parallel/build_commun.py

@@ -116 +116 @@
     for k in submesh["full_quan"]:
         pypar.send(submesh["full_quan"][k][p], p, use_buffer=True)
+        
     for k in submesh["ghost_quan"]:
         pypar.send(submesh["ghost_quan"][k][p], p, use_buffer=True)
@@ -252 +253 @@
     new_quan = zeros((no_full_triangles, 3), Float)
     submesh_cell["full_quan"]={}
+    
     for i in range(no_quantities):
-        submesh_cell["full_quan"][qkeys[i]] = pypar.receive(p, new_quan)
+        tmp = pypar.receive(p, new_quan)
+        submesh_cell["full_quan"][qkeys[i]]=zeros((no_full_triangles,3), Float)
+        submesh_cell["full_quan"][qkeys[i]][:] = tmp[:]
 
     new_quan = zeros((no_ghost_triangles, 3), Float)
     submesh_cell["ghost_quan"]={}
     for i in range(no_quantities):
-        submesh_cell["ghost_quan"][qkeys[i]] = pypar.receive(p, new_quan)
+        tmp = pypar.receive(p, new_quan)
+        submesh_cell["ghost_quan"][qkeys[i]]= zeros((no_ghost_triangles,3), Float)
+        submesh_cell["ghost_quan"][qkeys[i]][:] = tmp[:]
     
     # find the full triangles assigned to this processor

inundation/ga/storm_surge/parallel/run_parallel_sw_merimbula.py

@@ -88 +88 @@
 
 
-
-
 if myid == 0:
 
@@ -108 +106 @@
     nodes, triangles, boundary, triangles_per_proc, quantities = \
          pmesh_divide(domain_full, nx, ny)
-
+    
     rect = array(domain_full.xy_extent, Float)
 
     submesh = build_submesh(nodes, triangles, boundary,\
                             quantities, triangles_per_proc)
-
+    
     # send the mesh partition to the appropriate processor