 
 
 // ----------------------------------------------------------------------------
+/* Compute the norm of a table with 3 double */
 __device__ double norm(double vec[3])
 {
   return sqrt(vec[0] * vec[0] + vec[1] * vec[1] + vec[2] * vec[2]);
 }
 
 // ----------------------------------------------------------------------------
+/* Ray potential function which computes the increment to the current voxel */
 __device__ double cumulFunction(double diff, double currentVal)
 {
   double shift = 10 - 0.5 * std::abs(diff);



// ----------------------------------------------------------------------------
__host__ void doubleTableToVtkDoubleArray(double* table, vtkDoubleArray* output)
{
  int nbVoxels = output->GetNumberOfTuples();

                   int h_gridDims[3], // Dimensions of the output volume
                   double h_gridOrig[3], // Origin of the output volume
                   double h_gridSpacing[3], // Spacing of the output volume
                   vtkDoubleArray* io_outScalar) // It will be filled inside function
{
  if (h_dataList.size() == 0)
    return -1;
