 #include "../reconstruction/plmp_cuda.h"
 
 #ifdef DE //PRESSURE_DE
-#include "../hydro/hydro_cuda.h"
+#include "../utils/hydro_utilities.h"
 #endif
 
 

    E = dev_conserved[4*n_cells + id];
    E_kin = 0.5 * d_ipo * ( vx_ipo*vx_ipo + vy_ipo*vy_ipo + vz_ipo*vz_ipo );
    dge = dev_conserved[(n_fields-1)*n_cells + id];
    p_ipo = Get_Pressure_From_DE( E, E - E_kin, dge, gamma );
    #else
    p_ipo  = (dev_conserved[4*n_cells + id] - 0.5*d_ipo*(vx_ipo*vx_ipo + vy_ipo*vy_ipo + vz_ipo*vz_ipo)) * (gamma - 1.0);
    #endif //PRESSURE_DE

    E = dev_conserved[4*n_cells + id];
    E_kin = 0.5 * d_i * ( vx_i*vx_i + vy_i*vy_i + vz_i*vz_i );
    dge = dev_conserved[(n_fields-1)*n_cells + id];
    p_i = hydro_utilities::Get_Pressure_From_DE( E, E - E_kin, dge, gamma );
    #else
    p_i  = (dev_conserved[4*n_cells + id] - 0.5*d_i*(vx_i*vx_i + vy_i*vy_i + vz_i*vz_i)) * (gamma - 1.0);
    #endif //PRESSURE_DE
