Changeset 5847 for anuga_validation

Timestamp:

Oct 18, 2008, 9:13:18 PM (16 years ago)

Author:

steve

Message:

Changed parallel_api so that global mesh only needs to
be constructed on processor 0

Location:

anuga_validation

Files:

• okushiri_2005/run_okushiri_parallel.py (modified) (3 diffs)
• performance_tests/run_profile.py (modified) (1 diff)

anuga_validation/okushiri_2005/run_okushiri_parallel.py

@@ -39 +39 @@
 #-------------------------
 # Create Domain from mesh
+# on processor 0
 #-------------------------
-domain = Domain(project.mesh_filename, use_cache=use_cache, verbose=True)
-print domain.statistics()
+if myid == 0 :
+    domain = Domain(project.mesh_filename, use_cache=use_cache, verbose=True)
+    print domain.statistics()
 
-
-#-------------------------
-# Initial Conditions
-#-------------------------
-domain.set_quantity('friction', 0.0)
-domain.set_quantity('stage', 0.0)
-domain.set_quantity('elevation',
-                    filename=project.bathymetry_filename,
-                    alpha=0.02,                    
-                    verbose=True,
-                    use_cache=use_cache)
-
+    #-------------------------
+    # Initial Conditions
+    # At present need to do the
+    # fitting of the bathymetry
+    # on a global domain, ie before
+    # distributing the domain
+    #-------------------------
+    domain.set_quantity('friction', 0.0)
+    domain.set_quantity('stage', 0.0)
+    domain.set_quantity('elevation',
+                        filename=project.bathymetry_filename,
+                        alpha=0.02,                    
+                        verbose=True,
+                        use_cache=use_cache)
+else:
+    domain = None
 
 #-------------------------
@@ -86 +92 @@
 domain.set_boundary({'wave': Bts, 'wall': Br})
 
+# Select triangle containing ch5 for diagnostic output
+# around known gauge
+# This should get triangle id 32833 with centroid (4.5244, 1.1972) on one processor
+try: 
+    triangle_id = domain.get_triangle_containing_point([4.521, 1.196])
+except:
+    triangle_id = None
+
 
 #-------------------------
@@ -94 +108 @@
 
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
-    domain.write_time()
+    if myid == 0 :
+        domain.write_time()
+
+    if triangle_id is not None :
+        print domain.timestepping_statistics(track_speeds=False,
+                                             triangle_id=triangle_id)
 
 print 'That took %.2f seconds' %(time.time()-t0)

anuga_validation/performance_tests/run_profile.py

@@ -57 +57 @@
 s = 'for t in domain.evolve(yieldstep = 0.02, finaltime = 0.2): domain.write_time()'
 
+#for t in domain.evolve(yieldstep = 0.02, finaltime = 0.2): domain.write_time()
+
 
 import profile, pstats