  We give a survey of the known results on mixing time of Glauber dynamics for
the Ising model on the square lattice and present a technique that makes exact
sampling of the Ising model at all temperatures possible in polynomial time. At
high temperatures this is well-known and although this seems to be known also
in the low temperature case since Kramer and Waniers paper from the 1950s, we
did not found any reference that describes exact sampling for the Ising model
at low temperatures.
