         // 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);
-        if (timekernel){
-            cudaEventCreate(&start);
-            cudaEventRecord(start,0);
-        }
+
         kernelPixelBackprojectionFDK<<<grid,block>>>(geo,dimage,i,nalpha);
         cudaCheckErrors("Kernel fail");
-        if (timekernel)
-        {
-            cudaEventCreate(&stop);
-            cudaEventRecord(stop,0);
-            cudaEventSynchronize(stop);
-            cudaEventElapsedTime(&elapsedTime, start,stop);
-            mexPrintf("%f\n" ,elapsedTime);
-            cudaCheckErrors("cuda Timing fail");
-            
-        }
+
     }  // END for
     
     //////////////////////////////////////////////////////////////////////////////////////

     *
     *
     **/
    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
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

    cudaCheckErrors("cudaMalloc fail");
    
    // If we are going to time

    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).
