gistfile1.cu

cuh
    void *dDensSum;
    void *dDensCount;
cu:init()
    gpuErrchk(cudaMalloc((void**)&dDensSum, sizeof(float)));
    gpuErrchk(cudaMalloc((void**)&dDensCount, sizeof(int)));
cu:~Points()
    gpuErrchk(cudaFree(dDensSum))
    gpuErrchk(cudafree(dDensCount))
cu: //Copy back to CPU
    cudaMemcpy(GPUInt, dDensSum, sizeof(float), cudaMemcpyDeviceToHost);
cu:runKernel
    cudaMemeset(dDensSum, 0, sizeof(float))
    cudaMemset(dDensCount, 0, sizeof(int))
kernel.cu: computeColorGradsGlobal
    //Add densSum & densCount as parameters
        if (norm(colorGrad) < 0.1)
        {
            atomicAdd(densSum, densities[idx]);
            atomicAdd(densCount, 1);
        }