Changeset 2582

Timestamp:

Mar 23, 2006, 3:18:13 PM (17 years ago)

Author:

ole

Message:

Percentiles in mesh.py statistics

Location:

inundation

Files:

• examples/beach.py (modified) (1 diff)
• examples/bedslope.py (modified) (1 diff)
• pyvolution/mesh.py (modified) (2 diffs)

inundation/examples/beach.py

@@ -115 +115 @@
 print domain.get_extent()
 
+print domain.statistics()
 
 

inundation/examples/bedslope.py

@@ -21 +21 @@
 domain.set_quantities_to_be_stored('stage')  #See shallow_water.py
 
+print domain.statistics()
 
 #######################

inundation/pyvolution/mesh.py

@@ -642 +642 @@
         str += '  Areas:\n'
         str += '    A in [%f, %f]\n' %(min(areas), max(areas))
+        str += '    number of distinct areas: %d\n' %(len(areas))        
         str += '    Histogram:\n'
 
@@ -654 +655 @@
                 #Closed upper interval
                 hi = m
-                str += '      [%f, %f]: %d\n' %(lo, hi, count)                
+                str += '      [%f, %f]: %d\n' %(lo, hi, count)
+
+        N = len(areas)
+        if N > 10:
+            str += '    Percentiles (10%):\n'
+            areas = areas.tolist()
+            areas.sort()
+
+            k = 0
+            lower = min(areas)
+            for i, a in enumerate(areas):        
+                if i % (N/10) == 0 and i != 0: #For every 10% of the sorted areas                
+                    str += '      %d triangles in [%f, %f]\n' %(i-k, lower, a)
+                    lower = a
+                    k = i
+                    
+            str += '      %d triangles in [%f, %f]\n'\
+                   %(N-k, lower, max(areas))