X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FTools%2FVerima%2FStats%2Fmean.py;fp=src%2FTools%2FVerima%2FStats%2Fmean.py;h=0000000000000000000000000000000000000000;hp=cbf4f0a761ffe0587a48abd9db5fd18893bfde9b;hb=75d61804dc733ab5c810e416dec646406af11422;hpb=19a23c1d840a269deb7c4f6f549ac9bb942deb70 diff --git a/src/Tools/Verima/Stats/mean.py b/src/Tools/Verima/Stats/mean.py deleted file mode 100755 index cbf4f0a76..000000000 --- a/src/Tools/Verima/Stats/mean.py +++ /dev/null @@ -1,54 +0,0 @@ -#!/usr/bin/env python -# -*- coding: utf-8 -*- - -import os -import numpy, scipy -import scipy.stats - - -def getMean(fichierStatMaillage): - """ - """ - try : - monTableau=numpy.loadtxt(fichierStatMaillage) - except : - print "impossible de charger le fichier : ", fichierStatMaillage - - mesIntervalles=monTableau[ :, 0:2 ] - mesPoids=monTableau[ : ,2] - moyIntervalles=numpy.average(mesIntervalles, axis=1) - moyenne=numpy.average(moyIntervalles,weights=mesPoids) - freqCum=numpy.cumsum(mesPoids) - nbVal=freqCum[-1 ] - mesFrequences=mesPoids/nbVal - mesValeurs=scipy.stats.rv_discrete(values=(moyIntervalles,mesFrequences)) - Q1=mesValeurs.ppf(0.25) - M=mesValeurs.median() - Q3=mesValeurs.ppf(0.75) - - i=mesPoids.shape[0] -1 - while i > 0: - if mesPoids[i] > 0: - max=mesIntervalles[i][1] - break - i=i-1 - i=0 - while i < mesPoids.shape[0] -1: - if mesPoids[i] > 0: - min=mesIntervalles[i][0] - break - i=i+1 - i=0 - - return [max,min,Q1,M,Q3,moyenne] - - -if __name__ == "__main__": - from optparse import OptionParser - p=OptionParser() - p.add_option('-f',dest='fichier',default='tetra.taille',help='fichier a traiter') - #p.add_option('-f',dest='fichier',default='Mesh_1_aspect_ratio_3d.txt',help='fichier a traiter') - options, args = p.parse_args() - - getMean(options.fichier) -