  We present a multiscale thermalization algorithm for lattice gauge theory,
which enables efficient parallel generation of uncorrelated gauge field
configurations. The algorithm combines standard Monte Carlo techniques with
ideas drawn from real space renormalization group and multigrid methods. We
demonstrate the viability of the algorithm for pure Yang-Mills gauge theory for
both heat bath and hybrid Monte Carlo evolution, and show that it ameliorates
the problem of topological freezing up to controllable lattice spacing
artifacts.
