Ignore:
Timestamp:
Mar 23, 2006, 3:18:13 PM (18 years ago)
Author:
ole
Message:

Percentiles in mesh.py statistics

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/pyvolution/mesh.py

    r2535 r2582  
    642642        str += '  Areas:\n'
    643643        str += '    A in [%f, %f]\n' %(min(areas), max(areas))
     644        str += '    number of distinct areas: %d\n' %(len(areas))       
    644645        str += '    Histogram:\n'
    645646
     
    654655                #Closed upper interval
    655656                hi = m
    656                 str += '      [%f, %f]: %d\n' %(lo, hi, count)               
     657                str += '      [%f, %f]: %d\n' %(lo, hi, count)
     658
     659        N = len(areas)
     660        if N > 10:
     661            str += '    Percentiles (10%):\n'
     662            areas = areas.tolist()
     663            areas.sort()
     664
     665            k = 0
     666            lower = min(areas)
     667            for i, a in enumerate(areas):       
     668                if i % (N/10) == 0 and i != 0: #For every 10% of the sorted areas               
     669                    str += '      %d triangles in [%f, %f]\n' %(i-k, lower, a)
     670                    lower = a
     671                    k = i
     672                   
     673            str += '      %d triangles in [%f, %f]\n'\
     674                   %(N-k, lower, max(areas))                   
    657675               
    658676                     
Note: See TracChangeset for help on using the changeset viewer.