git.rot13.org Git - goodfet/blob - client/GoodFETGlitch.py

   1 #!/usr/bin/env python
   2 # GoodFET Client Library
   3 #
   4 # (C) 2009 Travis Goodspeed <travis at radiantmachines.com>
   5 #
   6 # This code is being rewritten and refactored.  You've been warned!
   7
   8 import sys, time, string, cStringIO, struct, glob, serial, os, random;
   9 import sqlite3;
  10
  11 from GoodFET import *;
  12
  13
  14 # After four million points, this kills 32-bit gnuplot.
  15 # Dumping to a bitmap might be preferable.
  16 script_timevcc="""
  17 plot "< sqlite3 glitch.db 'select time,vcc,glitchcount from glitches where count=0;'" \
  18 with dots \
  19 title "Scanned", \
  20 "< sqlite3 glitch.db 'select time,vcc,count from glitches where count>0;'" \
  21 with dots \
  22 title "Success", \
  23 "< sqlite3 glitch.db 'select time,vcc,count from glitches where count>0 and lock>0;'" \
  24 with dots \
  25 title "Exploited"
  26 """;
  27 script_timevccrange="""
  28 plot "< sqlite3 glitch.db 'select time,vcc,glitchcount from glitches where count=0;'" \
  29 with dots \
  30 title "Scanned", \
  31 "< sqlite3 glitch.db 'select time,vcc,count from glitches where count>0;'" \
  32 with dots \
  33 title "Success", \
  34 "< sqlite3 glitch.db 'select time,max(vcc),count from glitches where count=0 group by time ;'" with lines title "Max", \
  35 "< sqlite3 glitch.db 'select time,min(vcc),count from glitches where count>0 group by time ;'" with lines title "Min"
  36 """;
  37
  38 class GoodFETGlitch(GoodFET):
  39
  40     def __init__(self, *args, **kargs):
  41         print "Initializing GoodFET Glitcher."
  42         #Database connection w/ 30 second timeout.
  43         self.db=sqlite3.connect("glitch.db",30000);
  44
  45         #Training
  46         self.db.execute("create table if not exists glitches(time,vcc,gnd,trials,glitchcount,count,lock);");
  47         self.db.execute("create index if not exists glitchvcc on glitches(vcc);");
  48         self.db.execute("create index if not exists glitchtime on glitches(time);");
  49
  50         #Exploitation record, to be built from the training table.
  51         self.db.execute("create table if not exists exploits(time,vcc,gnd,trials,count);");
  52         self.db.execute("create index if not exists exploitvcc on exploits(vcc);");
  53         self.db.execute("create index if not exists exploittime on exploits(time);");
  54
  55         self.client=0;
  56     def setup(self,arch="avr"):
  57         self.client=getClient(arch);
  58         self.client.serInit();
  59
  60     def glitchvoltages(self,time):
  61         """Returns list of voltages to train at."""
  62         c=self.db.cursor();
  63         #c.execute("""select
  64         #             (select min(vcc) from glitches where time=? and count=1),
  65         #             (select max(vcc) from glitches where time=? and count=0);""",
  66         #          [time, time]);
  67         c.execute("select min,max from glitchrange where time=? and max-min>0;",[time]);
  68         rows=c.fetchall();
  69         for r in rows:
  70             min=r[0];
  71             max=r[1];
  72             if(min==None or max==None): return [];
  73
  74             spread=max-min;
  75             return range(min,max,1);
  76         #If we get here, there are no points.  Return empty set.
  77         return [];
  78     def buildglitchvoltages(self):
  79         """This builds tables for glitching voltage ranges from the training set."""
  80         print "Precomputing glitching ranges.  This might take a while.";
  81         sys.stdout.flush();
  82         self.db.execute("create temporary table glitchrange(time integer primary key asc,max,min);");
  83         self.db.execute("insert into glitchrange(time,max,min) select distinct time, 0, 0 from glitches;");
  84         self.db.execute("update glitchrange set max=(select max(vcc) from glitches where glitches.time=glitchrange.time and count=0);");
  85         self.db.execute("update glitchrange set min=(select min(vcc) from glitches where glitches.time=glitchrange.time and count>0);");
  86
  87     def graphx11(self):
  88         try:
  89             import Gnuplot, Gnuplot.PlotItems, Gnuplot.funcutils
  90         except ImportError:
  91             print "gnuplot-py is missing.  Can't graph."
  92             return;
  93         g = Gnuplot.Gnuplot(debug=1);
  94         g.clear();
  95
  96         g.title('Glitch Training Set');
  97         g.xlabel('Time (16MHz)');
  98         g.ylabel('VCC (DAC12)');
  99
 100         g('set datafile separator "|"');
 101
 102         g(script_timevcc);
 103         print "^C to exit.";
 104         while 1==1:
 105             time.sleep(30);
 106
 107
 108     def graph(self):
 109         import Gnuplot, Gnuplot.PlotItems, Gnuplot.funcutils
 110         g = Gnuplot.Gnuplot(debug=1);
 111
 112         g('\nset term png');
 113         g.title('Glitch Training Set');
 114         g.xlabel('Time (16MHz)');
 115         g.ylabel('VCC (DAC12)');
 116
 117         g('set datafile separator "|"');
 118         g('set term png');
 119         g('set output "timevcc.png"');
 120         g(script_timevcc);
 121     def explore(self,tstart=0,tstop=-1, trials=5):
 122         """Exploration phase.  Uses thresholds to find exploitable points."""
 123         gnd=0;
 124         self.scansetup(1); #Lock the chip, place key in eeprom.
 125         if tstop<0:
 126             tstop=self.client.glitchstarttime();
 127         times=range(tstart,tstop);
 128         random.shuffle(times);
 129         self.buildglitchvoltages();
 130         count=0.0;
 131         total=1.0*len(t);
 132         for t in times:
 133             voltages=self.glitchvoltages(t);
 134             count=count+1.0;
 135             print "%02.02f Exploring %04i points in t=%04i." % (count/total,len(voltages),t);
 136             sys.stdout.flush();
 137             for vcc in voltages:
 138                 self.scanat(1,trials,vcc,gnd,t);
 139     def learn(self):
 140         """Learning phase.  Finds thresholds at which the chip screws up."""
 141         trials=1;
 142         lock=0;  #1 locks, 0 unlocked
 143         vstart=0;
 144         vstop=1024;  #Could be as high as 0xFFF, but upper range is useless
 145         vstep=1;
 146         tstart=0;
 147         tstop=self.client.glitchstarttime();
 148         tstep=0x1; #Must be 1
 149         self.scan(lock,trials,range(vstart,vstop),range(tstart,tstop));
 150         print "Learning phase complete, beginning to expore.";
 151         self.explore();
 152
 153     def scansetup(self,lock):
 154         client=self.client;
 155         client.start();
 156         client.erase();
 157
 158         self.secret=0x69;
 159
 160         while(client.eeprompeek(0)!=self.secret):
 161             print "-- Setting secret";
 162             client.start();
 163
 164             #Flash the secret to the first two bytes of CODE memory.
 165             client.erase();
 166             client.eeprompoke(0,self.secret);
 167             client.eeprompoke(1,self.secret);
 168             sys.stdout.flush()
 169
 170         #Lock chip to unlock it later.
 171         if lock>0:
 172             client.lock();
 173
 174
 175     def scan(self,lock,trials,voltages,times):
 176         """Scan many voltages and times."""
 177         client=self.client;
 178         self.scansetup(lock);
 179         gnd=0;
 180         random.shuffle(voltages);
 181         #random.shuffle(times);
 182
 183         for vcc in voltages:
 184             if lock<0 and not self.vccexplored(vcc):
 185                 print "Exploring vcc=%i" % vcc;
 186                 sys.stdout.flush();
 187                 for time in times:
 188                     self.scanat(lock,trials,vcc,gnd,time)
 189                     sys.stdout.flush()
 190                 self.db.commit();
 191             else:
 192                 print "Voltage %i already explored." % vcc;
 193                 sys.stdout.flush();
 194
 195
 196     def vccexplored(self,vcc):
 197         c=self.db.cursor();
 198         c.execute("select vcc from glitches where vcc=? limit 1;",[vcc]);
 199         rows=c.fetchall();
 200         for a in rows:
 201             return True;
 202         return False;
 203     def scanat(self,lock,trials,vcc,gnd,time):
 204         client=self.client;
 205         db=self.db;
 206         client.glitchRate(time);
 207         client.glitchVoltages(gnd, vcc);  #drop voltage target
 208         gcount=0;
 209         scount=0;
 210         #print "-- (%5i,%5i)" % (time,vcc);
 211         #sys.stdout.flush();
 212         for i in range(0,trials):
 213             client.glitchstart();
 214
 215             #Try to read *0, which is secret if read works.
 216             a=client.eeprompeek(0x0);
 217             if lock>0: #locked
 218                 if(a!=0 and a!=0xFF and a!=self.secret):
 219                     gcount+=1;
 220                 if(a==self.secret):
 221                     print "-- %06i: %02x HELL YEAH! " % (time, a);
 222                     scount+=1;
 223             else: #unlocked
 224                 if(a!=self.secret):
 225                     gcount+=1;
 226                 if(a==self.secret):
 227                     scount+=1;
 228         #print "values (%i,%i,%i,%i,%i);" % (
 229         #    time,vcc,gnd,gcount,scount);
 230         if(lock>0):
 231             self.db.execute("insert into glitches(time,vcc,gnd,trials,glitchcount,count,lock)"
 232                    "values (%i,%i,%i,%i,%i,%i,%i);" % (
 233                 time,vcc,gnd,trials,gcount,scount,lock));
 234         else:
 235             self.db.execute("insert into exploits(time,vcc,gnd,trials,count)"
 236                    "values (%i,%i,%i,%i,%i);" % (
 237                 time,vcc,gnd,trials,scount));