 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

    // Note: cbi->img AND cbi->matrix are passed via texture memory
    //-------------------------------------
    kernel_fdk_optim <<< dimGrid, dimBlock >>> ( dev_vol,
                                                 make_int2(img_dim[0], img_dim[1]),
                                                 make_int3(vol_dim[0], vol_dim[1], vol_dim[2]) );
    }
  else

    // Note: cbi->img AND cbi->matrix are passed via texture memory
    //-------------------------------------
    kernel_fdk <<< dimGrid, dimBlock >>> ( dev_vol,
                                           make_int3(vol_dim[0], vol_dim[1], vol_dim[2]),
                                           blocksInY, 1.0f/(float)blocksInY );
    }
