Changeset 2535


Ignore:
Timestamp:
Mar 13, 2006, 2:46:33 PM (18 years ago)
Author:
ole
Message:

Histogram of mesh areas

Files:
3 edited

Legend:

Unmodified
Added
Removed
  • development/okushiri_2005/lwru2.py

    r2426 r2535  
    115115
    116116
     117
    117118import sys, os
    118119base = os.path.basename(sys.argv[0])
     
    124125print 'Number of triangles = ', len(domain)
    125126print 'The extent is ', domain.get_extent()
     127print domain.statistics()
    126128
    127129
  • development/stochastic_study/run_model.py

    r2484 r2535  
    5050print 'Number of triangles = ', len(domain)
    5151print 'The extent is ', domain.get_extent()
     52print domain.statistics()
     53
     54import sys; sys.exit()
     55
    5256
    5357domain.set_datadir('.')
  • inundation/pyvolution/mesh.py

    r2533 r2535  
    631631        #Setup 10 bins for area histogram
    632632        m = max(areas)
    633         bins = arange(0, m/10, m)
     633        bins = arange(0., m, m/10)
    634634        hist = histogram(areas, bins)
    635635
     
    643643        str += '    A in [%f, %f]\n' %(min(areas), max(areas))
    644644        str += '    Histogram:\n'
     645
     646        hi = bins[0]
    645647        for i, count in enumerate(hist):
    646             str += '      >=%f: %d\n' %(bins[i], count)
     648            lo = hi
     649            if i+1 < len(bins):
     650                #Open upper interval               
     651                hi = bins[i+1]
     652                str += '      [%f, %f[: %d\n' %(lo, hi, count)               
     653            else:
     654                #Closed upper interval
     655                hi = m
     656                str += '      [%f, %f]: %d\n' %(lo, hi, count)               
     657               
    647658                     
    648659        str += 'Boundary:\n'
Note: See TracChangeset for help on using the changeset viewer.