    //   logger.flush();
    //}
    // END DEBUG
-}
-
+
+   first = false;
+}
+

#include "particle_allocator.h"
#include "pic_utils.h"

//texture<float, 1, cudaReadModeElementType > exTex;
//texture<float, 1, cudaReadModeElementType > eyTex;

__global__
void checkBoundaryConditions(float2 location[], int numParticles, int height, int width, 

   checkForCudaError("Before moveParticles");
   moveParticles<<<numBlocks, blockSize, 0, stream>>>(
      partLoc.getPtr(), partVel.getPtr(),
      numParticles, mass, NX1, NY1);
   checkForCudaError("moveParticles");

