do

density = initial(density)