Changeset 6160 for anuga_validation/automated_validation_tests/okushiri_tank_validation

Timestamp:

Jan 14, 2009, 11:17:06 AM (16 years ago)

Author:

rwilson

Message:

Change Numeric imports to general form - ready to change to NumPy?.

Location:

anuga_validation/automated_validation_tests/okushiri_tank_validation

Files:

• compare_timeseries_with_measures.py (modified) (5 diffs)
• create_okushiri.py (modified) (4 diffs)

anuga_validation/automated_validation_tests/okushiri_tank_validation/compare_timeseries_with_measures.py

@@ -6 +6 @@
 """
 
-from Numeric import allclose, argmin, argmax
+import Numeric as num
 from Scientific.IO.NetCDF import NetCDFFile
 
@@ -113 +113 @@
 assert reference_time[0] == 0.0
 assert reference_time[-1] == finaltime
-assert allclose(reference_time, input_time)
+assert num.allclose(reference_time, input_time)
 
 for key in gauge_names:
@@ -152 +152 @@
     #                            rtol=rtol, atol=atol)
     if plotting is False:
-        assert allclose(reference_value, computed_value,
-                        rtol=rtol, atol=atol), msg
+        assert num.allclose(reference_value, computed_value,
+                            rtol=rtol, atol=atol), msg
     
 
@@ -198 +198 @@
     # Locations of extrema
     #i0 = argmax(observed_timeseries)
-    i1 = argmax(model)
+    i1 = num.argmax(model)
     res = reference_time[i1]
     report_difference('Location of maximum', res, expected_argmax[name], rtol, atol)    
@@ -205 +205 @@
     if not name in ['ch7', 'ch9']:
         # Minima of ch7 and ch9 are very flat and hard to pinpoint
-        i1 = argmin(model)
+        i1 = num.argmin(model)
         res = reference_time[i1]
         report_difference('Location of minimum', res, expected_argmin[name],

anuga_validation/automated_validation_tests/okushiri_tank_validation/create_okushiri.py

@@ -3 +3 @@
 
 
-from Numeric import array, zeros, Float, allclose
+import Numeric as num
 
 from anuga.pmesh.mesh import *
@@ -96 +96 @@
 
     from Scientific.IO.NetCDF import NetCDFFile
-    from Numeric import array
-
 
     print 'Creating', filename
@@ -113 +111 @@
 
     N = len(lines)
-    T = zeros(N, Float)  #Time
-    Q = zeros(N, Float)  #Values
+    T = num.zeros(N, num.Float)  #Time
+    Q = num.zeros(N, num.Float)  #Values
 
     for i, line in enumerate(lines):
@@ -130 +128 @@
     fid.starttime = 0.0
     fid.createDimension('number_of_timesteps', len(T))
-    fid.createVariable('time', Float, ('number_of_timesteps',))
+    fid.createVariable('time', num.Float, ('number_of_timesteps',))
     fid.variables['time'][:] = T
 
-    fid.createVariable('stage', Float, ('number_of_timesteps',))
+    fid.createVariable('stage', num.Float, ('number_of_timesteps',))
     fid.variables['stage'][:] = Q[:]
 
-    fid.createVariable('xmomentum', Float, ('number_of_timesteps',))
+    fid.createVariable('xmomentum', num.Float, ('number_of_timesteps',))
     fid.variables['xmomentum'][:] = 0.0
 
-    fid.createVariable('ymomentum', Float, ('number_of_timesteps',))
+    fid.createVariable('ymomentum', num.Float, ('number_of_timesteps',))
     fid.variables['ymomentum'][:] = 0.0