Changeset 4237

Timestamp:

Feb 7, 2007, 4:38:38 PM (18 years ago)

Author:

duncan

Message:

adding profiling

File:

• anuga_validation/automated_validation_tests/okushiri_tank_validation/loading_pts_test.py (modified) (1 diff)

anuga_validation/automated_validation_tests/okushiri_tank_validation/loading_pts_test.py

@@ -49 +49 @@
 bathymetry_filename=project.bathymetry_filename
 print 'Starting domain.set_quantity.  Loading ', bathymetry_filename
-domain.set_quantity('elevation',
-                    filename=bathymetry_filename,
-                    alpha=0.02,                    
-                    verbose=True,
-                    use_cache=False)
+s = "domain.set_quantity('elevation',filename=bathymetry_filename,alpha=0.02,verbose=True,use_cache=False)"
+
+
+import profile, pstats
+FN = 'profile.dat'
+
+profile.run(s, FN)
 
 print 'Set_quantity elevation took %.2f seconds' %(time.time()-t0)
-import sys; sys.exit() 
 
-#-------------------------
-# Distribute domain if run in parallel
-#-------------------------
-if numprocs > 1:
-    domain = distribute(domain)
+S = pstats.Stats(FN)
+#S.sort_stats('time').print_stats(20)
+s = S.sort_stats('cumulative').print_stats(30)
 
-
-#-------------------------
-# Set simulation parameters
-#-------------------------
-domain.set_name(project.output_filename)  # Name of output sww file
-domain.set_default_order(2)               # Apply second order scheme 
-domain.set_all_limiters(0.9)              # Max second order scheme (old lim)
-domain.set_minimum_storable_height(0.001) # Don't store w < 0.001m
-domain.set_maximum_allowed_speed(0.1)     # Allow a little runoff (0.1 is OK)
-
-
-#-------------------------
-# Boundary Conditions
-#-------------------------
-
-# Create boundary function from timeseries provided in file
-function = file_function(project.boundary_filename,
-                         domain, verbose=True)
-
-# Create and assign boundary objects
-Bts = Transmissive_Momentum_Set_Stage_boundary(domain, function)
-Br = Reflective_boundary(domain)
-domain.set_boundary({'wave': Bts, 'wall': Br})
-
-
-#-------------------------
-# Evolve through time
-#-------------------------
-import time
-t0 = time.time()
-
-for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
-    domain.write_time()
-
-print 'That took %.2f seconds' %(time.time()-t0)
+print s