Changeset 5004 for anuga_core/source/anuga/geospatial_data

Timestamp:

Feb 7, 2008, 10:56:03 AM (17 years ago)

Author:

nick

Message:

add some more progress statements "if verbose: print...."

File:

• anuga_core/source/anuga/geospatial_data/geospatial_data.py (modified) (3 diffs)

anuga_core/source/anuga/geospatial_data/geospatial_data.py

@@ -1552 +1552 @@
         #returns the predicted elevation of the points that were "split" out 
         #of the original data set for one particular alpha
+        if verbose: print 'get predicted elevation for location to be compared'
         elevation_predicted=domain.quantities[attribute_smoothed].\
                             get_values(interpolation_points=points_geo)
@@ -1564 +1565 @@
         #print 'memory usage during compare',mem_usage()
         
-
-        if verbose: print'cov',normal_cov[i][0],'= ',normal_cov[i][1]
-
-
+        if verbose: print'covariance for alpha ',normal_cov[i][0],'= ',normal_cov[i][1]
 
 #    if verbose: print 'Determine difference between predicted results and actual data'
@@ -1602 +1600 @@
         remove(mesh_file)
     
+    if verbose: 
+        print 'Final results:'
+        for i in len(alphas):
+            print'covariance for alpha %s = %s ' %(normal_cov[i][0],normal_cov[i][1])
+        print '\n Optimal alpha is: %s ' % normal_cov_new[(argmin(normal_cov_new,axis=0))[1],0]
+
     return min(normal_cov_new[:,1]) , normal_cov_new[(argmin(normal_cov_new,axis=0))[1],0]