Changeset 8499

Timestamp:

Aug 7, 2012, 4:46:14 PM (13 years ago)

Author:

steve

Message:

Playing with extra layer of ghost cells to reduce communication

File:

• trunk/anuga_core/source/anuga_parallel/distribute_mesh.py (modified) (2 diffs)

trunk/anuga_core/source/anuga_parallel/distribute_mesh.py

@@ -355 +355 @@
     # Find the first layer of boundary triangles
 
+    layer_width = 2
+
+
     trianglemap = num.zeros(ntriangles, 'i')
+
+
     for t in range(tlower, tupper):
         
@@ -372 +377 @@
                 trianglemap[n] = 1
 
-    # Find the second layer of boundary triangles
-
-    for t in range(len(trianglemap)):
-        if trianglemap[t]==1:
-            n = mesh.neighbours[t, 0]
-            if n >= 0:
-                if (n < tlower or n >= tupper) and trianglemap[n] == 0:
-                    trianglemap[n] = 2
-            n = mesh.neighbours[t, 1]
-            if n >= 0:
-                if (n < tlower or n >= tupper) and trianglemap[n] == 0:
-                    trianglemap[n] = 2
-            n = mesh.neighbours[t, 2]
-            if n >= 0:
-                if (n < tlower or n >= tupper) and trianglemap[n] == 0:
-                    trianglemap[n] = 2
+
+    # Find the subsequent layers of ghost triangles
+    for i in range(layer_width-1):
+
+        for t in range(len(trianglemap)):
+            if trianglemap[t]==i+1:
+                n = mesh.neighbours[t, 0]
+                if n >= 0:
+                    if (n < tlower or n >= tupper) and trianglemap[n] == 0:
+                        trianglemap[n] = i+2
+                n = mesh.neighbours[t, 1]
+                if n >= 0:
+                    if (n < tlower or n >= tupper) and trianglemap[n] == 0:
+                        trianglemap[n] = i+2
+                n = mesh.neighbours[t, 2]
+                if n >= 0:
+                    if (n < tlower or n >= tupper) and trianglemap[n] == 0:
+                        trianglemap[n] = i+2
 
     # Build the triangle list and make note of the vertices
+
+
+
+
 
     nodemap = num.zeros(ncoord, 'i')