Statistically independent equilibrium configurations of interacting systems can easily be generated with little computing effort by multigrid coarse-to-fine transforms. The very slow long-wavelength motions close to critical points are thereby accelerated by many orders of magnitude. This is demonstrated with a novel updating algorithm from small to large two-dimensional Ising systems specifying effective interactions according to renormalisation group considerations, corrected for non-linearities and saturation effects.