     vzr = mzr / dr;
     #ifdef DE //PRESSURE_DE
     E_kin = 0.5 * dr * ( vxr*vxr + vyr*vyr + vzr*vzr );
-    pr = Get_Pressure_From_DE( Er, Er - E_kin, dger, gamma );
+    pr = hydro_utilities::Get_Pressure_From_DE( Er, Er - E_kin, dger, gamma );
     #else
     pr  = (Er - 0.5*dr*(vxr*vxr + vyr*vyr + vzr*vzr)) * (gamma - 1.0);
     #endif //PRESSURE_DE

#include "../riemann_solvers/hllc_cuda.h"

#ifdef DE //PRESSURE_DE
#include "../hydro/hydro_cuda.h"
#endif



    vzl = mzl / dl;
    #ifdef DE //PRESSURE_DE
    E_kin = 0.5 * dl * ( vxl*vxl + vyl*vyl + vzl*vzl );
    pl = hydro_utilities::Get_Pressure_From_DE( El, El - E_kin, dgel, gamma );
    #else
    pl  = (El - 0.5*dl*(vxl*vxl + vyl*vyl + vzl*vzl)) * (gamma - 1.0);
    #endif //PRESSURE_DE
