The above smoothenings are directly carried over to the multi-dimensional densities obtained as direct products. Thus the multi-dimensional densities still do not need to be constructed, but instead computed when needed. Altogether, all required densities can in this manner be very efficiently evaluated from the one-mode density matrix which is in almost all cases anyway needed during the time-propagation of the wave function.