Changeset 4667

Timestamp:

Aug 9, 2007, 9:50:57 AM (18 years ago)

Author:

nick

Message:

removed old data file

Location:

anuga_work/development/alpha_validation

Files:

• find_alpha.py (modified) (5 diffs)
• pt_hedland_export.txt (deleted)

anuga_work/development/alpha_validation/find_alpha.py

@@ -7 +7 @@
 from anuga.shallow_water import Reflective_boundary
 from anuga.utilities.numerical_tools import cov
-from Numeric import array, resize,shape,Float,zeros,take,argsort
-from pylab import plot, ion, hold,savefig,semilogx,plotting
+from Numeric import array, resize,shape,Float,zeros,take,argsort,argmin
+from pylab import plot, ion, hold,savefig,semilogx,plotting,loglog
 
 
@@ -18 +18 @@
                          boundary_tags={'back': [2], 'side': [1,3],
                                         'ocean': [0]},
-                         maximum_triangle_area=100,
+                         maximum_triangle_area=1000,
                          filename=mesh_dir_name,
                          use_cache=True,
@@ -24 +24 @@
 
 
-topo_dir_name = 'pt_hedland_export.txt'
-sample = 'sample.txt'
+#topo_dir_name = 'pt_hedland_export.txt'
+topo_dir_name = 'pt_hedland_small.txt'
+#sample = 'pt_hedland_small.txt'
 remainder ='remainder.txt'
 #split topo data
@@ -31 +32 @@
 print 'start split'
 
-G_sample,G_loss= G.split(0.1, True)
+
+
+
+
+G_sample,G_loss= G.split(0.01, True)
+
+#G_sample.export_points_file(sample)
 
 G_small, G_other = G_sample.split(0.1,True)
 
-
+#import sys
+#sys.exit()
 
 alphas = [0.000001, 0.00001, 0.0001, 0.001, 0.01, 0.1, 1.0, 10.0, 100.0,1000.0,10000.0]
@@ -93 +101 @@
 print'data',data[100:200]
 print'normal cov',normal_cov
+#to sort array by column
 normal_cov0=normal_cov[:,0]
 normal_cov_new=take(normal_cov,argsort(normal_cov0))
 semilogx(normal_cov_new[:,0],normal_cov_new[:,1])
+#loglog(normal_cov_new[:,0],normal_cov_new[:,1])
+print 'normal_cov_new',normal_cov_new
 savefig("alphas",dpi=300)
+
+print argmin(normal_cov_new,axis=1)
+print min(normal_cov_new[:,1]) , normal_cov_new[(argmin(normal_cov_new,axis=0))[1],0]