 const REAL m_gamma = 0.4;
 const REAL dx = 0.5;
 
+//dbd is the boundary condition
 __constant__ REALfour dbd[2]; //0 is left 1 is right.
+//dimens has three fields x is dt/dx, y is gamma, z is gamma-1
 __constant__ REALthree dimens;
 
+//Calculates the pressure at the current node with the rho, u, e state variables.
 __device__
 __forceinline__
 REAL


}


__device__
REALfour
eulerStutterStep(REAL pfarLeft, REALfour stateLeft, REALfour stateCenter, REALfour stateRight, REAL pfarRight)

}

//Left and Center then Left and right.
//This is the meat of the flux calculation.  Fields: x is rho, y is u, z is e, w is p.
__device__
REALthree
eulerFlux(REALfour cvLeft, REALfour cvRight)
{

    //For the first calculation rho and p remain the same.
    REALthree flux;
    REAL uLeft = cvLeft.y/cvLeft.x;
    REAL uRight = cvRight.y/cvRight.x;
