#test 2: check time average of L = 32 is 53.9 p1 = 0.5 p2 = 0.5 for l in [16, 32]: i = 0 tc = [] h1 = [] t = [] sandbox = OS.box(l, N, p1, p2) while i <= N: #drive sandbox.drive() sandbox.relax() if sandbox.steady() == True: h = 0 tc.append(i)
for jg in range(M): # the number of realisations #open a file named l fileout = open(f"{l}_{jg}", "wb") i = 0 h1 = [] t = [] tc = [] sandbox = OS.box(l, N, p1, p2) while i <= N: #drive sandbox.drive() sandbox.relax() tt = i if sandbox.steady() == True: tc.append(i)