Mercurial > hg > medcouple

diff talk/code/show_mc.py @ 67:1aa38f4846e0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
showalgo: show the rank
author Jordi Gutiérrez Hermoso <jordigh@octave.org>
date Wed, 18 May 2016 17:43:45 -0400
parents 916349600c4b
children 5510ac95a3d3
line wrap: on
line diff
--- a/talk/code/show_mc.py
+++ b/talk/code/show_mc.py
@@ -95,6 +95,7 @@
     RemFile = open("remaining", "w")
     Ptotfile = open("Ptotal", "w")
     Qtotfile = open("Qtotal", "w")
+    medc_idxfile = open("medc_idx", "w")
     jackpot = open("jackpot", "w")
     
     # FIXME: Figure out what to do about NaNs.
@@ -166,6 +167,7 @@
     Ltot = 0
     Rtot = n_minus*n_plus
     medc_idx = Rtot//2
+    medc_idxfile.write("%s\n" % medc_idx)
 
     Lfile.write("%s\n" % L)
     Rfile.write("%s\n" % R)
@@ -259,7 +261,7 @@
 
 
 def main():
-    data = np.random.uniform(size=18)
+    data = np.random.uniform(size=26)
 
     mc = medcouple_1d(data)
     with open("mc", "w") as f: