   float3 ip;
   float  voxel_data = 0;
 
-  for (unsigned int proj = 0; proj<c_proj_size.z; proj++)
+  for (unsigned int proj = 0; proj<c_projSize.z; proj++)
     {
     // matrix multiply
     ip = matrix_multiply(make_float3(i,j,k), &(c_matrices[12*proj]));


// Constant memory
__constant__ float c_matrices[1024 * 12]; //Can process stacks of at most 1024 projections
__constant__ int3 c_proj_size;
__constant__ int3 c_vol_size;

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


  // Place it into the volume
  dev_vol_out[vol_idx] = dev_vol_in[vol_idx] + voxel_data;
}

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