Exemplo n.º 1
0
import subprocess
import os.path

# SETUP ###################################################################################################
print '==Setup=='

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 100000

params.T_Start = 1
params.T_Stop = 3
params.T_StepSize = 0.1

Ising_L = [4,8,16,32,64]
Binning_K = [50,100,200,400,800]

useCache = True

p = Plotter(show = True, pdf = False, pgf = False, latex=False, name='3beta')

Tc = 2.27
v = -1

# FUNCTIONS ###############################################################################################
#===Error analysis===
def binning(allValues,k):
Exemplo n.º 2
0
# SETUP ###################################################################################################
print '==Setup=='

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 10000

#chose between cases 1, 2 and 3 and comment the others out

# standard parameters (used for tasks 2, 3)
params.T_Start = 1
params.T_Stop = 5
params.T_StepSize = 0.1

Ising_L = [4,16,64]
Binning_K = [50,200,800]

Ising_L = [4,16,64]
Binning_K = [50,200,800]

# Choose the right parameters according to task

# standard parameters
OverParam = 0

if OverParam != 0:
    # parameters to determine T_c and calculate the Binder parameter
Exemplo n.º 3
0
import os.path
from matplotlib.colors import LogNorm

# SETUP ###################################################################################################
print '==Setup=='

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 10000

params.T_Start = 1
params.T_Stop = 5
params.T_StepSize = 0.1

Ising_L = [4, 16, 64]
Binning_K = [50, 200, 800]

useCache = True

p = Plotter(show=True, pdf=False, pgf=False, latex=False, name='1simulation')

# FUNCTIONS ###############################################################################################
#==Exact===
itmp = np.arange(-1, 4 - 1)


def generateAllStates(n=4):
Exemplo n.º 4
0
import subprocess
import os.path

# SETUP ###################################################################################################
print '==Setup=='

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 10000

params.T_Start = 2.27 # << estimated critical Temperature
params.T_Stop = params.T_Start
params.T_StepSize = 0.01 # only 0.02 would be needed

Ising_L = [4,16,32,64,128]
Binning_K = [50,200,800,800,800]

useCache = True

p = Plotter(show = True, pdf = False, pgf = False, latex=False, name='3beta')

v=-1

# FUNCTIONS ###############################################################################################
#===Error analysis===
def binning(allValues,k):
    nBlocks=len(allValues)//k
Exemplo n.º 5
0
import subprocess
import os.path

# SETUP ###################################################################################################
print '==Setup=='

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 10000

params.T_Start = 2
params.T_Stop = 2.4
params.T_StepSize = 0.005 # only 0.02 would be needed

Ising_L = [4,16,64]
Binning_K = [50,200,800]

useCache = True

p = Plotter(show = True, pdf = False, pgf = False, latex=False, name='2binder')

# FUNCTIONS ###############################################################################################
#===Error analysis===
def binning(allValues,k):
    nBlocks=len(allValues)//k

    allBlocks = allValues[:nBlocks*k].reshape((-1,k))
Exemplo n.º 6
0
# SETUP ###################################################################################################
print "==Setup=="

params = Params()

params.Ising_J = 1
params.Ising_H = 0

params.MC_Seed = 42
params.MC_Sweeps = 10000

# chose between cases 1, 2 and 3 and comment the others out

# standard parameters (used for tasks 2, 3)
params.T_Start = 1
params.T_Stop = 5
params.T_StepSize = 0.1

Ising_L = [4, 16, 64]
Binning_K = [50, 200, 800]

Ising_L = [4, 16, 64]
Binning_K = [50, 200, 800]

# Choose the right parameters according to task

# standard parameters
OverParam = 0

if OverParam != 0:
    # parameters to determine T_c and calculate the Binder parameter