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):
# 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
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):
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
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))
# 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