The variables for the equations in conservation form are not yet optimised. In fact we could define: y = ( ρ ρ j 0 ) ; F = ( j 1 ρ ρ j ⊗ ⊗ j + p I j ρ ρ ) . where j = ρ ρ u is the momentum density, a conservation variable.