Changeset 4725

Timestamp:

Sep 11, 2007, 12:02:31 PM (17 years ago)

Author:

ole

Message:

Investigated mesh (and timestepping) differences across platforms

Location:

anuga_work/development/demos

Files:

• island.py (modified) (2 diffs)
• island_linux.msh (added)
• island_windows.msh (added)

anuga_work/development/demos/island.py

@@ -30 +30 @@
 #------------------------------------------------------------------------------
 
-#Create basic mesh
-create_mesh_from_regions( [[0,0], [100,0], [100,100], [0,100]],
-                          boundary_tags = {'bottom': [0],
-                                           'right': [1],
-                                           'top': [2],
-                                           'left': [3]},
-                          maximum_triangle_area = 2,
-                          filename = 'island.msh' ,
-                          interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 10.0)]
-                          #interior_holes=[[[50,25], [70,25], [70,75], [50,75]]],
-                          )
+if False:
+    #Create basic mesh
+    create_mesh_from_regions( [[0,0], [100,0], [100,100], [0,100]],
+                              boundary_tags = {'bottom': [0],
+                                               'right': [1],
+                                               'top': [2],
+                                               'left': [3]},
+                              maximum_triangle_area = 2,
+                              filename = 'island.msh' ,
+                              interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 10.0)]
+                              #interior_holes=[[[50,25], [70,25], [70,75], [50,75]]],
+                              )
 
 
 
+
+
+# Select precomputed meshes here. In changeset:4725 there is a huge difference
+# in timestepping characteristics depending on the mesh used.
+# The mesh generated on windows takes 10 times as many timesteps.
+
+name = 'island_linux'    # Mesh generated on a Linux box (here nautilus) but similar on cyclone
+#name = 'island_windows' # Mesh generated on a Windows box (asset 0031136)
+
+
 #Create shallow water domain
-domain = Domain(mesh_filename = 'island.msh')
+domain = Domain(mesh_filename = name + '.msh')
 domain.smooth = False #True
-domain.set_name('island_unique')
+domain.set_name(name)
 #domain.set_name('island_not_unique')
 domain.default_order = 2
 
+print domain.statistics()
 
 #I tried to introduce this parameter top control the h-limiter,
@@ -89 +101 @@
 domain.maximum_allowed_speed = 0 
 domain.beta_h = 0.0
-#domain.H0 = 0.01
 domain.tight_slope_limiters = 1
-
-