Add function limitVelocityAndEnergy().
Make the Lax-Friedrichs flux a bit more efficient.
Some amendments to parameters in DoubleMachRefl.