astronomix._fluid_equations._enforce_positivity

Here we protect the density and pressure from going negative.

In my view this is a bit of a shady practice, hiding unphysical updates under the rug. However, it is common practice.