GPU based Monte-Carlo simulation of the Ising model on the square 512×512 grid. it is the number of MC iterations per frame. The thermostat algorithm and linear congruential random number generator are used.
For low temperature you can watch formation of clusters of ordered phase from random state. For T > TC = 2.269 you will get desordered phase, for T < 0 the antiferromagnetic model. See Phase transitions on lattices.