Multiple routines currently do their own job of calculating the equation of state P(e, rho); These need to all be cleaned up to use the macro definitions in cudaCommon.h, or to use the cuda kernel that returns an array of pressure.
Otherwise any attempt to alter the equation of state used in simulations is suicidal.