Sampling the thermal Wigner density via a generalized Langevin dynamics.