Standard 2D Ising model, multi-spin coding, Metropolis algorithm, written in CUDA