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

  float3 ip;
  float  voxel_data = 0;

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


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

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