In [1]:

    

include("../src/Wang_Landau.jl")


    

    
Out[1]:





explore_phase_space (generic function with 4 methods)

In [2]:

    

using PyPlot


    

    




INFO: Loading help data...

In [2]:

    

@time xs, ys, thetas, free_paths, free_path_distribution = explore_phase_space(0.3, 0.1, 1000000);


    

    




elapsed time: 5.708427309 seconds (3763151600 bytes allocated, 33.75% gc time)

In [38]:

    

plot(xs, ys, alpha=0.5)
axis("image")


    

    













    
Out[38]:





(-0.4999977560526506,0.4999974815262893,-0.49999777387548505,0.49999276906109147)

In [48]:

    




    

    
Out[48]:





(0.0:1.0:100.0,[639757,218668,60575,24485,13138,7843,4810,3351,2593,2118  …  15,26,27,28,28,25,20,30,14,14])

In [11]:

    

h = hist(free_paths[free_paths.<1000], 100)
h1, h2 = h
loglog(h1[1:end-1], h2)


    

    













    
Out[11]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x1149002d0>

In [8]:

    

plot(free_path_distribution, "o-")
yscale("log")
xscale("log")


    

In [41]:

    

thetas


    

    
Out[41]:





99819-element Array{Float64,1}:
 -1.8368 
 -1.83799
 -1.83409
 -1.83057
 -1.83382
 -1.83162
 -1.829  
 -1.832  
 -1.83142
 -1.83076
 -1.82999
 -1.83209
 -1.83348
  ⋮      
 -1.85933
 -1.85702
 -1.85762
 -1.85277
 -1.85107
 -1.85129
 -1.85342
 -1.85339
 -1.85762
 -1.85557
 -1.86031
 -1.86319

In [42]:

    

plot(thetas)


    

    













    
Out[42]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x11bd57110>

In [ ]: