+ self.client.serInit();
+
+ def glitchvoltages(self,time):
+ """Returns list of voltages to train at."""
+ c=self.db.cursor();
+ #c.execute("""select
+ # (select min(vcc) from glitches where time=? and count=1),
+ # (select max(vcc) from glitches where time=? and count=0);""",
+ # [time, time]);
+ c.execute("select min,max from glitchrange where time=? and max-min>0;",[time]);
+ rows=c.fetchall();
+ for r in rows:
+ min=r[0];
+ max=r[1];
+ if(min==None or max==None): return [];
+
+ spread=max-min;
+ return range(min,max,1);
+ #If we get here, there are no points. Return empty set.
+ return [];
+ def crunch(self):
+ """This builds tables for glitching voltage ranges from the training set."""
+ print "Precomputing glitching ranges. This might take a long while.";
+ print "Times...";
+ sys.stdout.flush();
+ self.db.execute("drop table if exists glitchrange;");
+ self.db.execute("create table glitchrange(time integer primary key asc,max,min);");
+ self.db.execute("insert into glitchrange(time,max,min) select distinct time, 0, 0 from glitches;");
+ print "Maximums...";
+ sys.stdout.flush();
+ self.db.execute("update glitchrange set max=(select max(vcc) from glitches where glitches.time=glitchrange.time and count=0);");
+ print "Minimums...";
+ sys.stdout.flush();
+ self.db.execute("update glitchrange set min=(select min(vcc) from glitches where glitches.time=glitchrange.time and count>0);");
+ print "Ranges calculated.";
+ def graphx11(self):
+ try:
+ import Gnuplot, Gnuplot.PlotItems, Gnuplot.funcutils
+ except ImportError:
+ print "gnuplot-py is missing. Can't graph."
+ return;
+ g = Gnuplot.Gnuplot(debug=1);
+ g.clear();
+
+ g.title('Glitch Training Set');
+ g.xlabel('Time (16MHz)');
+ g.ylabel('VCC (DAC12)');
+
+ g('set datafile separator "|"');
+
+ g(script_timevcc);
+ print "^C to exit.";
+ while 1==1:
+ time.sleep(30);
+
+
+ def graph(self):
+ import Gnuplot, Gnuplot.PlotItems, Gnuplot.funcutils
+ g = Gnuplot.Gnuplot(debug=1);
+
+ g('\nset term png');
+ g.title('Glitch Training Set');
+ g.xlabel('Time (16MHz)');
+ g.ylabel('VCC (DAC12)');
+
+ g('set datafile separator "|"');
+ g('set term png');
+ g('set output "timevcc.png"');
+ g(script_timevcc);
+ def explore(self,tstart=0,tstop=-1, trials=5):
+ """Exploration phase. Uses thresholds to find exploitable points."""
+ gnd=0;
+ self.scansetup(1); #Lock the chip, place key in eeprom.
+ if tstop<0:
+ tstop=self.client.glitchstarttime();
+ times=range(tstart,tstop);
+ random.shuffle(times);
+ #self.crunch();
+ count=0.0;
+ total=1.0*len(times);
+ for t in times:
+ voltages=self.glitchvoltages(t);
+ count=count+1.0;
+ print "%02.02f Exploring %04i points in t=%04i." % (count/total,len(voltages),t);
+ sys.stdout.flush();
+ for vcc in voltages:
+ self.scanat(1,trials,vcc,gnd,t);