This is simple python code for generating samples from Ising model using Gibbs sampling.
Call get_sample(W, u, n) to get n samples from Ising:
For thining and applying burn-in you can change gibbs_sampling function in sample_ising.py.
For W as below, and
n=10,000, energy of samples plotted: