experiments/qian/analyseReport.py~

   1 lf = ['finalAnalysisReport_jf_oldci_stable_b_4bit.txt',  'finalAnalysisReport_jf_oldci_stable_b_6bit.txt','finalAnalysisReport_jf_oldci_stable_b_8bit.txt']
   2
   3 fd =[open(f,"r") for f in lf]
   4
   5 lignes=[fdi.readlines()[7:] for fdi in fd]
   6
   7 label = ["Frequency (Monobit)","Frequency  within a Block","Runs","Longest Run of Ones in a Block","Binary Matrix Rank","Discrete Fourier Transform (Spectral)","Non-overlapping Template Matching*","Overlapping Template Matching","Maurer's \"Universal Statistical\"","Linear Complexity","Serial* (m=10)","Approximate Entropy (m=10)","Cumulative Sums (Cusum) *","Random Excursions *","Random Excursions Variant *"]
   8
   9 endl="\\\\ \\hline \n"
  10
  11
  12 st=""
  13
  14 #freq
  15
  16 def traitenlignes(num1,num2):
  17     sti=""
  18     for l in lignes:
  19         #calcul des moyennes mv ms
  20         mv,ms=0,0
  21         for j in range(num1,num2+1):
  22
  23             stl=[]
  24             for x in l[j].split(" "):
  25                 if x !='':
  26                     stl += [x]
  27
  28             idxdiv = 0
  29             for i in range(len(stl)):
  30                 if stl[i].find('/') != -1:
  31                     idxdiv = i
  32             mv += float(stl[idxdiv-1])
  33             rate = stl[idxdiv].split("/")
  34             ms += float(rate[0])/float(rate[1])
  35
  36         #sti += "& " + str(mv/(num2-num1+1))[:5] + " (" + str(ms/(num2-num1+1))[0:5] +")"
  37         sti += "& " + str(mv/(num2-num1+1))[:5] + "/" + str(round(100*ms/(num2-num1+1)))[0:2]
  38
  39     return sti + endl
  40
  41
  42
  43 def traite1ligne(numl):
  44     sti=""
  45     for l in lignes:
  46         stl=[]
  47         splitee = l[numl].split(" ")
  48         for x in l[numl].split(" "):
  49             if x !='':
  50                 stl += [x]
  51
  52         idxdiv = 0
  53         for i in range(len(stl)):
  54             #print splitee[i],i
  55             if stl[i].find('/') != -1:
  56                 idxdiv = i
  57         #print idxdiv
  58
  59         r = str(round(float(stl[idxdiv].split("/")[0])))
  60         #ti += "& " + stl[idxdiv-1][:5] + " (" + r[0:5]+")"
  61         sti += "& " + stl[idxdiv-1][:5] + "/" + r[0:2]
  62     return sti + endl
  63
  64 st += label[0] + traite1ligne(0)
  65 st += label[1] + traite1ligne(1)
  66 st += label[2] + traite1ligne(4)
  67 st += label[3] + traite1ligne(5)
  68 st += label[4] + traite1ligne(6)
  69 st += label[5] + traite1ligne(7)
  70 st += label[6] + traitenlignes(8,155)
  71 st += label[7] + traite1ligne(156)
  72 st += label[8] + traite1ligne(157)
  73 st += label[9] + traite1ligne(185)
  74 st += label[10] + traitenlignes(185,186)
  75 st += label[11] + traite1ligne(187)
  76 st += label[12] + traitenlignes(2,3)
  77 st += label[13] + traitenlignes(159,166)
  78 st += label[14] + traitenlignes(167,184)
  79
  80
  81 print st
  82
  83
  84
  85
  86
  87