The simulation step is based on the cluster MC method, described in detail in [22]. Figure 2 schematically shows the single iteration step including the cluster building procedure as the modification of the Metropolis algorithm.