scan/parseDat.py

   1 #!/usr/bin/env python2
   2
   3 import os, subprocess
   4
   5 def listtofile(listing, file):
   6         fd = open(file, 'r')
   7         for i in listing:
   8                 fd.write(i+"\n")
   9         fd.close()
  10
  11 def makeGraphs(graphs, graphtitle = "Title", xlabel = "Numbers calculated", ylabel = "Time needed for calculation"):
  12         gendir = "stats/generated/"
  13         subprocess.call(["mkdir", "-p", gendir])
  14         p = subprocess.Popen(['gnuplot'],stdout=subprocess.PIPE,stdin=subprocess.PIPE)
  15         # now do the plots
  16         for graphfile in graphs:
  17                 # set output 'filename'
  18                 p.stdin.write("set terminal png font arial 8 #size 600,300\n")
  19                 p.stdin.write("set grid\n")
  20                 p.stdin.write("set datafile separator ';'\n")
  21                 p.stdin.write("set title 'Timing of each "+graphtitle+"'\n")
  22                 p.stdin.write("set xlabel '"+xlabel+"'\n")
  23                 p.stdin.write("set ylabel '"+ylabel+"'\n")
  24                 p.stdin.write("set output '"+gendir+graphfile+"'\n")
  25                 # plot 'filename' using 1:2 title 'algorithm' with lines, 'filename2' using 1:2 title 'algo2' with lines
  26                 plots = []
  27                 for title, file, cols in graphs[graphfile]:
  28                         plots.append("'"+gendir+file+"' using "+str(cols)+" title '"+str(title)+"' with lines")
  29                 plotcommand = ", ".join(plots)
  30                 #print gendir+graphfile
  31                 p.stdin.write("plot "+plotcommand+"\n")
  32         p.communicate()[0]
  33         p.stdin.close()
  34
  35 statsdir = 'stats/raw/'
  36 runs = os.listdir(statsdir)
  37 #print runs
  38
  39 bigrunlist = []
  40 for i in runs:
  41         runfiles = os.listdir(statsdir+i)
  42 #       print runfiles
  43         for file in runfiles:
  44 #               print "File: %s%s/%s" % (statsdir, i, file)
  45                 splitfilename = file.split('_')
  46                 dict = {}
  47                 dict['no'] = int(splitfilename[0].lstrip('no'))
  48                 dict['nnp'] = int(splitfilename[1].lstrip('nnp'))
  49                 dict['n'] = int(splitfilename[2].lstrip('n'))
  50                 fd = open(statsdir+i+'/'+file, "r")
  51                 fd_text = fd.read()
  52                 fd.close()
  53                 fd_lines = fd_text.split("\n")
  54                 for line in fd_lines:
  55                         if line == "":
  56                                 continue
  57                         elements = line.split(" ")
  58                         if len(elements) < 2:
  59                                 continue
  60                         if elements[1] != "timings:":
  61                                 continue
  62                         for elem in elements:
  63                                 (key, value) = elem.split(":", 1)
  64                                 if key == 'prep':
  65                                         dict['prep'] = float(value)
  66                                 elif key == 'algo':
  67                                         dict['algo'] = float(value)
  68                                 elif key == 'postp':
  69                                         dict['postp'] = float(value)
  70                                 elif key == 'end':
  71                                         dict['end'] = float(value)
  72                         bigrunlist.append(dict)
  73 nos = {}
  74 nnps = {}
  75 ns = {}
  76 foolist = {}
  77 for i in bigrunlist:
  78         if not (i['no'], i['nnp'], i['n']) in foolist:
  79                 foolist[(i['no'], i['nnp'], i['n'])] = [i['end']]
  80         else:
  81                 foolist[(i['no'], i['nnp'], i['n'])].append(i['end'])
  82         if not i['no'] in nos:
  83                 nos[i['no']] = {}
  84         if not i['nnp'] in nnps:
  85                 nnps[i['nnp']] = {}
  86         if not i['n'] in ns:
  87                 ns[i['n']] = {}
  88 #print "NOs: %s" % nos
  89 #print "NNPs: %s" % nnps
  90 #print "Ns: %s" % ns
  91 uniquelist = []
  92 for key in foolist:
  93         minval = min(foolist[key])
  94         maxval = max(foolist[key])
  95         avgval = sum(foolist[key])/len(foolist[key])
  96         uniquelist.append((key[0], key[1], key[2], minval))
  97 uniquelist.sort()
  98 alg_dict = {}
  99 sched_dict = {}
 100 for i in uniquelist:
 101         if i[0] not in alg_dict:
 102                 alg_dict[i[0]] = {}
 103         if i[1] not in alg_dict[i[0]]:
 104                 alg_dict[i[0]][i[1]] = []
 105         alg_dict[i[0]][i[1]].append(str(i[2])+";"+str(i[3]))
 106         if i[1] not in sched_dict:
 107                 sched_dict[i[1]] = {}
 108         if i[0] not in sched_dict[i[1]]:
 109                 sched_dict[i[1]][i[0]] = []
 110         sched_dict[i[1]][i[0]].append(str(i[2])+";"+str(i[3]))
 111         #print "od["+i[0]+"]["+i[1]+"] = "+str(i[2])+":"+str(i[3])
 112
 113 graphs_alg = {}
 114 graphs_sched = {}
 115 gendir = "stats/generated/"
 116 subprocess.call(["mkdir", "-p", gendir])
 117 for alg in alg_dict:
 118         pngfile = 'no.'+str(alg)+'.png'
 119         graphs_alg[pngfile] = []
 120         for sched in alg_dict[alg]:
 121                 datfile = 'data.no'+str(alg)+'.nnp'+str(sched)+'.dat'
 122                 subprocess.call(["rm", "-f", gendir+datfile])
 123                 fh = open(gendir+datfile, 'a+')
 124                 graphs_alg[pngfile].append((sched, datfile, '1:2'))
 125                 for dot in alg_dict[alg][sched]:
 126                         fh.write(dot+"\n")
 127                 fh.close()
 128 for sched in sched_dict:
 129         pngfile = 'nnp.'+str(sched)+'.png'
 130         graphs_sched[pngfile] = []
 131         for alg in sched_dict[sched]:
 132                 datfile = 'data.no'+str(alg)+'.nnp'+str(sched)+'.dat'
 133                 graphs_sched[pngfile].append((alg, datfile, '1:2'))
 134
 135 #graphs['test1.png'] = [('linetitle1', 'prefix.hillis.dat', '1:2')]
 136
 137 makeGraphs(graphs_alg, "no")
 138 makeGraphs(graphs_sched, "nnp")