  A new Monte Carlo algorithm is introduced for the simulation of supercooled
liquids and glass formers, and tested in two model glasses. The algorithm is
shown to thermalize well below the Mode Coupling temperature and to outperform
other optimized Monte Carlo methods. Using the algorithm, we obtain finite size
effects in the specific heat. This effect points to the existence of a large
correlation length measurable in equal time correlation functions.
