     E = dev_conserved[4*n_cells + id];
     GE = dev_conserved[(n_fields-1)*n_cells + id];
     E_kin = 0.5 * d * ( vx*vx + vy*vy + vz*vz );
-    P = Get_Pressure_From_DE( E, E - E_kin, GE, gamma );
+    P = hydro_utilities::Get_Pressure_From_DE( E, E - E_kin, GE, gamma );
     P  = fmax(P, (Real) TINY_NUMBER);
 
     imo = xid-1;

    E = dev_conserved[4*n_cells + id];
    GE = dev_conserved[(n_fields-1)*n_cells + id];
    E_kin = 0.5 * d * ( vx*vx + vy*vy + vz*vz );
    P = Get_Pressure_From_DE( E, E - E_kin, GE, gamma );
    P  = fmax(P, (Real) TINY_NUMBER);

    imo = xid-1 + yid*nx;

    E = dev_conserved[4*n_cells + id];
    GE = dev_conserved[(n_fields-1)*n_cells + id];
    E_kin = 0.5 * d * ( vx*vx + vy*vy + vz*vz );
    P = hydro_utilities::Get_Pressure_From_DE( E, E - E_kin, GE, gamma );
    P  = fmax(P, (Real) TINY_NUMBER);

    imo = xid-1 + yid*nx + zid*nx*ny;
