blob: aa458c21195f6e98cceab039011e4d5729f39bc5 [file] [log] [blame]
giuros0192fd9432018-12-03 17:30:00 +00001/*
George Wortd88590f2018-12-12 17:39:58 +00002 * Copyright (c) 2018-2019 ARM Limited.
giuros0192fd9432018-12-03 17:30:00 +00003 *
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#include "arm_compute/core/NEON/kernels/NEElementwiseOperationKernel.h"
25
26#include "arm_compute/core/CPP/Validate.h"
27#include "arm_compute/core/Error.h"
28#include "arm_compute/core/Helpers.h"
29#include "arm_compute/core/IAccessWindow.h"
30#include "arm_compute/core/ITensor.h"
31#include "arm_compute/core/NEON/NEAsymm.h"
32#include "arm_compute/core/NEON/NEFixedPoint.h"
33#include "arm_compute/core/NEON/wrapper/wrapper.h"
34#include "arm_compute/core/TensorInfo.h"
35#include "arm_compute/core/Validate.h"
36
37#include <algorithm>
38#include <arm_neon.h>
39#include <cstdint>
40#include <map>
41#include <string>
42
43namespace arm_compute
44{
45class Coordinates;
46
47namespace
48{
49float32x4x4_t load_quantized(const uint8_t *input1_ptr, const int32x4_t &offset, const float32x4_t &scale)
50{
51 qasymm8x16_t x = vld1q_u8(input1_ptr);
52 const float32x4x4_t out =
53 {
54 {
55 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_low_u16(vmovl_u8(vget_low_u8(x))))), offset)), scale),
56 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_high_u16(vmovl_u8(vget_low_u8(x))))), offset)), scale),
57 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_low_u16(vmovl_u8(vget_high_u8(x))))), offset)), scale),
58 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_high_u16(vmovl_u8(vget_high_u8(x))))), offset)), scale),
59 }
60 };
61 return out;
62}
63
George Wortd88590f2018-12-12 17:39:58 +000064void store_quantized(uint8_t *output_ptr, const uint32x4x4_t &out)
65{
66 const uint8x8_t pa = vqmovn_u16(vcombine_u16(vqmovn_u32(out.val[0]), vqmovn_u32(out.val[1])));
67 const uint8x8_t pb = vqmovn_u16(vcombine_u16(vqmovn_u32(out.val[2]), vqmovn_u32(out.val[3])));
68 vst1q_u8(output_ptr, vcombine_u8(pa, pb));
69}
70
71void store_quantized(uint8_t *output_ptr, const int32x4x4_t &out)
72{
73 const uint8x8_t pa = vqmovun_s16(vcombine_s16(vqmovn_s32(out.val[0]), vqmovn_s32(out.val[1])));
74 const uint8x8_t pb = vqmovun_s16(vcombine_s16(vqmovn_s32(out.val[2]), vqmovn_s32(out.val[3])));
75 vst1q_u8(output_ptr, vcombine_u8(pa, pb));
76}
77
giuros0192fd9432018-12-03 17:30:00 +000078void store_quantized(uint8_t *output_ptr, const float32x4x4_t &rf, const float32x4_t &offset, const float32x4_t &invscale)
79{
80 int32x4x4_t out =
81 {
Georgios Pinitasd57891a2019-02-19 18:10:03 +000082 {
83 vcvtq_s32_f32(vmlaq_f32(offset, rf.val[0], invscale)),
84 vcvtq_s32_f32(vmlaq_f32(offset, rf.val[1], invscale)),
85 vcvtq_s32_f32(vmlaq_f32(offset, rf.val[2], invscale)),
86 vcvtq_s32_f32(vmlaq_f32(offset, rf.val[3], invscale)),
87 }
giuros0192fd9432018-12-03 17:30:00 +000088 };
George Wortd88590f2018-12-12 17:39:58 +000089 store_quantized(output_ptr, out);
giuros0192fd9432018-12-03 17:30:00 +000090}
91
92float32x4x4_t dup_quantized(qasymm8_t broadcast_value, int offset, float scale)
93{
94 const qasymm8x16_t broadcast_value_vec = vdupq_n_u8(broadcast_value);
95 const int32x4_t voffset = vdupq_n_s32(offset);
96 const float32x4_t vscale = vdupq_n_f32(scale);
97
98 const float32x4x4_t broadcast_vector =
99 {
100 {
101 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_low_u16(vmovl_u8(vget_low_u8(broadcast_value_vec))))), voffset)), vscale),
102 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_high_u16(vmovl_u8(vget_low_u8(broadcast_value_vec))))), voffset)), vscale),
103 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_low_u16(vmovl_u8(vget_high_u8(broadcast_value_vec))))), voffset)), vscale),
104 vmulq_f32(vcvtq_f32_s32(vsubq_s32(vreinterpretq_s32_u32(vmovl_u16(vget_high_u16(vmovl_u8(vget_high_u8(broadcast_value_vec))))), voffset)), vscale),
105 }
106 };
107 return broadcast_vector;
108}
109
110template <ArithmeticOperation op, typename ScalarType>
George Wortd88590f2018-12-12 17:39:58 +0000111inline ScalarType elementwise_arithm_op_scalar(const ScalarType &a, const ScalarType &b)
giuros0192fd9432018-12-03 17:30:00 +0000112{
113 auto res = ScalarType(0);
114
115 switch(op)
116 {
117 case ArithmeticOperation::MAX:
118 res = std::max(a, b);
119 break;
120 case ArithmeticOperation::MIN:
121 res = std::min(a, b);
122 break;
123 case ArithmeticOperation::SQUARED_DIFF:
124 {
125 res = (a - b) * (a - b);
126 break;
127 }
George Worta1e7e282019-01-15 11:00:29 +0000128 case ArithmeticOperation::DIV:
129 {
130 res = a / b;
131 break;
132 }
giuros0192fd9432018-12-03 17:30:00 +0000133 default:
134 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
135 }
136 return res;
137}
138
George Wortd88590f2018-12-12 17:39:58 +0000139template <ArithmeticOperation op>
140inline uint8_t elementwise_arithm_op_quantized_scalar(const float &a, const float &b, QuantizationInfo qinfo)
141{
142 return qinfo.quantize(elementwise_arithm_op_scalar<op>(a, b), RoundingPolicy::TO_NEAREST_UP);
143}
144
giuros0192fd9432018-12-03 17:30:00 +0000145template <ArithmeticOperation op, typename VectorType>
George Wortd88590f2018-12-12 17:39:58 +0000146inline VectorType elementwise_arithm_op(const VectorType &a, const VectorType &b)
giuros0192fd9432018-12-03 17:30:00 +0000147{
148 VectorType res = { 0, 0, 0, 0 };
149
150 switch(op)
151 {
152 case ArithmeticOperation::MAX:
153 res = wrapper::vmax(a, b);
154 break;
155 case ArithmeticOperation::MIN:
156 res = wrapper::vmin(a, b);
157 break;
158 case ArithmeticOperation::SQUARED_DIFF:
159 {
160 const VectorType tmp = wrapper::vsub(a, b);
161 res = wrapper::vmul(tmp, tmp);
162 break;
163 }
giuros0192fd9432018-12-03 17:30:00 +0000164 default:
165 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
166 }
167
168 return res;
169}
170
George Worta1e7e282019-01-15 11:00:29 +0000171template <>
172inline float32x4_t elementwise_arithm_op<ArithmeticOperation::DIV, float32x4_t>(const float32x4_t &a, const float32x4_t &b)
173{
174 return wrapper::vdiv(a, b);
175}
176
177#ifdef __ARM_FEATURE_FP16_VECTOR_ARITHMETIC
178template <>
179inline float16x8_t elementwise_arithm_op<ArithmeticOperation::DIV, float16x8_t>(const float16x8_t &a, const float16x8_t &b)
180{
181 return wrapper::vdiv(a, b);
182}
183#endif // __ARM_FEATURE_FP16_VECTOR_ARITHMETIC
184
giuros0192fd9432018-12-03 17:30:00 +0000185template <ArithmeticOperation op>
George Wortd88590f2018-12-12 17:39:58 +0000186inline float32x4x4_t elementwise_arithm_op(const float32x4x4_t &a, const float32x4x4_t &b)
giuros0192fd9432018-12-03 17:30:00 +0000187{
188 float32x4x4_t out =
189 {
Georgios Pinitasd57891a2019-02-19 18:10:03 +0000190 {
191 elementwise_arithm_op<op>(a.val[0], b.val[0]),
192 elementwise_arithm_op<op>(a.val[1], b.val[1]),
193 elementwise_arithm_op<op>(a.val[2], b.val[2]),
194 elementwise_arithm_op<op>(a.val[3], b.val[3]),
195 }
giuros0192fd9432018-12-03 17:30:00 +0000196 };
197 return out;
198}
199
George Wortd88590f2018-12-12 17:39:58 +0000200template <ArithmeticOperation op, typename ScalarType, typename VectorType>
201inline VectorType elementwise_arithm_op_broadcast(const VectorType &a, const ScalarType &broadcast_value, const bool reorder)
202{
203 VectorType broadcast_vector = wrapper::vdup_n(broadcast_value, wrapper::traits::vector_128_tag());
204 return elementwise_arithm_op<op>(reorder ? broadcast_vector : a, reorder ? a : broadcast_vector);
205}
206
207template <ComparisonOperation op, typename InputScalarType>
208inline uint8_t elementwise_comp_op_scalar(const InputScalarType &a, const InputScalarType &b)
209{
210 bool res = false;
211
212 switch(op)
213 {
214 case ComparisonOperation::Equal:
215 res = (a == b);
216 break;
217 case ComparisonOperation::NotEqual:
218 res = (a != b);
219 break;
220 case ComparisonOperation::Greater:
221 res = (a > b);
222 break;
223 case ComparisonOperation::GreaterEqual:
224 res = (a >= b);
225 break;
226 case ComparisonOperation::Less:
227 res = (a < b);
228 break;
229 case ComparisonOperation::LessEqual:
230 res = (a <= b);
231 break;
232 default:
233 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
234 }
235 return res ? ~static_cast<uint8_t>(0) : static_cast<uint8_t>(0);
236}
237
238template <ComparisonOperation op>
239inline uint8_t elementwise_comp_op_quantized_scalar(const float &a, const float &b, QuantizationInfo qinfo)
240{
241 ARM_COMPUTE_UNUSED(qinfo);
242 return elementwise_comp_op_scalar<op>(a, b);
243}
244
245template <ComparisonOperation op, typename InputVectorType, typename OutputVectorType>
246inline OutputVectorType elementwise_comp_op(const InputVectorType &a, const InputVectorType &b)
247{
248 OutputVectorType res = { 0, 0, 0, 0 };
249
250 switch(op)
251 {
252 case ComparisonOperation::Equal:
253 res = wrapper::vceq(a, b);
254 break;
255 case ComparisonOperation::NotEqual:
256 res = wrapper::vnot(wrapper::vceq(a, b));
257 break;
258 case ComparisonOperation::Greater:
259 res = wrapper::vcgt(a, b);
260 break;
261 case ComparisonOperation::GreaterEqual:
262 res = wrapper::vcge(a, b);
263 break;
264 case ComparisonOperation::Less:
265 res = wrapper::vcgt(b, a);
266 break;
267 case ComparisonOperation::LessEqual:
268 res = wrapper::vcge(b, a);
269 break;
270 default:
271 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
272 }
273
274 return res;
275}
276
277template <ComparisonOperation op>
278inline uint32x4x4_t elementwise_comp_op(const float32x4x4_t &a, const float32x4x4_t &b)
279{
280 uint32x4x4_t out =
281 {
Georgios Pinitasd57891a2019-02-19 18:10:03 +0000282 {
283 elementwise_comp_op<op, float32x4_t, uint32x4_t>(a.val[0], b.val[0]),
284 elementwise_comp_op<op, float32x4_t, uint32x4_t>(a.val[1], b.val[1]),
285 elementwise_comp_op<op, float32x4_t, uint32x4_t>(a.val[2], b.val[2]),
286 elementwise_comp_op<op, float32x4_t, uint32x4_t>(a.val[3], b.val[3])
287 }
George Wortd88590f2018-12-12 17:39:58 +0000288 };
289 return out;
290}
291
292template <ComparisonOperation op, typename InputScalarType, typename InputVectorType, typename OutputVectorType>
293inline OutputVectorType elementwise_comp_op_broadcast(const InputVectorType &a, const InputScalarType &broadcast_value, const bool reorder)
294{
295 InputVectorType broadcast_vector = wrapper::vdup_n(broadcast_value, wrapper::traits::vector_128_tag());
296 return elementwise_comp_op<op, InputVectorType, OutputVectorType>(reorder ? broadcast_vector : a, reorder ? a : broadcast_vector);
297}
298
299template <ArithmeticOperation op, typename ScalarType, typename VectorType>
300inline int elementwise_arithm_op_loop(int window_start_x, int window_end_x, int window_step_x,
301 const ScalarType *input1_ptr, const ScalarType *input2_ptr, ScalarType *output_ptr)
302{
303 int x = window_start_x;
304 for(; x <= (window_end_x - window_step_x); x += window_step_x)
305 {
306 const auto a = wrapper::vloadq(input1_ptr + x);
307 const auto b = wrapper::vloadq(input2_ptr + x);
308 wrapper::vstore(output_ptr + x, elementwise_arithm_op<op>(a, b));
309 }
310 return x;
311}
312
313template <ArithmeticOperation op>
314inline int elementwise_arithm_op_quantized_loop(int window_start_x, int window_end_x, int window_step_x,
315 const uint8_t *input1_ptr, const uint8_t *input2_ptr, uint8_t *output_ptr,
316 int32x4_t voffset1, int32x4_t voffset2, float32x4_t vscale1, float32x4_t vscale2,
317 float32x4_t voffseto, float32x4_t invvscaleo)
318{
319 int x = window_start_x;
320 for(; x <= (window_end_x - window_step_x); x += window_step_x)
321 {
322 // Get inputs and compute output
323 const float32x4x4_t af = load_quantized(input1_ptr + x, voffset1, vscale1);
324 const float32x4x4_t bf = load_quantized(input2_ptr + x, voffset2, vscale2);
325 const float32x4x4_t rf = elementwise_arithm_op<op>(af, bf);
326 store_quantized(output_ptr + x, rf, voffseto, invvscaleo);
327 }
328 return x;
329}
330
331template <ArithmeticOperation op, typename ScalarType, typename VectorType>
332inline int elementwise_arithm_op_broadcast_loop(int window_start_x, int window_end_x, int window_step_x,
333 const ScalarType *non_broadcast_input_ptr, const ScalarType &broadcast_value, ScalarType *output_ptr, const bool reorder)
334{
335 int x = window_start_x;
336 for(; x <= (window_end_x - window_step_x); x += window_step_x)
337 {
338 const auto a = wrapper::vloadq((non_broadcast_input_ptr + x));
339 wrapper::vstore(output_ptr + x, elementwise_arithm_op_broadcast<op>(a, broadcast_value, reorder));
340 }
341 return x;
342}
343
344template <ArithmeticOperation op>
345inline int elementwise_arithm_op_quantized_broadcast_loop(int window_start_x, int window_end_x, int window_step_x,
346 const uint8_t *non_broadcast_input_ptr, float32x4x4_t broadcast_vector, uint8_t *output_ptr,
347 int32x4_t voffset_non_broadcast, float32x4_t vscale_non_broadcast,
348 float32x4_t voffseto, float32x4_t invvscaleo, bool reorder)
349{
350 int x = window_start_x;
351 for(; x <= (window_end_x - window_step_x); x += window_step_x)
352 {
353 const float32x4x4_t af = load_quantized(non_broadcast_input_ptr + x, voffset_non_broadcast, vscale_non_broadcast);
354 const float32x4x4_t rf = elementwise_arithm_op<op>(reorder ? broadcast_vector : af, reorder ? af : broadcast_vector);
355 store_quantized(output_ptr + x, rf, voffseto, invvscaleo);
356 }
357 return x;
358}
359
360template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
361inline int elementwise_comp_op_16_loop(int window_start_x, int window_end_x, int window_step_x,
362 const InputScalarType *input1_ptr, const InputScalarType *input2_ptr, uint8_t *output_ptr)
363{
364 int x = window_start_x;
365 for(; x <= (window_end_x - window_step_x); x += window_step_x)
366 {
367 const auto a = wrapper::vloadq(input1_ptr + x);
368 const auto b = wrapper::vloadq(input2_ptr + x);
369 const auto res = elementwise_comp_op<op, InputVectorType, uint16x8_t>(a, b);
370 wrapper::vstore(output_ptr + x, wrapper::vmovn(res));
371 }
372 return x;
373}
374
375template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
376inline int elementwise_comp_op_32_loop(int window_start_x, int window_end_x, int window_step_x,
377 const InputScalarType *input1_ptr, const InputScalarType *input2_ptr, uint8_t *output_ptr)
378{
379 int x = window_start_x;
380 for(; x <= (window_end_x - window_step_x); x += window_step_x)
381 {
382 auto a = wrapper::vloadq(input1_ptr + x);
383 auto b = wrapper::vloadq(input2_ptr + x);
384 const auto res = elementwise_comp_op<op, InputVectorType, uint32x4_t>(a, b);
385 a = wrapper::vloadq(input1_ptr + x + 4);
386 b = wrapper::vloadq(input2_ptr + x + 4);
387 const auto res2 = elementwise_comp_op<op, InputVectorType, uint32x4_t>(a, b);
388 wrapper::vstore(output_ptr + x, wrapper::vmovn(wrapper::vcombine(wrapper::vmovn(res), wrapper::vmovn(res2))));
389 }
390 if(x <= window_end_x - 4)
391 {
392 const auto a = wrapper::vloadq(input1_ptr + x);
393 const auto b = wrapper::vloadq(input2_ptr + x);
394 const auto res = elementwise_comp_op<op, InputVectorType, uint32x4_t>(a, b);
395 for(int i = 0; i < 4; i++)
396 {
397 *(output_ptr + x + i) = wrapper::vgetlane(res, i);
398 }
399 x = +4;
400 }
401 return x;
402}
403
404template <ComparisonOperation op>
405inline int elementwise_comp_op_quantized_loop(int window_start_x, int window_end_x, int window_step_x,
406 const uint8_t *input1_ptr, const uint8_t *input2_ptr, uint8_t *output_ptr,
407 int32x4_t voffset1, int32x4_t voffset2, float32x4_t vscale1, float32x4_t vscale2,
408 float32x4_t voffseto, float32x4_t invvscaleo)
409{
410 ARM_COMPUTE_UNUSED(voffseto, invvscaleo);
411 int x = window_start_x;
412 for(; x <= (window_end_x - window_step_x); x += window_step_x)
413 {
414 const float32x4x4_t af = load_quantized(input1_ptr + x, voffset1, vscale1);
415 const float32x4x4_t bf = load_quantized(input2_ptr + x, voffset2, vscale2);
416 const uint32x4x4_t rf = elementwise_comp_op<op>(af, bf);
417 store_quantized(output_ptr + x, rf);
418 }
419 return x;
420}
421
422template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
423inline int elementwise_comp_op_broadcast_16_loop(int window_start_x, int window_end_x, int window_step_x,
424 const InputScalarType *non_broadcast_input_ptr, const InputScalarType &broadcast_value, uint8_t *output_ptr, const bool reorder)
425{
426 int x = window_start_x;
427 for(; x <= (window_end_x - window_step_x); x += window_step_x)
428 {
429 const auto a = elementwise_comp_op_broadcast<op, InputScalarType, InputVectorType, uint16x8_t>(wrapper::vloadq((non_broadcast_input_ptr + x)), broadcast_value, reorder);
430 wrapper::vstore(output_ptr + x, wrapper::vmovn(a));
431 }
432 return x;
433}
434
435template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
436inline int elementwise_comp_op_broadcast_32_loop(int window_start_x, int window_end_x, int window_step_x,
437 const InputScalarType *non_broadcast_input_ptr, const InputScalarType &broadcast_value, uint8_t *output_ptr, const bool reorder)
438{
439 int x = window_start_x;
440 for(; x <= (window_end_x - window_step_x); x += window_step_x)
441 {
442 const auto a = elementwise_comp_op_broadcast<op, InputScalarType, InputVectorType, uint32x4_t>(wrapper::vloadq(non_broadcast_input_ptr + x), broadcast_value, reorder);
443 const auto b = elementwise_comp_op_broadcast<op, InputScalarType, InputVectorType, uint32x4_t>(wrapper::vloadq(non_broadcast_input_ptr + x + 4), broadcast_value, reorder);
444 wrapper::vstore(output_ptr + x, wrapper::vmovn(wrapper::vcombine(wrapper::vmovn(a), wrapper::vmovn(b))));
445 }
446 if(x <= window_end_x - 4)
447 {
448 const auto a = elementwise_comp_op_broadcast<op, InputScalarType, InputVectorType, uint32x4_t>(wrapper::vloadq((non_broadcast_input_ptr + x)), broadcast_value, reorder);
449 for(int i = 0; i < 4; i++)
450 {
451 *(output_ptr + x + i) = wrapper::vgetlane(a, i);
452 }
453 x = +4;
454 }
455 return x;
456}
457
458template <ComparisonOperation op>
459inline int elementwise_comp_op_quantized_broadcast_loop(int window_start_x, int window_end_x, int window_step_x,
460 const uint8_t *non_broadcast_input_ptr, float32x4x4_t broadcast_vector, uint8_t *output_ptr,
461 int32x4_t voffset_non_broadcast, float32x4_t vscale_non_broadcast,
462 float32x4_t voffseto, float32x4_t invvscaleo, bool reorder)
463{
464 ARM_COMPUTE_UNUSED(voffseto, invvscaleo);
465 int x = window_start_x;
466 for(; x <= (window_end_x - window_step_x); x += window_step_x)
467 {
468 const float32x4x4_t af = load_quantized(non_broadcast_input_ptr + x, voffset_non_broadcast, vscale_non_broadcast);
469 const uint32x4x4_t rf = elementwise_comp_op<op>(reorder ? broadcast_vector : af, reorder ? af : broadcast_vector);
470 store_quantized(output_ptr + x, rf);
471 }
472 return x;
473}
474
475template <typename InputScalarType, typename OutputScalarType, typename InputVectorType>
476void elementwise_op(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window,
477 OutputScalarType (*scalar_func)(const InputScalarType &, const InputScalarType &),
478 int (*broadcast_func)(int, int, int, const InputScalarType *, const InputScalarType &, OutputScalarType *, const bool),
479 int (*neon_func)(int, int, int, const InputScalarType *, const InputScalarType *, OutputScalarType *))
giuros0192fd9432018-12-03 17:30:00 +0000480{
481 // Create input windows
482 Window input1_win = window.broadcast_if_dimension_le_one(in1->info()->tensor_shape());
483 Window input2_win = window.broadcast_if_dimension_le_one(in2->info()->tensor_shape());
484
485 // Clear X Dimension on execution window as we handle manually
486 Window win = window;
487 win.set(Window::DimX, Window::Dimension(0, 1, 1));
488
Michalis Spyroue8c0c432019-01-22 11:08:31 +0000489 const int window_step_x = std::min(16 / static_cast<int>(sizeof(OutputScalarType)), 8);
giuros0192fd9432018-12-03 17:30:00 +0000490 const auto window_start_x = static_cast<int>(window.x().start());
491 const auto window_end_x = static_cast<int>(window.x().end());
492 const bool is_broadcast_across_x = (input1_win.x().step() == 0) || (input2_win.x().step() == 0);
493
494 if(is_broadcast_across_x)
495 {
giuros0192fd9432018-12-03 17:30:00 +0000496 const bool is_broadcast_input_2 = input2_win.x().step() == 0;
497 Window broadcast_win = is_broadcast_input_2 ? input2_win : input1_win;
498 Window non_broadcast_win = !is_broadcast_input_2 ? input2_win : input1_win;
499 const ITensor *broadcast_tensor = is_broadcast_input_2 ? in2 : in1;
500 const ITensor *non_broadcast_tensor = !is_broadcast_input_2 ? in2 : in1;
501
502 // Clear X Dimension on execution window as we handle manually
503 non_broadcast_win.set(Window::DimX, Window::Dimension(0, 1, 1));
504
505 Iterator broadcast_input(broadcast_tensor, broadcast_win);
506 Iterator non_broadcast_input(non_broadcast_tensor, non_broadcast_win);
507 Iterator output(out, win);
508
509 execute_window_loop(win, [&](const Coordinates & id)
510 {
George Wortd88590f2018-12-12 17:39:58 +0000511 auto output_ptr = reinterpret_cast<OutputScalarType *>(output.ptr());
512 const auto non_broadcast_input_ptr = reinterpret_cast<const InputScalarType *>(non_broadcast_input.ptr());
513 const InputScalarType broadcast_value = *reinterpret_cast<const InputScalarType *>(broadcast_input.ptr());
giuros0192fd9432018-12-03 17:30:00 +0000514
George Wortd88590f2018-12-12 17:39:58 +0000515 int x = (*broadcast_func)(window_start_x, window_end_x, window_step_x, non_broadcast_input_ptr, broadcast_value, output_ptr, !is_broadcast_input_2);
giuros0192fd9432018-12-03 17:30:00 +0000516 for(; x < window_end_x; ++x)
517 {
518 const auto a = *(non_broadcast_input_ptr + x);
George Wortd88590f2018-12-12 17:39:58 +0000519 *(output_ptr + x) = (*scalar_func)(!is_broadcast_input_2 ? broadcast_value : a, !is_broadcast_input_2 ? a : broadcast_value);
giuros0192fd9432018-12-03 17:30:00 +0000520 }
521 },
522 broadcast_input, non_broadcast_input, output);
523 }
524 else
525 {
526 // Clear X Dimension on execution window as we handle manually
527 input1_win.set(Window::DimX, Window::Dimension(0, 1, 1));
528 input2_win.set(Window::DimX, Window::Dimension(0, 1, 1));
529
530 Iterator input1(in1, input1_win);
531 Iterator input2(in2, input2_win);
532 Iterator output(out, win);
533
534 execute_window_loop(win, [&](const Coordinates & id)
535 {
George Wortd88590f2018-12-12 17:39:58 +0000536 auto output_ptr = reinterpret_cast<OutputScalarType *>(output.ptr());
537 const auto input1_ptr = reinterpret_cast<const InputScalarType *>(input1.ptr());
538 const auto input2_ptr = reinterpret_cast<const InputScalarType *>(input2.ptr());
giuros0192fd9432018-12-03 17:30:00 +0000539
George Wortd88590f2018-12-12 17:39:58 +0000540 int x = (*neon_func)(window_start_x, window_end_x, window_step_x, input1_ptr, input2_ptr, output_ptr);
giuros0192fd9432018-12-03 17:30:00 +0000541 for(; x < window_end_x; ++x)
542 {
543 const auto a = *(input1_ptr + x);
544 const auto b = *(input2_ptr + x);
George Wortd88590f2018-12-12 17:39:58 +0000545 *(output_ptr + x) = (*scalar_func)(a, b);
giuros0192fd9432018-12-03 17:30:00 +0000546 }
giuros0192fd9432018-12-03 17:30:00 +0000547 },
548 input1, input2, output);
549 }
550}
551
George Wortd88590f2018-12-12 17:39:58 +0000552void elementwise_op_quantized(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window,
553 uint8_t (*scalar_func)(const float &, const float &, QuantizationInfo),
554 int (*broadcast_func)(int, int, int, const uint8_t *, float32x4x4_t, uint8_t *, int32x4_t, float32x4_t,
555 float32x4_t, float32x4_t, const bool),
556 int (*neon_func)(int, int, int, const uint8_t *, const uint8_t *, uint8_t *,
557 int32x4_t, int32x4_t, float32x4_t, float32x4_t,
558 float32x4_t, float32x4_t))
giuros0192fd9432018-12-03 17:30:00 +0000559{
560 // Create input windows
561 Window input1_win = window.broadcast_if_dimension_le_one(in1->info()->tensor_shape());
562 Window input2_win = window.broadcast_if_dimension_le_one(in2->info()->tensor_shape());
563
564 // Clear X Dimension on execution window as we handle manually
565 Window win = window;
566 win.set(Window::DimX, Window::Dimension(0, 1, 1));
567
568 const int window_step_x = 16;
569 const auto window_start_x = static_cast<int>(window.x().start());
570 const auto window_end_x = static_cast<int>(window.x().end());
571 const bool is_broadcast_across_x = (input1_win.x().step() == 0) || (input2_win.x().step() == 0);
572
573 const float output_scale = out->info()->quantization_info().scale;
574 const int output_offset = out->info()->quantization_info().offset;
575
576 // Output quantization info (add 0.5 to round toward the nearest integer - 0.5 rounds away from zero)
577 const float32x4_t voffseto = vdupq_n_f32(output_offset + 0.5f);
578 const float32x4_t invvscaleo = vdupq_n_f32(1.f / output_scale);
579
580 if(is_broadcast_across_x)
581 {
582 // Select the broadcast input on the X axis
583 const bool is_broadcast_input_2 = input2_win.x().step() == 0;
584 Window broadcast_win = is_broadcast_input_2 ? input2_win : input1_win;
585 Window non_broadcast_win = !is_broadcast_input_2 ? input2_win : input1_win;
586 const ITensor *broadcast_tensor = is_broadcast_input_2 ? in2 : in1;
587 const ITensor *non_broadcast_tensor = !is_broadcast_input_2 ? in2 : in1;
588
589 const QuantizationInfo broadcast_qinfo = broadcast_tensor->info()->quantization_info();
590 const QuantizationInfo non_broadcast_qinfo = non_broadcast_tensor->info()->quantization_info();
591
592 const int32x4_t voffset_non_broadcast = vdupq_n_s32(non_broadcast_qinfo.offset);
593 const float32x4_t vscale_non_broadcast = vdupq_n_f32(non_broadcast_qinfo.scale);
594
595 // Clear X Dimension on execution window as we handle manually
596 non_broadcast_win.set(Window::DimX, Window::Dimension(0, 1, 1));
597
598 Iterator broadcast_input(broadcast_tensor, broadcast_win);
599 Iterator non_broadcast_input(non_broadcast_tensor, non_broadcast_win);
600 Iterator output(out, win);
601
602 execute_window_loop(win, [&](const Coordinates & id)
603 {
604 const auto non_broadcast_input_ptr = reinterpret_cast<const uint8_t *>(non_broadcast_input.ptr());
605 const auto output_ptr = reinterpret_cast<uint8_t *>(output.ptr());
606
607 const uint8_t broadcast_value = *reinterpret_cast<const uint8_t *>(broadcast_input.ptr());
608 const float32x4x4_t broadcast_vector = dup_quantized(broadcast_value, broadcast_qinfo.offset, broadcast_qinfo.scale);
609
George Wortd88590f2018-12-12 17:39:58 +0000610 int x = (*broadcast_func)(window_start_x, window_end_x, window_step_x, non_broadcast_input_ptr, broadcast_vector, output_ptr,
611 voffset_non_broadcast, vscale_non_broadcast, voffseto, invvscaleo, !is_broadcast_input_2);
giuros0192fd9432018-12-03 17:30:00 +0000612 for(; x < window_end_x; ++x)
613 {
George Wortd88590f2018-12-12 17:39:58 +0000614 const float afs = scvt_f32_qasymm8(*(non_broadcast_input_ptr + x), non_broadcast_qinfo.scale, non_broadcast_qinfo.offset);
615 const float bfs = scvt_f32_qasymm8(broadcast_value, broadcast_qinfo.scale, broadcast_qinfo.offset);
616 *(output_ptr + x) = (*scalar_func)(!is_broadcast_input_2 ? bfs : afs, !is_broadcast_input_2 ? afs : bfs,
617 out->info()->quantization_info());
giuros0192fd9432018-12-03 17:30:00 +0000618 }
619 },
620 broadcast_input, non_broadcast_input, output);
621 }
622 else
623 {
624 // Input1 quantization info
625 const int32x4_t voffset1 = vdupq_n_s32(in1->info()->quantization_info().offset);
626 const float32x4_t vscale1 = vdupq_n_f32(in1->info()->quantization_info().scale);
627
628 // Input2 quantization info
629 const int32x4_t voffset2 = vdupq_n_s32(in2->info()->quantization_info().offset);
630 const float32x4_t vscale2 = vdupq_n_f32(in2->info()->quantization_info().scale);
631
632 // Clear X Dimension on execution window as we handle manually
633 input1_win.set(Window::DimX, Window::Dimension(0, 1, 1));
634 input2_win.set(Window::DimX, Window::Dimension(0, 1, 1));
635
636 const QuantizationInfo input1_qinfo = in1->info()->quantization_info();
637 const QuantizationInfo input2_qinfo = in2->info()->quantization_info();
638
639 Iterator input1(in1, input1_win);
640 Iterator input2(in2, input2_win);
641 Iterator output(out, win);
642
643 execute_window_loop(win, [&](const Coordinates & id)
644 {
645 const auto input1_ptr = reinterpret_cast<const uint8_t *>(input1.ptr());
646 const auto input2_ptr = reinterpret_cast<const uint8_t *>(input2.ptr());
647 const auto output_ptr = reinterpret_cast<uint8_t *>(output.ptr());
648
George Wortd88590f2018-12-12 17:39:58 +0000649 int x = (*neon_func)(window_start_x, window_end_x, window_step_x, input1_ptr, input2_ptr, output_ptr, voffset1, voffset2,
650 vscale1, vscale2, voffseto, invvscaleo);
giuros0192fd9432018-12-03 17:30:00 +0000651 for(; x < window_end_x; ++x)
652 {
George Wortd88590f2018-12-12 17:39:58 +0000653 const float afs = scvt_f32_qasymm8(*(input1_ptr + x), input1_qinfo.scale, input1_qinfo.offset);
654 const float bfs = scvt_f32_qasymm8(*(input2_ptr + x), input2_qinfo.scale, input2_qinfo.offset);
655 *(output_ptr + x) = (*scalar_func)(afs, bfs, out->info()->quantization_info());
giuros0192fd9432018-12-03 17:30:00 +0000656 }
657 },
658 input1, input2, output);
659 }
660}
661
George Wortd88590f2018-12-12 17:39:58 +0000662template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
663void elementwise_comp_op_16(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window)
giuros0192fd9432018-12-03 17:30:00 +0000664{
George Wortd88590f2018-12-12 17:39:58 +0000665 elementwise_op<InputScalarType, uint8_t, InputVectorType>(in1, in2, out, window,
666 &elementwise_comp_op_scalar<op, InputScalarType>,
667 &elementwise_comp_op_broadcast_16_loop<op, InputScalarType, InputVectorType>,
668 &elementwise_comp_op_16_loop<op, InputScalarType, InputVectorType>);
669}
670
671template <ComparisonOperation op, typename InputScalarType, typename InputVectorType>
672void elementwise_comp_op_32(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window)
673{
674 elementwise_op<InputScalarType, uint8_t, InputVectorType>(in1, in2, out, window,
675 &elementwise_comp_op_scalar<op, InputScalarType>,
676 &elementwise_comp_op_broadcast_32_loop<op, InputScalarType, InputVectorType>,
677 &elementwise_comp_op_32_loop<op, InputScalarType, InputVectorType>);
678}
679
680template <ArithmeticOperation op, typename ScalarType, typename VectorType>
681void elementwise_arithm_op(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window)
682{
683 elementwise_op<ScalarType, ScalarType, VectorType>(in1, in2, out, window,
684 &elementwise_arithm_op_scalar<op, ScalarType>,
685 &elementwise_arithm_op_broadcast_loop<op, ScalarType, VectorType>,
686 &elementwise_arithm_op_loop<op, ScalarType, VectorType>);
687}
688
689template <ArithmeticOperation op>
690void elementwise_arithm_op_quantized(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window)
691{
692 elementwise_op_quantized(in1, in2, out, window, &elementwise_arithm_op_quantized_scalar<op>,
693 &elementwise_arithm_op_quantized_broadcast_loop<op>,
694 &elementwise_arithm_op_quantized_loop<op>);
695}
696
697template <ComparisonOperation op>
698void elementwise_comp_op_quantized(const ITensor *in1, const ITensor *in2, ITensor *out, const Window &window)
699{
700 elementwise_op_quantized(in1, in2, out, window, &elementwise_comp_op_quantized_scalar<op>,
701 &elementwise_comp_op_quantized_broadcast_loop<op>,
702 &elementwise_comp_op_quantized_loop<op>);
703}
704
705std::function<void(const ITensor *, const ITensor *, ITensor *, const Window &)>
706configure_func(const ITensor *input1, const ITensor *input2, ITensor *output,
707 std::map<std::string, NEElementwiseOperationKernel::ElementwiseFunction *> map_function)
708{
709 std::string function_to_call("op_");
710 function_to_call += string_from_data_type(input1->info()->data_type()) + "_";
711 function_to_call += string_from_data_type(input2->info()->data_type()) + "_";
712 function_to_call += string_from_data_type(output->info()->data_type());
713
714 auto it = map_function.find(function_to_call);
715
716 if(it != map_function.end())
717 {
718 auto func = it->second;
719 return [func](const ITensor * input1, const ITensor * input2, ITensor * output, const Window & window)
720 {
721 func(input1, input2, output, window);
722 };
723 }
724 return nullptr;
725}
726
727template <ArithmeticOperation op>
728std::function<void(const ITensor *, const ITensor *, ITensor *, const Window &)>
729configure_arithm_func(const ITensor *input1, const ITensor *input2, ITensor *output)
730{
731 static std::map<std::string, NEElementwiseOperationKernel::ElementwiseFunction *> map_function =
732 {
733 { "op_F32_F32_F32", &elementwise_arithm_op<op, float, float32x4_t> },
734 { "op_S16_S16_S16", &elementwise_arithm_op<op, int16_t, int16x8_t> },
735 { "op_S32_S32_S32", &elementwise_arithm_op<op, int32_t, int32x4_t> },
736 { "op_QASYMM8_QASYMM8_QASYMM8", &elementwise_arithm_op_quantized<op> }
737 };
738#ifdef __ARM_FEATURE_FP16_VECTOR_ARITHMETIC
739 map_function["op_F16_F16_F16"] = &elementwise_arithm_op<op, float16_t, float16x8_t>;
740#endif /* __ARM_FEATURE_FP16_VECTOR_ARITHMETIC */
741
742 return configure_func(input1, input2, output, map_function);
743}
744
745template <ComparisonOperation op>
746std::function<void(const ITensor *input1, const ITensor *input2, ITensor *output, const Window &window)>
747configure_comp_func(const ITensor *input1, const ITensor *input2, ITensor *output)
748{
749 static std::map<std::string, NEElementwiseOperationKernel::ElementwiseFunction *> map_function =
750 {
751 { "op_F32_F32_U8", &elementwise_comp_op_32<op, float, float32x4_t> },
752 { "op_S16_S16_U8", &elementwise_comp_op_16<op, int16_t, int16x8_t> },
753 { "op_S32_S32_U8", &elementwise_comp_op_32<op, int32_t, int32x4_t> },
754 { "op_QASYMM8_QASYMM8_U8", &elementwise_comp_op_quantized<op> }
755 };
756#ifdef __ARM_FEATURE_FP16_VECTOR_ARITHMETIC
757 map_function["op_F16_F16_U8"] = &elementwise_comp_op_16<op, float16_t, float16x8_t>;
758#endif /* __ARM_FEATURE_FP16_VECTOR_ARITHMETIC */
759
760 return configure_func(input1, input2, output, map_function);
761}
762} // namespace
763
764NEElementwiseOperationKernel::NEElementwiseOperationKernel()
765 : _function(nullptr), _input1(nullptr), _input2(nullptr), _output(nullptr)
766{
767}
768
769Status NEElementwiseOperationKernel::validate_arguments_common(const ITensorInfo &input1, const ITensorInfo &input2, const ITensorInfo &output)
770{
giuros0192fd9432018-12-03 17:30:00 +0000771 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&input1, 1, DataType::QASYMM8, DataType::S16, DataType::F16, DataType::S32, DataType::F32);
772 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&input2, 1, DataType::QASYMM8, DataType::S16, DataType::F16, DataType::S32, DataType::F32);
George Wortd88590f2018-12-12 17:39:58 +0000773 ARM_COMPUTE_RETURN_ERROR_ON_CPU_F16_UNSUPPORTED(&input1);
giuros0192fd9432018-12-03 17:30:00 +0000774 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(&input1, &input2);
775
776 const TensorShape out_shape = TensorShape::broadcast_shape(input1.tensor_shape(), input2.tensor_shape());
777
778 ARM_COMPUTE_RETURN_ERROR_ON_MSG(out_shape.total_size() == 0, "Inputs are not broadcast compatible");
779
780 // Validate in case of configured output
781 if(output.total_size() > 0)
782 {
giuros0192fd9432018-12-03 17:30:00 +0000783 ARM_COMPUTE_RETURN_ERROR_ON_MSG(detail::have_different_dimensions(out_shape, output.tensor_shape(), 0),
784 "Wrong shape for output");
785 }
786
787 return Status{};
788}
giuros0192fd9432018-12-03 17:30:00 +0000789
giuros0192fd9432018-12-03 17:30:00 +0000790void NEElementwiseOperationKernel::configure_common(const ITensor *input1, const ITensor *input2, ITensor *output)
791{
792 ARM_COMPUTE_ERROR_ON_NULLPTR(input1, input2, output);
giuros0192fd9432018-12-03 17:30:00 +0000793
794 // Configure kernel window
795 const std::pair<TensorShape, ValidRegion> broadcast_pair = ITensorInfo::broadcast_shape_and_valid_region(*input1->info(), *input2->info());
796 const TensorShape &out_shape = broadcast_pair.first;
797 const ValidRegion &valid_region = broadcast_pair.second;
798
799 // Auto initialize output if not initialized
800 auto_init_if_empty(*output->info(), out_shape, 1, input1->info()->data_type());
801
802 Window win = calculate_max_window(valid_region);
803
giuros0192fd9432018-12-03 17:30:00 +0000804 _input1 = input1;
805 _input2 = input2;
806 _output = output;
807
giuros0192fd9432018-12-03 17:30:00 +0000808 INEKernel::configure(win);
809}
810
811void NEElementwiseOperationKernel::run(const Window &window, const ThreadInfo &info)
812{
George Wortd88590f2018-12-12 17:39:58 +0000813 ARM_COMPUTE_UNUSED(info, window);
giuros0192fd9432018-12-03 17:30:00 +0000814 ARM_COMPUTE_ERROR_ON_UNCONFIGURED_KERNEL(this);
815 ARM_COMPUTE_ERROR_ON_INVALID_SUBWINDOW(INEKernel::window(), window);
George Wortd88590f2018-12-12 17:39:58 +0000816 ARM_COMPUTE_ERROR_ON(_function == nullptr);
817 _function(_input1, _input2, _output, window);
giuros0192fd9432018-12-03 17:30:00 +0000818}
819
820/** Arithmetic operators (min, max, squared_diff) */
821
822void NEArithmeticOperationKernel::configure(ArithmeticOperation op, const ITensor *input1, const ITensor *input2, ITensor *output)
823{
George Wortd88590f2018-12-12 17:39:58 +0000824 ARM_COMPUTE_ERROR_THROW_ON(validate_arguments(*input1->info(), *input2->info(), *output->info()));
825 configure_common(input1, input2, output);
giuros0192fd9432018-12-03 17:30:00 +0000826 switch(op)
827 {
828 case ArithmeticOperation::MAX:
George Wortd88590f2018-12-12 17:39:58 +0000829 _function = configure_arithm_func<ArithmeticOperation::MAX>(input1, input2, output);
giuros0192fd9432018-12-03 17:30:00 +0000830 break;
831 case ArithmeticOperation::MIN:
George Wortd88590f2018-12-12 17:39:58 +0000832 _function = configure_arithm_func<ArithmeticOperation::MIN>(input1, input2, output);
giuros0192fd9432018-12-03 17:30:00 +0000833 break;
834 case ArithmeticOperation::SQUARED_DIFF:
George Wortd88590f2018-12-12 17:39:58 +0000835 _function = configure_arithm_func<ArithmeticOperation::SQUARED_DIFF>(input1, input2, output);
giuros0192fd9432018-12-03 17:30:00 +0000836 break;
837 default:
838 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
839 }
840}
841
George Wortd88590f2018-12-12 17:39:58 +0000842Status NEArithmeticOperationKernel::validate_arguments(const ITensorInfo &input1, const ITensorInfo &input2, const ITensorInfo &output)
843{
844 // Validate in case of configured output
845 if(output.total_size() > 0)
846 {
847 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(&input1, &output);
848 }
849 return validate_arguments_common(input1, input2, output);
850}
851
giuros0192fd9432018-12-03 17:30:00 +0000852Status NEArithmeticOperationKernel::validate(ArithmeticOperation op, const ITensorInfo *input1, const ITensorInfo *input2, const ITensorInfo *output)
853{
854 ARM_COMPUTE_UNUSED(op);
855 ARM_COMPUTE_RETURN_ERROR_ON_NULLPTR(input1, input2, output);
George Wortd88590f2018-12-12 17:39:58 +0000856 ARM_COMPUTE_RETURN_ON_ERROR(validate_arguments(*input1, *input2, *output));
giuros0192fd9432018-12-03 17:30:00 +0000857 return Status{};
858}
859
George Worta1e7e282019-01-15 11:00:29 +0000860/** The division operator */
861
862void NEDivisionOperationKernel::configure(const ITensor *input1, const ITensor *input2, ITensor *output)
863{
864 ARM_COMPUTE_ERROR_THROW_ON(validate_arguments(*input1->info(), *input2->info(), *output->info()));
865 configure_common(input1, input2, output);
866 _function = configure_arithm_func<ArithmeticOperation::DIV>(input1, input2, output);
867}
868
869Status NEDivisionOperationKernel::validate_arguments(const ITensorInfo &input1, const ITensorInfo &input2, const ITensorInfo &output)
870{
871 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&input1, 1, DataType::F16, DataType::F32);
872 return NEArithmeticOperationKernel::validate_arguments(input1, input2, output);
873}
874
875Status NEDivisionOperationKernel::validate(const ITensorInfo *input1, const ITensorInfo *input2, const ITensorInfo *output)
876{
877 ARM_COMPUTE_RETURN_ERROR_ON_NULLPTR(input1, input2, output);
878 ARM_COMPUTE_RETURN_ON_ERROR(validate_arguments(*input1, *input2, *output));
879 return Status{};
880}
881
George Wortd88590f2018-12-12 17:39:58 +0000882/** Comparison operators (equal, not equal, less than, greater than, less than or equal, greater than or equal) */
883
884void NEComparisonOperationKernel::configure(ComparisonOperation op, const ITensor *input1, const ITensor *input2, ITensor *output)
giuros0192fd9432018-12-03 17:30:00 +0000885{
George Wortd88590f2018-12-12 17:39:58 +0000886 ARM_COMPUTE_ERROR_THROW_ON(validate_arguments(*input1->info(), *input2->info(), *output->info()));
887 configure_common(input1, input2, output);
888 switch(op)
889 {
890 case ComparisonOperation::Equal:
891 _function = configure_comp_func<ComparisonOperation::Equal>(input1, input2, output);
892 break;
893 case ComparisonOperation::NotEqual:
894 _function = configure_comp_func<ComparisonOperation::NotEqual>(input1, input2, output);
895 break;
896 case ComparisonOperation::Greater:
897 _function = configure_comp_func<ComparisonOperation::Greater>(input1, input2, output);
898 break;
899 case ComparisonOperation::GreaterEqual:
900 _function = configure_comp_func<ComparisonOperation::GreaterEqual>(input1, input2, output);
901 break;
902 case ComparisonOperation::Less:
903 _function = configure_comp_func<ComparisonOperation::Less>(input1, input2, output);
904 break;
905 case ComparisonOperation::LessEqual:
906 _function = configure_comp_func<ComparisonOperation::LessEqual>(input1, input2, output);
907 break;
908 default:
909 ARM_COMPUTE_ERROR("NOT_SUPPORTED!");
910 }
911}
912
913Status NEComparisonOperationKernel::validate_arguments(const ITensorInfo &input1, const ITensorInfo &input2, const ITensorInfo &output)
914{
915 // Validate in case of configured output
916 if(output.total_size() > 0)
917 {
918 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&output, 1, DataType::U8);
919 }
920 return validate_arguments_common(input1, input2, output);
921}
922
923Status NEComparisonOperationKernel::validate(ComparisonOperation op, const ITensorInfo *input1, const ITensorInfo *input2, const ITensorInfo *output)
924{
925 ARM_COMPUTE_UNUSED(op);
926 ARM_COMPUTE_RETURN_ERROR_ON_NULLPTR(input1, input2, output);
927 ARM_COMPUTE_RETURN_ON_ERROR(validate_arguments(*input1, *input2, *output));
928 return Status{};
giuros0192fd9432018-12-03 17:30:00 +0000929}
930} // namespace arm_compute