     cudaCheckErrors("cudaMalloc fail");
     
     // If we are going to time
-    bool timekernel=false;
-    cudaEvent_t start, stop;
-    float elapsedTime;
-    
-    
+
     int divx,divy,divz;
     
     // RB: Use the optimal (in their tests) block size from paper by Zinsser and Keck (16 in x and 32 in y).

     *
     *
     **/
    texture<float, cudaTextureType2DLayered , cudaReadModeElementType> tex;
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// RB, 10/31/2016: Add constant memory arrays to store parameters for all projections to be analyzed during a single kernel call
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

        // Copy the prepared parameter arrays to constant memory to make it available for the kernel
        cudaMemcpyToSymbol(projSinCosArrayDev, projSinCosArrayHost, sizeof(float)*5*PROJ_PER_KERNEL);
        cudaMemcpyToSymbol(projParamsArrayDev, projParamsArrayHost, sizeof(Point3D)*6*PROJ_PER_KERNEL);

        kernelPixelBackprojectionFDK<<<grid,block>>>(geo,dimage,i,nalpha);
        cudaCheckErrors("Kernel fail");

    }  // END for
    
    //////////////////////////////////////////////////////////////////////////////////////
