"""
from numpy import vstack, linspace, array
from pylab import plot, legend, show
from latgas import lattice
from latgas.lattice import Lattice, Potential, System
M_t = None
E_t = None
sizes = [1000, 500, 200, 100]
for sz in sizes:
inter = lambda r: -4
lattice = Lattice(sz, dim = 2)
potential = Potential(inter, 1.0, 5)
sys = System(lattice, potential)
sys.set_mu(8.0)
E = []
N = []
temp = linspace(2.5, 2.0, 26)
for i in temp:
print i
sys.set_T(i)
ee, nn = sys.run(1000)
ee, nn = sys.run(1000)
E.append(ee)
N.append(nn)
E = array(E)
N = array(N)
M = -(2*N - sz**2)/(sz**2)
if E_t == None and M_t == None:
measure to compare times.
"""
#!/usr/bin/env python
from latgas import lattice
from latgas.lattice import Lattice, Potential, System
from time import time
from pylab import plot, legend, xlabel, ylabel, show
from numpy import array
sz = 30
inter = lambda r: -4
lattice = Lattice(sz, dim = 2)
potential = Potential(inter, 1.0, 5)
sys = System(lattice, potential)
sys.set_mu(8.0)
sys.set_T(20.0)
steps = range(1, 201)
Trun = []
Tpy = []
Tc = []
for i in steps:
print i
t1 = time()
for j in xrange(i):
sys.whole_lattice()
t2 = time()
Tpy.append(1000*(t2-t1))
t1 = time()
