These functions works if rho (density of particles) and u (specific energy of particles) are correctely defined.