 
 }
 
-
+//This is the predictor step of the finite volume scheme.
 __device__
 REALfour
 eulerStutterStep(REAL pfarLeft, REALfour stateLeft, REALfour stateCenter, REALfour stateRight, REAL pfarRight)

    REALthree fluxL, fluxR, pR;
    REALfour tempStateLeft, tempStateRight;

    pR = make_float3(pressureRatio(pfarLeft,stateLeft.w,stateCenter.w),
        pressureRatio(stateLeft.w,stateCenter.w,stateRight.w),
        pressureRatio(stateCenter.w,stateRight.w,pfarRight));

    tempStateLeft = limitor(make_float3(stateLeft), make_float3(stateCenter), pR.x);
    tempStateRight = limitor(make_float3(stateCenter), make_float3(stateLeft), 1.0/pR.y);
    tempStateLeft.w = pressure(tempStateLeft);
    tempStateRight.w = pressure(tempStateRight);
    fluxL = eulerFlux(tempStateLeft,tempStateRight);

    tempStateLeft = limitor(make_float3(stateCenter), make_float3(stateRight), pR.y);
    tempStateRight = limitor(make_float3(stateRight), make_float3(stateCenter), 1.0/pR.z);
    tempStateLeft.w = pressure(tempStateLeft);
    tempStateRight.w = pressure(tempStateRight);
    fluxR = eulerFlux(tempStateLeft,tempStateRight);

    stateCenter += make_float4(0.5 * dimens.x * (fluxL-fluxR));
    stateCenter.w = pressure(stateCenter);

    return stateCenter;
}

__device__
REALfour
eulerFinalStep(REAL pfarLeft, REALfour stateLeft, REALfour stateCenter, REALfour stateCenter_orig, REALfour stateRight, REAL pfarRight)

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
