Simulation code for a 2D Ising ferromagnet with periodic boundary conditions using the Wolff algorithm