1 /*
2 * Copyright (c) 2018-2019, 2022 Arm Limited.
3 *
4 * SPDX-License-Identifier: MIT
5 *
6 * Permission is hereby granted, free of charge, to any person obtaining a copy
7 * of this software and associated documentation files (the "Software"), to
8 * deal in the Software without restriction, including without limitation the
9 * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
10 * sell copies of the Software, and to permit persons to whom the Software is
11 * furnished to do so, subject to the following conditions:
12 *
13 * The above copyright notice and this permission notice shall be included in all
14 * copies or substantial portions of the Software.
15 *
16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
17 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
18 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
19 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
20 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
21 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
22 * SOFTWARE.
23 */
24
25 #include "Gather.h"
26
27 #include "arm_compute/core/Types.h"
28 #include "arm_compute/core/utils/misc/ShapeCalculator.h"
29 #include "tests/validation/Helpers.h"
30
31 namespace arm_compute
32 {
33 namespace test
34 {
35 namespace validation
36 {
37 namespace reference
38 {
39 template <typename T>
gather(const SimpleTensor<T> & src,const SimpleTensor<uint32_t> & indices,uint32_t actual_axis)40 SimpleTensor<T> gather(const SimpleTensor<T> &src, const SimpleTensor<uint32_t> &indices, uint32_t actual_axis)
41 {
42 const auto *indices_ptr = static_cast<const uint32_t *>(indices.data());
43 const TensorShape dst_shape = arm_compute::misc::shape_calculator::compute_gather_shape(src.shape(), indices.shape(), actual_axis);
44 SimpleTensor<T> dst(dst_shape, src.data_type());
45
46 Window win;
47 win.use_tensor_dimensions(dst_shape);
48 if(indices.shape().num_dimensions() == 1u)
49 {
50 execute_window_loop(win, [&](const Coordinates & id)
51 {
52 Coordinates offset;
53 for(unsigned int dim = 0; dim < id.num_dimensions(); ++dim)
54 {
55 if(dim == actual_axis)
56 {
57 offset.set(dim, indices_ptr[id[dim]]);
58 }
59 else
60 {
61 offset.set(dim, id[dim]);
62 }
63 }
64 *reinterpret_cast<T *>(dst(id)) = *reinterpret_cast<const T *>(src(offset));
65 });
66 }
67 else
68 {
69 if(actual_axis == 1)
70 {
71 win.set(Window::DimX, Window::Dimension(0, 1, 1));
72 execute_window_loop(win, [&](const Coordinates & id)
73 {
74 auto *dst_ptr = dst(id);
75 Coordinates index_offset;
76 for(uint32_t k = 0; k < indices.shape().num_dimensions(); ++k)
77 {
78 index_offset.set(k, id[k + 1]);
79 }
80 const uint32_t row = *reinterpret_cast<const uint32_t *>(indices(index_offset));
81 Coordinates src_offset;
82 src_offset.set(0, 0);
83 src_offset.set(1, row);
84 for(uint32_t j = 2; j < src.shape().num_dimensions(); ++j)
85 {
86 src_offset.set(j, id[1 + indices.shape().num_dimensions() + (j - 2)]);
87 }
88 const auto in_ptr_row = src(src_offset);
89 memcpy(dst_ptr, in_ptr_row, src.shape()[0] * src.element_size());
90 });
91 }
92 else
93 {
94 ARM_COMPUTE_ERROR("Not implemented.");
95 }
96 }
97
98 return dst;
99 }
100
101 template SimpleTensor<float> gather(const SimpleTensor<float> &src, const SimpleTensor<uint32_t> &indices, uint32_t actual_axis);
102 template SimpleTensor<half> gather(const SimpleTensor<half> &src, const SimpleTensor<uint32_t> &indices, uint32_t actual_axis);
103 template SimpleTensor<uint16_t> gather(const SimpleTensor<uint16_t> &src, const SimpleTensor<uint32_t> &indices, uint32_t actual_axis);
104 template SimpleTensor<uint8_t> gather(const SimpleTensor<uint8_t> &src, const SimpleTensor<uint32_t> &indices, uint32_t actual_axis);
105 } // namespace reference
106 } // namespace validation
107 } // namespace test
108 } // namespace arm_compute
109