In [1]:

    

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
from IsingModel.Ising1 import IsingSimple


    

In [15]:

    

model = IsingSimple(n=256,a=0,b=10)
model.randomflips(0.3)
print(model.energy())
model.savegrid("hello")


    

    




0






    

In [16]:

    

for a in [.1,.3,.5,.7,.9,2]:
    model = IsingSimple(n=100,a=a,b=0)
    for percent in [0.1,0.3,.5,.7,.9]


    

    
Out[16]:





array([[ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546],
       [ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546],
       [ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546],
       ..., 
       [ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546],
       [ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546],
       [ 0.9999546,  0.9999546,  0.9999546, ...,  0.9999546,  0.9999546,
         0.9999546]])

In [ ]:

    




    

In [ ]: