After calculating the update to the pressure (conservation of mass) and temperature (conservation of energy) the density is updated:
is obtained from and :
whereas satisfies