 float3 ip;
 float  voxel_data = 0;
 
-for (unsigned int proj = 0; proj<c_projSize[2]; 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]));

  float3 ip;
  float  voxel_data = 0;

  for (unsigned int proj = 0; proj<c_projSize[2]; proj++)
    {
    // matrix multiply
    ip = matrix_multiply(make_float3(i,j,k), &(c_matrices[12*proj]));

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

__constant__ float c_matrices[1024 * 12]; //Can process stacks of at most 1024 projections
__constant__ int3 c_projSize;

//_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
// K E R N E L S -_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_-_
