Ising_square_Ice. samples the partition function of a square ice ising model H=J\sum_v Q_v, where Q_v=\sum_{i \in v} \sigma^{z}_i using single spin flips and short-loop updates