A general question to the Monte Carlo experts. When I use Wolff algorithm for global updates, say for Ising 2d, I always flip at least one spin (the initial random spin in the cluster). So, near the critical temperature when coming from below - where this should converge much faster than Metropolis' single flip - to a finite magnetization, we will always have strong osicillations because in each iteration we flip at least one spin with probability one. So we can't get rid of oscillations in this procedure, although they get weaker with higher temperature, where we can't flip more than one spin because of the low acceptance probability. Indeed in my matlab code I always see oscillations. Also near the critical temperature. What am I missing?