2 # -*- coding: utf-8 -*-
9 def getMean(fichierStatMaillage):
13 monTableau=numpy.loadtxt(fichierStatMaillage)
15 print "impossible de charger le fichier : ", fichierStatMaillage
17 mesIntervalles=monTableau[ :, 0:2 ]
18 mesPoids=monTableau[ : ,2]
19 moyIntervalles=numpy.average(mesIntervalles, axis=1)
20 moyenne=numpy.average(moyIntervalles,weights=mesPoids)
21 freqCum=numpy.cumsum(mesPoids)
23 mesFrequences=mesPoids/nbVal
24 mesValeurs=scipy.stats.rv_discrete(values=(moyIntervalles,mesFrequences))
25 Q1=mesValeurs.ppf(0.25)
27 Q3=mesValeurs.ppf(0.75)
29 i=mesPoids.shape[0] -1
32 max=mesIntervalles[i][1]
36 while i < mesPoids.shape[0] -1:
38 min=mesIntervalles[i][0]
43 return [max,min,Q1,M,Q3,moyenne]
46 if __name__ == "__main__":
47 from optparse import OptionParser
49 p.add_option('-f',dest='fichier',default='tetra.taille',help='fichier a traiter')
50 #p.add_option('-f',dest='fichier',default='Mesh_1_aspect_ratio_3d.txt',help='fichier a traiter')
51 options, args = p.parse_args()
53 getMean(options.fichier)