Lanczos method for the numerical propagation of quantum densities in phase space with an application to the kicked harmonic oscillator