 CUDA_reconstruct_conebeam_init(
   int img_dim[2],
   int vol_dim[3],
-  float *&dev_vol,             // Holds voxels on device
+  float *&dev_vol,         // Holds voxels on device
   float *&dev_img,         // Holds image pixels on device
-  float *&dev_matrix           // Holds matrix on device
+  float *&dev_matrix       // Holds matrix on device
   )
 {
   // Size of volume Malloc

#include <cuda.h>

// P R O T O T Y P E S ////////////////////////////////////////////////////
__global__ void kernel_fdk(float *dev_vol, int2 img_dim, int3 vol_dim, unsigned int Blocks_Y, float invBlocks_Y);

///////////////////////////////////////////////////////////////////////////


//_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_

__global__
void kernel_fdk(float *dev_vol, int3 vol_dim, unsigned int Blocks_Y, float invBlocks_Y)
{
  // CUDA 2.0 does not allow for a 3D grid, which severely
  // limits the manipulation of large 3D arrays of data.  The
