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: E_t = E M_t = M else: E_t = vstack([E_t, E]) M_t = vstack([M_t, M]) for i in M_t:
#!/usr/bin/env python """ test that everything works This should be superseeded by an eventual unittesting """ from latgas import lattice from latgas.lattice import Lattice, Potential, System sz = 30 def inter(r): if r <= 1.0: return -4 if r <= 3.0: return 1 lattice = Lattice(sz, dim = 2) potential = Potential(inter, 3.0, 20) sys = System(lattice, potential) sys.lattice.random(1) sys.set_mu(8.0) sys.set_T(20.0) sys.N = sys.tot_population() sys.E = sys.tot_energy() sys.run(3000) print "E = {0}, N = {1}".format(sys.E, sys.N)
""" #!/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() sys.run(i)
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: E_t = E M_t = M else: E_t = vstack([E_t, E]) M_t = vstack([M_t, M]) for i in M_t: