Data Analysis for Ising model, 2-dim.

In [1]:

    

import numpy
import numpy as np


    

In [2]:

    

avgdatperT = np.fromfile("./data/IsingMetroCPU.bin",dtype=np.float64) # averages data for each temperature T


    

In [3]:

    

print(type(avgdatperT))  
# from output.h, expect to output each temperature T, 
# temperature T, Eavg, Evariance, Mavg, Mvariance, |M|_avg, so 6 parameters  
print(avgdatperT.size)


    

    




<type 'numpy.ndarray'>
246

In [7]:

    

avgdatperT1 = avgdatperT.reshape( avgdatperT.size/6,6)
#print(avgdatperT1);


    

In [8]:

    

import matplotlib.pyplot as plt


    

In [9]:

    

fig = plt.figure()
ax = fig.add_subplot(1,1,1)
T = avgdatperT1[:,0]
E_avg = avgdatperT1[:,1]
ax.scatter( T, E_avg)
plt.show()


    

In [10]:

    

Evar_avg = avgdatperT1[:,2]
plt.scatter( T, Evar_avg)
plt.show()


    

In [11]:

    

M_avg = avgdatperT1[:,3]
Mvar_avg = avgdatperT1[:,4]
absM_avg = avgdatperT1[:,5]
plt.scatter( T, M_avg)
plt.show()


    

In [ ]: