/usr/share/ufo/backproject.cl is in ufo-filters-data 0.14.1+dfsg1-1.1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 | /*
* Copyright (C) 2011-2013 Karlsruhe Institute of Technology
*
* This file is part of Ufo.
*
* This library is free software: you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation, either
* version 3 of the License, or (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with this library. If not, see <http://www.gnu.org/licenses/>.
*/
constant sampler_t volumeSampler = CLK_NORMALIZED_COORDS_FALSE |
CLK_ADDRESS_CLAMP |
CLK_FILTER_LINEAR;
kernel void
backproject_nearest (global float *sinogram,
global float *slice,
constant float *sin_lut,
constant float *cos_lut,
const unsigned int x_offset,
const unsigned int y_offset,
const unsigned int angle_offset,
const unsigned n_projections,
const float axis_pos)
{
const int idx = get_global_id(0);
const int idy = get_global_id(1);
const int width = get_global_size(0);
const float bx = idx - axis_pos + x_offset + 0.5f;
const float by = idy - axis_pos + y_offset + 0.5f;
float sum = 0.0f;
for(int proj = 0; proj < n_projections; proj++) {
float h = axis_pos + bx * cos_lut[angle_offset + proj] + by * sin_lut[angle_offset + proj];
sum += sinogram[(int)(proj * width + h)];
}
slice[idy * width + idx] = sum * M_PI_F / n_projections;
}
kernel void
backproject_tex (read_only image2d_t sinogram,
global float *slice,
constant float *sin_lut,
constant float *cos_lut,
const unsigned int x_offset,
const unsigned int y_offset,
const unsigned int angle_offset,
const unsigned int n_projections,
const float axis_pos)
{
const int idx = get_global_id(0);
const int idy = get_global_id(1);
const float bx = idx - axis_pos + x_offset + 0.5f;
const float by = idy - axis_pos + y_offset + 0.5f;
float sum = 0.0f;
#ifdef DEVICE_GEFORCE_GTX_TITAN_BLACK
#pragma unroll 8
#endif
#ifdef DEVICE_GEFORCE_GTX_TITAN
#pragma unroll 14
#endif
for(int proj = 0; proj < n_projections; proj++) {
float h = by * sin_lut[angle_offset + proj] + bx * cos_lut[angle_offset + proj] + axis_pos;
sum += read_imagef (sinogram, volumeSampler, (float2)(h, proj + 0.5f)).x;
}
slice[idy * get_global_size(0) + idx] = sum * M_PI_F / n_projections;
}
|