\frac{\partial p}{\partial t} = -\nabla_x \cdot (f p) + \frac{1}{2} \nabla_x \nabla_x : (D p)