Mercurial > hg > medcouple

comparison pymedcouple @ 21:d63e763f8ac1

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
wmedian: don't infinite loop if wleft == wright
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Sat, 17 Jan 2015 19:54:53 -0500
parents 8e02346789b2
children f5b4a2ab6204
comparison
equal deleted inserted replaced
20:8e02346789b2 21:d63e763f8ac1
63 # weight. 63 # weight.
64 wright += w 64 wright += w
65 65
66 if 2*wleft > wtot: 66 if 2*wleft > wtot:
67 end = mid 67 end = mid
68 elif 2*wright <= wtot: 68 elif 2*wright < wtot:
69 beg = mid 69 beg = mid
70 else: 70 else:
71 return trial 71 return trial
72 72
73 def medcouple_1d(X, eps1 = 2**-52, eps2 = 2**-1022): 73 def medcouple_1d(X, eps1 = 2**-52, eps2 = 2**-1022):