diff options
Diffstat (limited to 'silx/resources/opencl/medfilt.cl')
-rw-r--r-- | silx/resources/opencl/medfilt.cl | 141 |
1 files changed, 0 insertions, 141 deletions
diff --git a/silx/resources/opencl/medfilt.cl b/silx/resources/opencl/medfilt.cl deleted file mode 100644 index 0680230..0000000 --- a/silx/resources/opencl/medfilt.cl +++ /dev/null @@ -1,141 +0,0 @@ -/* - * Project: Azimuthal regroupping OpenCL kernel for PyFAI. - * Median filter for 1D, 2D and 3D datasets, only 2D for now - * - * - * Copyright (C) 2017-2017 European Synchrotron Radiation Facility - * Grenoble, France - * - * Principal authors: J. Kieffer (kieffer@esrf.fr) - * Last revision: 07/02/2017 - * - * Permission is hereby granted, free of charge, to any person obtaining a copy - * of this software and associated documentation files (the "Software"), to deal - * in the Software without restriction, including without limitation the rights - * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell - * copies of the Software, and to permit persons to whom the Software is - * furnished to do so, subject to the following conditions: - * - * The above copyright notice and this permission notice shall be included in - * all copies or substantial portions of the Software. - * - * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR - * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, - * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE - * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER - * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, - * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN - * THE SOFTWARE. - * - */ - -/* - * Needs to be concatenated with bitonic.cl prior to compilation -*/ - -/* - * Perform the 2D median filtering of an image 2D image - * - * dim0 => wg=number_of_element in the tile /8 - * dim1 = x: wg=1 - - * - * Actually the workgoup size is a bit more complicated: - * if window size = 1,3,5,7: WG1=8 - * if window size = 9,11,13,15: WG1=32 - * if window size = 17, ...,21: WG1=64 - * - * More Generally the workgroup size must be: at least: kfs2*(kfs1+7)/8 - * Each thread treats 8 values aligned vertically, this allows (almost) - * coalesced reading between threads in one line of the tile. - * - * Later on it will be more efficient to re-use the same tile - * and slide it vertically by one line. - * The additionnal need for shared memory will be kfs2 floats and a float8 as register. - * - * Theoritically, it should be possible to handle up to windows-size 83x83 - */ -__kernel void medfilt2d(__global float *image, // input image - __global float *result, // output array - __local float4 *l_data,// local storage 4x the number of threads - int khs1, // Kernel half-size along dim1 (nb lines) - int khs2, // Kernel half-size along dim2 (nb columns) - int height, // Image size along dim1 (nb lines) - int width) // Image size along dim2 (nb columns) -{ - int threadid = get_local_id(0); - //int wg = get_local_size(0); - int x = get_global_id(1); - - if (x < width) - { - union - { - float ary[8]; - float8 vec; - } output, input; - input.vec = (float8)(MAXFLOAT, MAXFLOAT, MAXFLOAT, MAXFLOAT, MAXFLOAT, MAXFLOAT, MAXFLOAT, MAXFLOAT); - int kfs1 = 2 * khs1 + 1; //definition of kernel full size - int kfs2 = 2 * khs2 + 1; - int nbands = (kfs1 + 7) / 8; // 8 elements per thread, aligned vertically in 1 column - for (int y=0; y<height; y++) - { - //Select only the active threads, some may remain inactive - int nb_threads = (nbands * kfs2); - int band_nr = threadid / kfs2; - int band_id = threadid % kfs2; - int pos_x = clamp((int)(x + band_id - khs2), (int) 0, (int) width-1); - int max_vec = clamp(kfs1 - 8 * band_nr, 0, 8); - if (y == 0) - { - for (int i=0; i<max_vec; i++) - { - if (threadid<nb_threads) - { - int pos_y = clamp((int)(y + 8 * band_nr + i - khs1), (int) 0, (int) height-1); - input.ary[i] = image[pos_x + width * pos_y]; - } - } - } - else - { - //store storage.s0 to some shared memory to retrieve it from another thread. - l_data[threadid].s0 = input.vec.s0; - - //Offset to the bottom - input.vec = (float8)(input.vec.s1, - input.vec.s2, - input.vec.s3, - input.vec.s4, - input.vec.s5, - input.vec.s6, - input.vec.s7, - MAXFLOAT); - - barrier(CLK_LOCAL_MEM_FENCE); - - int read_from = threadid + kfs2; - if (read_from < nb_threads) - input.vec.s7 = l_data[read_from].s0; - else if (threadid < nb_threads) //we are on the last band - { - int pos_y = clamp((int)(y + 8 * band_nr + max_vec - 1 - khs1), (int) 0, (int) height-1); - input.ary[max_vec - 1] = image[pos_x + width * pos_y]; - } - - } - - //This function is defined in bitonic.cl - output.vec = my_sort_file(get_local_id(0), get_group_id(0), get_local_size(0), - input.vec, l_data); - - size_t target = (kfs1 * kfs2) / 2; - if (threadid == (target / 8)) - { - result[y * width + x] = output.ary[target % 8]; - } - - } - } -} - |