from math import sqrt from time import clock from sys import exit import explore from explore import * print trials,size,perturb,eps,fill_level,bound print bigdata = eval(open("sle.dat").read()) big = int ( 0.5 + 13.5 * (size*bound)**2 ) rough = int ( 0.5 + 13.5 * (size*eps)**2 ) def work(): explore.clr = {} boundary() if fill_level > 0: fill_all() explor1 = exploration() forget_middle() if fill_level > 0: fill_middle() explor2 = exploration() explor1 = shorten(explor1,rough) explor2 = shorten(explor2,rough) l1,l2 = len(explor1),len(explor2) explor1 = dict(zip(range(l1),explor1)) explor2 = dict(zip(range(l2),explor2)) d = distance(explor1,explor2,big) return d for trial in range(trials[0],trials[1]): sd = (size_code*16+perturb_code+1, (1-fill_level)*16+eps_code+1, trial+1) whrandom.seed(sd[0],sd[1],sd[2]) d = work() print trial,d,clock() bigdata[(size_code,perturb_code,eps_code,trial)] = d print >>open("sle.dat","w"), bigdata print "clock: ", clock()