From 6ff3b19ee6120edf015fad8caab2991faa3070af Mon Sep 17 00:00:00 2001 From: Anthony Barbier Date: Mon, 4 Sep 2017 18:44:23 +0100 Subject: COMPMID-344 Updated doxygen Change-Id: I32f7b84daa560e460b77216add529c8fa8b327ae --- src/core/NEON/kernels/NEGaussian5x5Kernel.cpp | 203 ++++++++++++++++++++++++++ 1 file changed, 203 insertions(+) create mode 100644 src/core/NEON/kernels/NEGaussian5x5Kernel.cpp (limited to 'src/core/NEON/kernels/NEGaussian5x5Kernel.cpp') diff --git a/src/core/NEON/kernels/NEGaussian5x5Kernel.cpp b/src/core/NEON/kernels/NEGaussian5x5Kernel.cpp new file mode 100644 index 0000000000..f872cc2f0a --- /dev/null +++ b/src/core/NEON/kernels/NEGaussian5x5Kernel.cpp @@ -0,0 +1,203 @@ +/* + * Copyright (c) 2016, 2017 ARM Limited. + * + * SPDX-License-Identifier: MIT + * + * 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. + */ +#include "arm_compute/core/NEON/kernels/NEGaussian5x5Kernel.h" + +#include "arm_compute/core/Coordinates.h" +#include "arm_compute/core/Helpers.h" +#include "arm_compute/core/ITensor.h" +#include "arm_compute/core/NEON/INEKernel.h" +#include "arm_compute/core/Types.h" +#include "arm_compute/core/Validate.h" +#include "arm_compute/core/Window.h" + +#include +#include +#include + +using namespace arm_compute; + +NEGaussian5x5HorKernel::NEGaussian5x5HorKernel() + : _border_size(0) +{ +} + +BorderSize NEGaussian5x5HorKernel::border_size() const +{ + return _border_size; +} + +void NEGaussian5x5HorKernel::configure(const ITensor *input, ITensor *output, bool border_undefined) +{ + ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(input, 1, DataType::U8); + ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(output, 1, DataType::S16); + + _input = input; + _output = output; + _border_size = BorderSize(border_undefined ? 0 : 2, 2); + + // Configure kernel window + constexpr unsigned int num_elems_processed_per_iteration = 8; + constexpr unsigned int num_elems_read_per_iteration = 16; + constexpr unsigned int num_elems_written_per_iteration = 8; + + Window win = calculate_max_window_horizontal(*input->info(), Steps(num_elems_processed_per_iteration), border_undefined, border_size()); + AccessWindowHorizontal output_access(output->info(), 0, num_elems_written_per_iteration); + + update_window_and_padding(win, + AccessWindowHorizontal(input->info(), -border_size().left, num_elems_read_per_iteration), + output_access); + + output_access.set_valid_region(win, input->info()->valid_region(), border_undefined, border_size()); + + INEKernel::configure(win); +} + +void NEGaussian5x5HorKernel::run(const Window &window) +{ + ARM_COMPUTE_ERROR_ON_UNCONFIGURED_KERNEL(this); + ARM_COMPUTE_ERROR_ON_INVALID_SUBWINDOW(INEKernel::window(), window); + + Window win_in(window); + win_in.shift(Window::DimX, -2); + + Iterator input(_input, win_in); + Iterator output(_output, window); + + static const int16x8_t six = vdupq_n_s16(6); + static const int16x8_t four = vdupq_n_s16(4); + + execute_window_loop(window, [&](const Coordinates & id) + { + uint8x16_t data = vld1q_u8(input.ptr()); + + const int16x8x2_t data_s16 = + { + { + vreinterpretq_s16_u16(vmovl_u8(vget_low_u8(data))), + vreinterpretq_s16_u16(vmovl_u8(vget_high_u8(data))) + } + }; + + int16x8_t out = vaddq_s16(data_s16.val[0], vextq_s16(data_s16.val[0], data_s16.val[1], 4)); + out = vmlaq_s16(out, vextq_s16(data_s16.val[0], data_s16.val[1], 1), four); + out = vmlaq_s16(out, vextq_s16(data_s16.val[0], data_s16.val[1], 2), six); + out = vmlaq_s16(out, vextq_s16(data_s16.val[0], data_s16.val[1], 3), four); + + vst1q_s16(reinterpret_cast(output.ptr()), out); + }, + input, output); +} + +BorderSize NEGaussian5x5VertKernel::border_size() const +{ + return BorderSize(2, 0); +} + +void NEGaussian5x5VertKernel::configure(const ITensor *input, ITensor *output, bool border_undefined) +{ + ARM_COMPUTE_ERROR_ON_FORMAT_NOT_IN(input, Format::S16); + ARM_COMPUTE_ERROR_ON_FORMAT_NOT_IN(output, Format::U8); + + _input = input; + _output = output; + + // Configure kernel window + constexpr unsigned int num_elems_processed_per_iteration = 16; + constexpr unsigned int num_elems_read_per_iteration = 32; + constexpr unsigned int num_elems_written_per_iteration = 16; + constexpr unsigned int num_rows_read_per_iteration = 5; + + Window win = calculate_max_window(*input->info(), Steps(num_elems_processed_per_iteration), border_undefined, border_size()); + AccessWindowHorizontal output_access(output->info(), 0, num_elems_written_per_iteration); + + update_window_and_padding(win, + AccessWindowRectangle(input->info(), 0, -border_size().top, num_elems_read_per_iteration, num_rows_read_per_iteration), + output_access); + + output_access.set_valid_region(win, input->info()->valid_region(), border_undefined, border_size()); + + INEKernel::configure(win); +} + +void NEGaussian5x5VertKernel::run(const Window &window) +{ + ARM_COMPUTE_ERROR_ON_UNCONFIGURED_KERNEL(this); + ARM_COMPUTE_ERROR_ON_INVALID_SUBWINDOW(INESimpleKernel::window(), window); + + Iterator input(_input, window); + Iterator output(_output, window); + + const uint8_t *input_top2_ptr = _input->ptr_to_element(Coordinates(0, -2)); + const uint8_t *input_top_ptr = _input->ptr_to_element(Coordinates(0, -1)); + const uint8_t *input_mid_ptr = _input->ptr_to_element(Coordinates(0, 0)); + const uint8_t *input_low_ptr = _input->ptr_to_element(Coordinates(0, 1)); + const uint8_t *input_low2_ptr = _input->ptr_to_element(Coordinates(0, 2)); + + const uint16x8_t six = vdupq_n_u16(6); + const uint16x8_t four = vdupq_n_u16(4); + + execute_window_loop(window, [&](const Coordinates & id) + { + const size_t input_offset_high_s16 = input.offset(); + const size_t input_offset_low_s16 = input.offset() + 16; + + //HIGH DATA + //top2 + uint16x8_t data_high = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_top2_ptr + input_offset_high_s16))); + uint16x8_t out_high = data_high; + //top + data_high = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_top_ptr + input_offset_high_s16))); + out_high = vmlaq_u16(out_high, data_high, four); + //mid + data_high = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_mid_ptr + input_offset_high_s16))); + out_high = vmlaq_u16(out_high, data_high, six); + //low + data_high = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_low_ptr + input_offset_high_s16))); + out_high = vmlaq_u16(out_high, data_high, four); + //low2 + data_high = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_low2_ptr + input_offset_high_s16))); + out_high = vaddq_u16(out_high, data_high); + + //LOW DATA + //top2 + uint16x8_t data_low = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_top2_ptr + input_offset_low_s16))); + uint16x8_t out_low = data_low; + //top + data_low = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_top_ptr + input_offset_low_s16))); + out_low = vmlaq_u16(out_low, data_low, four); + //mid + data_low = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_mid_ptr + input_offset_low_s16))); + out_low = vmlaq_u16(out_low, data_low, six); + //low + data_low = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_low_ptr + input_offset_low_s16))); + out_low = vmlaq_u16(out_low, data_low, four); + //low2 + data_low = vreinterpretq_u16_s16(vld1q_s16(reinterpret_cast(input_low2_ptr + input_offset_low_s16))); + out_low = vaddq_u16(out_low, data_low); + + vst1q_u8(output.ptr(), vcombine_u8(vqshrn_n_u16(out_high, 8), + vqshrn_n_u16(out_low, 8))); + }, + input, output); +} -- cgit v1.2.1