blob: d1dfcd9bb1dd130691600d831899ea20558f932a [file] [log] [blame]
steniu0127b386c2017-07-18 17:37:43 +01001/*
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +00002 * Copyright (c) 2017-2020 ARM Limited.
steniu0127b386c2017-07-18 17:37:43 +01003 *
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/CL/kernels/CLDirectConvolutionLayerKernel.h"
25
26#include "arm_compute/core/AccessWindowStatic.h"
27#include "arm_compute/core/CL/CLHelpers.h"
28#include "arm_compute/core/CL/CLKernelLibrary.h"
Vidhya Sudhan Loganathanf1f49062018-05-25 13:21:26 +010029#include "arm_compute/core/CL/CLValidate.h"
steniu0127b386c2017-07-18 17:37:43 +010030#include "arm_compute/core/CL/ICLTensor.h"
31#include "arm_compute/core/Error.h"
32#include "arm_compute/core/Helpers.h"
33#include "arm_compute/core/IAccessWindow.h"
34#include "arm_compute/core/ITensor.h"
35#include "arm_compute/core/Types.h"
Gian Marco Iodice5cb4d6a2017-08-08 10:53:00 +010036#include "arm_compute/core/Utils.h"
Giorgio Arenac0f54432018-03-16 14:02:34 +000037#include "arm_compute/core/utils/misc/ShapeCalculator.h"
Chunosovd621bca2017-11-03 17:33:15 +070038#include "arm_compute/core/utils/quantization/AsymmHelpers.h"
Matthew Bentham758b5ba2020-03-05 23:37:48 +000039#include "support/StringSupport.h"
steniu0127b386c2017-07-18 17:37:43 +010040
Michele Di Giorgio14cbfb22019-10-23 10:53:10 +010041namespace arm_compute
42{
Georgios Pinitas30902ed2017-11-14 15:32:57 +000043namespace
44{
Georgios Pinitas631c41a2017-12-06 11:53:03 +000045Status validate_arguments(const ITensorInfo *input, const ITensorInfo *weights, const ITensorInfo *biases, const ITensorInfo *output, const PadStrideInfo &conv_info)
Georgios Pinitas30902ed2017-11-14 15:32:57 +000046{
Vidhya Sudhan Loganathanf1f49062018-05-25 13:21:26 +010047 ARM_COMPUTE_RETURN_ERROR_ON_F16_UNSUPPORTED(input);
Sheri Zhang681f2d42020-02-20 11:23:08 +000048 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(input, 1, DataType::QASYMM8_SIGNED, DataType::QASYMM8, DataType::F16, DataType::F32);
Georgios Pinitas30902ed2017-11-14 15:32:57 +000049 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(input, weights);
Pablo Tello3d319462018-06-21 15:13:17 +010050
51 const DataLayout data_layout = input->data_layout();
52 const int width_idx = get_data_layout_dimension_index(data_layout, DataLayoutDimension::WIDTH);
53 const int height_idx = get_data_layout_dimension_index(data_layout, DataLayoutDimension::HEIGHT);
54 const int channel_idx = get_data_layout_dimension_index(data_layout, DataLayoutDimension::CHANNEL);
55
56 ARM_COMPUTE_RETURN_ERROR_ON_MSG(weights->dimension(width_idx) != weights->dimension(height_idx), "Weights should have same width and height");
Michalis Spyrou45091732019-05-13 17:41:01 +010057 ARM_COMPUTE_RETURN_ERROR_ON_MSG(weights->dimension(width_idx) != 1 && weights->dimension(width_idx) != 3 && weights->dimension(width_idx) != 5 && weights->dimension(width_idx) != 9,
58 "Kernel sizes other than 1x1, 3x3, 5x5 or 9x9 are not supported");
Pablo Tello3d319462018-06-21 15:13:17 +010059 ARM_COMPUTE_RETURN_ERROR_ON_MSG(weights->dimension(channel_idx) != input->dimension(channel_idx),
Georgios Pinitas30902ed2017-11-14 15:32:57 +000060 "Weights feature map dimension should match the respective input's one");
Pablo Tello3d319462018-06-21 15:13:17 +010061 ARM_COMPUTE_RETURN_ERROR_ON_MSG(weights->num_dimensions() > 4, "Weights can be at most 4 dimensional");
62 ARM_COMPUTE_RETURN_ERROR_ON_MSG((weights->dimension(width_idx) == 1) && std::get<0>(conv_info.stride()) > 3, "Strides larger than 3 not supported for 1x1 convolution.");
63 ARM_COMPUTE_RETURN_ERROR_ON_MSG((weights->dimension(width_idx) == 3 || weights->dimension(width_idx) == 5) && std::get<0>(conv_info.stride()) > 2,
Georgios Pinitas30902ed2017-11-14 15:32:57 +000064 "Strides larger than 2 not supported for 3x3 convolution.");
Sang-Hoon Parkab5b1a22019-10-15 09:29:13 +010065
66 const auto data_type = input->data_type();
67
68 if(weights->dimension(width_idx) == 9)
69 {
70 const auto supported_data_layout = is_data_type_quantized(data_type) ? DataLayout::NCHW : DataLayout::NHWC;
71 const auto error_message = std::string("Only " + string_from_data_layout(supported_data_layout) + " layout is supported for 9x9 convolution with " + string_from_data_type(
Michele Di Giorgio14cbfb22019-10-23 10:53:10 +010072 data_type)
73 + " type");
Sang-Hoon Parkab5b1a22019-10-15 09:29:13 +010074
75 ARM_COMPUTE_RETURN_ERROR_ON_MSG((supported_data_layout != data_layout), error_message.c_str());
76 }
Georgios Pinitas30902ed2017-11-14 15:32:57 +000077
78 if(biases != nullptr)
79 {
80 if(is_data_type_quantized_asymmetric(input->data_type()))
81 {
82 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(biases, 1, DataType::S32);
83 }
84 else
85 {
86 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(weights, biases);
87 }
88 ARM_COMPUTE_RETURN_ERROR_ON_MSG(biases->dimension(0) != weights->dimension(3),
89 "Biases size and number of input feature maps should match");
90 ARM_COMPUTE_RETURN_ERROR_ON_MSG(biases->num_dimensions() > 1,
91 "Biases should be one dimensional");
92 }
93
94 // Checks performed when output is configured
95 if(output->total_size() != 0)
96 {
97 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DIMENSIONS(output->tensor_shape(),
Giorgio Arenac0f54432018-03-16 14:02:34 +000098 misc::shape_calculator::compute_deep_convolution_shape(*input, *weights, conv_info));
Georgios Pinitas30902ed2017-11-14 15:32:57 +000099 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(input, output);
Georgios Pinitas30902ed2017-11-14 15:32:57 +0000100 }
101
Michele Di Giorgio14cbfb22019-10-23 10:53:10 +0100102 if(is_data_type_quantized(data_type))
103 {
104 const UniformQuantizationInfo iqinfo = input->quantization_info().uniform();
105 const UniformQuantizationInfo wqinfo = weights->quantization_info().uniform();
106 const UniformQuantizationInfo oqinfo = output->quantization_info().uniform();
107
108 float multiplier = iqinfo.scale * wqinfo.scale / oqinfo.scale;
109 int output_multiplier = 0;
110 int output_shift = 0;
111 ARM_COMPUTE_RETURN_ON_ERROR(quantization::calculate_quantized_multiplier(multiplier, &output_multiplier, &output_shift));
112 }
Georgios Pinitas631c41a2017-12-06 11:53:03 +0000113 return Status{};
Georgios Pinitas30902ed2017-11-14 15:32:57 +0000114}
115
Pablo Tello3d319462018-06-21 15:13:17 +0100116inline bool can_run_optimized_kernel_for_bifrost(GPUTarget gpu_target, unsigned int conv_stride_x, unsigned int conv_stride_y, unsigned int kernel_size,
117 DataType data_type, DataLayout data_layout)
Giorgio Arena59486342017-12-01 10:42:47 +0000118{
Georgios Pinitasa34286e2018-09-04 12:18:50 +0100119 return gpu_target_is_in(gpu_target,
120 GPUTarget::G71, GPUTarget::G72, GPUTarget::G76,
121 GPUTarget::G51, GPUTarget::G51BIG, GPUTarget::G51LIT,
122 GPUTarget::G52, GPUTarget::G52LIT)
123 && (kernel_size <= 5)
124 && (conv_stride_x == 1) && (conv_stride_y == 1)
125 && (data_type == DataType::F32)
126 && (data_layout == DataLayout::NCHW);
Pablo Tello3d319462018-06-21 15:13:17 +0100127}
Giorgio Arena59486342017-12-01 10:42:47 +0000128
Michalis Spyrou45091732019-05-13 17:41:01 +0100129inline bool can_run_optimized_kernel_for_bifrost_nhwc(GPUTarget gpu_target, unsigned int conv_stride_x, unsigned int conv_stride_y, unsigned int kernel_size,
130 DataType data_type, DataLayout data_layout)
131{
132 return gpu_target_is_in(gpu_target,
133 GPUTarget::G71, GPUTarget::G72, GPUTarget::G76,
134 GPUTarget::G51, GPUTarget::G51BIG, GPUTarget::G51LIT,
135 GPUTarget::G52, GPUTarget::G52LIT)
136 && (kernel_size == 9)
137 && (conv_stride_x == 1) && (conv_stride_y == 1)
138 && (data_type == DataType::F32)
139 && (data_layout == DataLayout::NHWC);
140}
141
Pablo Tello3d319462018-06-21 15:13:17 +0100142inline void setup_num_elems(unsigned int &num_elems_read_per_iteration_x, unsigned int &num_elems_read_per_iteration_y,
143 unsigned int &num_elems_written_per_iteration_x, unsigned int &num_elems_written_per_iteration_y,
144 unsigned int kernel_size, const PadStrideInfo &conv_info, const GPUTarget target, ITensorInfo *input)
145{
146 const DataType data_type = input->data_type();
147 const DataLayout data_layout = input->data_layout();
148 unsigned int conv_stride_x = std::get<0>(conv_info.stride());
149 unsigned int conv_stride_y = std::get<1>(conv_info.stride());
Giorgio Arena59486342017-12-01 10:42:47 +0000150
Pablo Tello3d319462018-06-21 15:13:17 +0100151 const bool run_optimized_bifrost = can_run_optimized_kernel_for_bifrost(target, conv_stride_x, conv_stride_y, kernel_size, data_type, data_layout);
Giorgio Arena59486342017-12-01 10:42:47 +0000152
Pablo Tello3d319462018-06-21 15:13:17 +0100153 if(run_optimized_bifrost)
Giorgio Arena59486342017-12-01 10:42:47 +0000154 {
155 // Configure kernel window
Giorgio Arena59486342017-12-01 10:42:47 +0000156 switch(kernel_size)
157 {
158 case 1:
159 {
160 num_elems_read_per_iteration_x = 4;
161 num_elems_read_per_iteration_y = 4;
162 num_elems_written_per_iteration_x = 4;
163 num_elems_written_per_iteration_y = 4;
164 break;
165 }
166 case 3:
167 {
168 num_elems_read_per_iteration_x = 6;
169 num_elems_read_per_iteration_y = 5;
170 num_elems_written_per_iteration_x = 4;
171 num_elems_written_per_iteration_y = 3;
172 break;
173 }
174 case 5:
175 {
176 num_elems_read_per_iteration_x = 8;
177 num_elems_read_per_iteration_y = 6;
178 num_elems_written_per_iteration_x = 4;
179 num_elems_written_per_iteration_y = 2;
180 break;
181 }
182 default:
183 {
184 ARM_COMPUTE_ERROR("Kernel size not optimized for Bifrost");
185 }
186 }
187 }
Michalis Spyrou45091732019-05-13 17:41:01 +0100188 else if(data_layout == DataLayout::NCHW)
Giorgio Arena59486342017-12-01 10:42:47 +0000189 {
Giorgio Arena59486342017-12-01 10:42:47 +0000190 num_elems_read_per_iteration_y = kernel_size;
191 num_elems_written_per_iteration_x = 8;
192 num_elems_written_per_iteration_y = 1;
Anthony Barbiercc9fed52017-12-13 10:46:00 +0000193 switch(kernel_size)
194 {
195 case 1:
196 switch(conv_stride_x)
197 {
198 case 1:
199 num_elems_read_per_iteration_x = 8;
200 break;
201 case 2:
202 num_elems_read_per_iteration_x = 16;
203 break;
204 case 3:
205 switch(input->element_size())
206 {
207 case 1:
208 num_elems_read_per_iteration_x = 28;
209 break;
210 case 2:
211 num_elems_read_per_iteration_x = 24;
212 break;
213 case 4:
214 num_elems_read_per_iteration_x = 22;
215 break;
216 default:
217 ARM_COMPUTE_ERROR("Invalid data size");
218 }
219 break;
220 default:
221 ARM_COMPUTE_ERROR("Invalid convolution stride X");
222 }
223 break;
224 case 3:
225 switch(conv_stride_x)
226 {
227 case 1:
228 num_elems_read_per_iteration_x = 10;
229 break;
230 case 2:
231 num_elems_read_per_iteration_x = 17;
232 break;
233 default:
234 ARM_COMPUTE_ERROR("Invalid convolution stride X");
235 }
236 break;
237 case 5:
238 switch(conv_stride_x)
239 {
240 case 1:
241 num_elems_read_per_iteration_x = 12;
242 break;
243 case 2:
244 num_elems_read_per_iteration_x = 20;
245 break;
246 default:
247 ARM_COMPUTE_ERROR("Invalid convolution stride X");
248 }
249 break;
Sang-Hoon Parkab5b1a22019-10-15 09:29:13 +0100250 case 9:
251 switch(conv_stride_x)
252 {
253 case 1:
254 num_elems_read_per_iteration_x = 16;
255 break;
256 case 2:
257 num_elems_read_per_iteration_x = 24;
258 break;
259 default:
260 ARM_COMPUTE_ERROR("Invalid convolution stride X");
261 }
262 break;
Anthony Barbiercc9fed52017-12-13 10:46:00 +0000263 default:
264 ARM_COMPUTE_ERROR("Invalid direct convolution size");
265 }
Giorgio Arena59486342017-12-01 10:42:47 +0000266 }
Michalis Spyrou45091732019-05-13 17:41:01 +0100267 else // data_layout == NHWC
Pablo Tello3d319462018-06-21 15:13:17 +0100268 {
Michalis Spyrou45091732019-05-13 17:41:01 +0100269 const bool run_optimized_bifrost_nhwc = can_run_optimized_kernel_for_bifrost_nhwc(target, conv_stride_x, conv_stride_y, kernel_size, data_type, data_layout);
270
Pablo Tello3d319462018-06-21 15:13:17 +0100271 num_elems_written_per_iteration_x = 1;
Michalis Spyrou45091732019-05-13 17:41:01 +0100272
273 if(run_optimized_bifrost_nhwc)
274 {
275 num_elems_read_per_iteration_x = 4;
276 }
giuros01c878f1f2019-07-09 11:01:34 +0100277 else
278 {
279 num_elems_read_per_iteration_x = 1;
280 }
Michalis Spyrou45091732019-05-13 17:41:01 +0100281
Pablo Tello3d319462018-06-21 15:13:17 +0100282 switch(kernel_size)
283 {
284 case 1:
285 switch(conv_stride_x)
286 {
287 case 1:
288 num_elems_read_per_iteration_y = 8;
289 num_elems_written_per_iteration_y = 8;
290 break;
291 case 2:
292 num_elems_read_per_iteration_y = 16;
293 num_elems_written_per_iteration_y = 8;
294 break;
295 default:
296 ARM_COMPUTE_ERROR("Invalid convolution stride X");
297 }
298 break;
299 case 3:
300 switch(conv_stride_x)
301 {
302 case 1:
303 num_elems_read_per_iteration_y = 10;
304 num_elems_written_per_iteration_y = 8;
305 break;
306 case 2:
307 num_elems_read_per_iteration_y = 17;
308 num_elems_written_per_iteration_y = 8;
309 break;
310 default:
311 ARM_COMPUTE_ERROR("Invalid convolution stride X");
312 }
313 break;
314 case 5:
315 switch(conv_stride_x)
316 {
317 case 1:
318 num_elems_read_per_iteration_y = 12;
319 num_elems_written_per_iteration_y = 8;
320 break;
321 case 2:
322 num_elems_read_per_iteration_y = 20;
323 num_elems_written_per_iteration_y = 8;
324 break;
325 default:
326 ARM_COMPUTE_ERROR("Invalid convolution stride X");
327 }
328 break;
Michalis Spyrou45091732019-05-13 17:41:01 +0100329 case 9:
330 switch(conv_stride_x)
331 {
332 case 1:
333 num_elems_read_per_iteration_y = 16;
334 num_elems_written_per_iteration_y = 8;
335 break;
336 case 2:
337 num_elems_read_per_iteration_y = 24;
338 num_elems_written_per_iteration_y = 8;
339 break;
340 default:
341 ARM_COMPUTE_ERROR("Invalid convolution stride X");
342 }
343 break;
Pablo Tello3d319462018-06-21 15:13:17 +0100344 default:
345 ARM_COMPUTE_ERROR("Not implemented.");
346 break;
347 }
348 }
349}
350
351std::pair<Status, Window> validate_and_configure_window(ITensorInfo *input, ITensorInfo *weights, ITensorInfo *output, const PadStrideInfo &conv_info, const GPUTarget target)
352{
353 const DataLayout data_layout = input->data_layout();
354 const int width_idx = get_data_layout_dimension_index(data_layout, DataLayoutDimension::WIDTH);
355 const unsigned int kernel_size = weights->dimension(width_idx);
356
357 // Get convolved dimensions
358 TensorShape output_shape = misc::shape_calculator::compute_deep_convolution_shape(*input, *weights, conv_info);
359
360 // Output auto inizialitation if not yet initialized
Georgios Pinitasf52cd782019-03-25 14:06:14 +0000361 // TODO(COMPMID-2078): input->clone()->set_tensor_shape(output_shape) doesn't work with subtensors for grouped direct convolutions (AlexNet).
Pablo Tello3d319462018-06-21 15:13:17 +0100362 auto_init_if_empty(*output, output_shape,
363 1,
364 input->data_type(),
365 input->quantization_info());
366
367 unsigned int num_elems_read_per_iteration_x = 0;
368 unsigned int num_elems_read_per_iteration_y = 0;
369 unsigned int num_elems_written_per_iteration_x = 0;
370 unsigned int num_elems_written_per_iteration_y = 0;
371
372 unsigned int conv_pad_left = conv_info.pad_left();
373 unsigned int conv_pad_top = conv_info.pad_top();
374 unsigned int conv_stride_x = std::get<0>(conv_info.stride());
375 unsigned int conv_stride_y = std::get<1>(conv_info.stride());
376
377 setup_num_elems(num_elems_read_per_iteration_x, num_elems_read_per_iteration_y,
378 num_elems_written_per_iteration_x, num_elems_written_per_iteration_y,
379 kernel_size, conv_info, target, input);
380
Giorgio Arena59486342017-12-01 10:42:47 +0000381 // Create window and update padding
Anthony Barbiercc9fed52017-12-13 10:46:00 +0000382 bool window_changed = false;
383 Window win = calculate_max_window(*output, Steps(num_elems_written_per_iteration_x, num_elems_written_per_iteration_y));
Giorgio Arena59486342017-12-01 10:42:47 +0000384
Pablo Tello3d319462018-06-21 15:13:17 +0100385 if(data_layout == DataLayout::NHWC)
386 {
387 AccessWindowStatic input_access(input, 0, -conv_pad_left,
giuros01c878f1f2019-07-09 11:01:34 +0100388 ceil_to_multiple(input->dimension(0), num_elems_read_per_iteration_x),
Pablo Tello3d319462018-06-21 15:13:17 +0100389 ceil_to_multiple(input->dimension(1) + conv_info.pad_right(), num_elems_read_per_iteration_y));
390 AccessWindowStatic weights_access(weights, 0, 0, weights->dimension(0), weights->dimension(1));
391 AccessWindowRectangle output_access(output, 0, 0, num_elems_written_per_iteration_x, num_elems_written_per_iteration_y);
392 window_changed = update_window_and_padding(win, input_access, weights_access, output_access);
393 output_access.set_valid_region(win, ValidRegion(Coordinates(), output->tensor_shape()));
394 Status err = (window_changed) ? ARM_COMPUTE_CREATE_ERROR(ErrorCode::RUNTIME_ERROR, "Insufficient Padding!") : Status{};
395 return std::make_pair(err, win);
396 }
397 else if(data_layout == DataLayout::NCHW)
398 {
399 AccessWindowRectangle input_access(input, -conv_pad_left, -conv_pad_top, num_elems_read_per_iteration_x, num_elems_read_per_iteration_y, conv_stride_x, conv_stride_y);
400 AccessWindowStatic weights_access(weights, 0, 0, kernel_size, kernel_size);
401 AccessWindowRectangle output_access(output, 0, 0, num_elems_written_per_iteration_x, num_elems_written_per_iteration_y);
402 window_changed = update_window_and_padding(win, input_access, weights_access, output_access);
403 output_access.set_valid_region(win, ValidRegion(Coordinates(), output->tensor_shape()));
404 Status err = (window_changed) ? ARM_COMPUTE_CREATE_ERROR(ErrorCode::RUNTIME_ERROR, "Insufficient Padding!") : Status{};
405 return std::make_pair(err, win);
406 }
407 else
408 {
409 ARM_COMPUTE_ERROR("Not supported");
410 }
Giorgio Arena59486342017-12-01 10:42:47 +0000411}
412} // namespace
413
414CLDirectConvolutionLayerKernel::CLDirectConvolutionLayerKernel()
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000415 : _input(nullptr), _biases(nullptr), _weights(nullptr), _output(nullptr), _data_layout(DataLayout::UNKNOWN), _border_size(0), _conv_stride_x(0), _conv_stride_y(0)
Giorgio Arena59486342017-12-01 10:42:47 +0000416{
417}
418
419BorderSize CLDirectConvolutionLayerKernel::border_size() const
420{
421 return _border_size;
422}
423
424void CLDirectConvolutionLayerKernel::configure(const ICLTensor *input, const ICLTensor *weights, const ICLTensor *biases, ICLTensor *output, const PadStrideInfo &conv_info)
425{
426 ARM_COMPUTE_ERROR_ON_NULLPTR(input, weights, output);
427
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000428 _data_layout = input->info()->data_layout();
429 const int width_idx = get_data_layout_dimension_index(_data_layout, DataLayoutDimension::WIDTH);
430 const int height_idx = get_data_layout_dimension_index(_data_layout, DataLayoutDimension::HEIGHT);
431 const int channel_idx = get_data_layout_dimension_index(_data_layout, DataLayoutDimension::CHANNEL);
Pablo Tello3d319462018-06-21 15:13:17 +0100432
433 const unsigned int kernel_size = weights->info()->dimension(width_idx);
Giorgio Arena59486342017-12-01 10:42:47 +0000434 const DataType data_type = input->info()->data_type();
435
436 // Get convolved dimensions
Giorgio Arenac0f54432018-03-16 14:02:34 +0000437 TensorShape output_shape = misc::shape_calculator::compute_deep_convolution_shape(*input->info(), *weights->info(), conv_info);
Giorgio Arena59486342017-12-01 10:42:47 +0000438
439 // Output auto inizialitation if not yet initialized
Georgios Pinitasf52cd782019-03-25 14:06:14 +0000440 // TODO(COMPMID-2078): input->clone()->set_tensor_shape(output_shape) doesn't work with subtensors for grouped direct convolutions (AlexNet).
Giorgio Arena59486342017-12-01 10:42:47 +0000441 auto_init_if_empty(*output->info(),
442 output_shape,
443 1,
444 input->info()->data_type(),
Giorgio Arena59486342017-12-01 10:42:47 +0000445 input->info()->quantization_info());
446
447 // Perform validation step
448 ARM_COMPUTE_ERROR_THROW_ON(validate_arguments(input->info(),
449 weights->info(),
450 (biases != nullptr) ? biases->info() : nullptr,
451 output->info(),
452 conv_info));
453
454 _conv_stride_x = std::get<0>(conv_info.stride());
455 _conv_stride_y = std::get<1>(conv_info.stride());
Pablo Tello3d319462018-06-21 15:13:17 +0100456
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000457 if(_data_layout == DataLayout::NHWC)
Pablo Tello3d319462018-06-21 15:13:17 +0100458 {
459 _border_size = BorderSize(conv_info.pad_left(), 0, conv_info.pad_right(), 0);
460 }
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000461 else if(_data_layout == DataLayout::NCHW)
Pablo Tello3d319462018-06-21 15:13:17 +0100462 {
463 _border_size = BorderSize(conv_info.pad_top(), conv_info.pad_right(), conv_info.pad_bottom(), conv_info.pad_left());
464 }
465 else
466 {
467 ARM_COMPUTE_ERROR("Not supported");
468 }
Giorgio Arena59486342017-12-01 10:42:47 +0000469
470 _input = input;
471 _weights = weights;
472 _output = output;
473 _biases = biases;
474
Michalis Spyroua9676112018-02-22 18:07:43 +0000475 const GPUTarget gpu_target = get_target();
Giorgio Arena59486342017-12-01 10:42:47 +0000476
477 std::stringstream kernel_name;
478 kernel_name << "direct_convolution" << kernel_size << "x" << kernel_size;
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000479 if(_data_layout == DataLayout::NHWC)
Pablo Tello3d319462018-06-21 15:13:17 +0100480 {
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000481 kernel_name << "_" << lower_string(string_from_data_layout(_data_layout));
Pablo Tello3d319462018-06-21 15:13:17 +0100482 }
Giorgio Arena59486342017-12-01 10:42:47 +0000483
484 CLBuildOptions build_options;
485 build_options.add_option_if(_biases != nullptr, std::string("-DHAS_BIAS"));
486
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000487 const bool run_optimized_for_bifrost = can_run_optimized_kernel_for_bifrost(gpu_target, _conv_stride_x, _conv_stride_y, kernel_size, data_type, _data_layout);
Pablo Tello3d319462018-06-21 15:13:17 +0100488
489 if(run_optimized_for_bifrost)
Giorgio Arena59486342017-12-01 10:42:47 +0000490 {
Pablo Tello3d319462018-06-21 15:13:17 +0100491 build_options.add_option(std::string("-DWEIGHTS_DEPTH=" + support::cpp11::to_string(_weights->info()->dimension(channel_idx))));
Giorgio Arena59486342017-12-01 10:42:47 +0000492
493 kernel_name << "_f32_bifrost";
494 _kernel = static_cast<cl::Kernel>(CLKernelLibrary::get().create_kernel(kernel_name.str(), build_options.options()));
Giorgio Arena59486342017-12-01 10:42:47 +0000495 }
496 else
497 {
Giorgio Arena59486342017-12-01 10:42:47 +0000498 build_options.add_option(std::string("-DDATA_TYPE=" + get_cl_type_from_data_type(data_type)));
499 build_options.add_option(std::string("-DDATA_SIZE=" + get_data_size_from_data_type(data_type)));
Pablo Tello3d319462018-06-21 15:13:17 +0100500 build_options.add_option(std::string("-DWEIGHTS_DEPTH=" + support::cpp11::to_string(_weights->info()->dimension(channel_idx))));
Giorgio Arena59486342017-12-01 10:42:47 +0000501 build_options.add_option(std::string("-DSTRIDE_X=" + support::cpp11::to_string(_conv_stride_x)));
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000502 if(_data_layout == DataLayout::NHWC)
Pablo Tello3d319462018-06-21 15:13:17 +0100503 {
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000504 const bool run_optimized_for_bifrost_nhwc = can_run_optimized_kernel_for_bifrost_nhwc(gpu_target, _conv_stride_x, _conv_stride_y, kernel_size, data_type, _data_layout);
Pablo Tello3d319462018-06-21 15:13:17 +0100505 build_options.add_option(std::string("-DDATA_LAYOUT_NHWC=1"));
506 build_options.add_option(std::string("-DDST_HEIGHT=" + support::cpp11::to_string(_output->info()->dimension(height_idx))));
507 build_options.add_option(std::string("-DDST_WIDTH=" + support::cpp11::to_string(_output->info()->dimension(width_idx))));
508 build_options.add_option(std::string("-DSRC_HEIGHT=" + support::cpp11::to_string(_input->info()->dimension(height_idx))));
509 build_options.add_option(std::string("-DSRC_WIDTH=" + support::cpp11::to_string(_input->info()->dimension(width_idx))));
510 build_options.add_option(std::string("-DPAD_LEFT=" + support::cpp11::to_string(conv_info.pad_left())));
511 build_options.add_option(std::string("-DPAD_TOP=" + support::cpp11::to_string(conv_info.pad_top())));
Giorgio Arena3c4bf0c2020-03-02 09:49:29 +0000512 build_options.add_option(std::string("-DPAD_BOTTOM=" + support::cpp11::to_string(conv_info.pad_bottom())));
Pablo Tello3d319462018-06-21 15:13:17 +0100513 build_options.add_option(std::string("-DSTRIDE_Y=" + support::cpp11::to_string(_conv_stride_y)));
Michalis Spyrou45091732019-05-13 17:41:01 +0100514 if(run_optimized_for_bifrost_nhwc)
515 {
516 const unsigned int num_elems_read_per_iteration_x = 4;
517 _border_size.right = num_elems_read_per_iteration_x;
518 build_options.add_option("-DVEC_SIZE=" + support::cpp11::to_string(num_elems_read_per_iteration_x));
519 }
Pablo Tello3d319462018-06-21 15:13:17 +0100520 }
Vidhya Sudhan Loganathan7485d5a2018-07-04 09:34:00 +0100521 build_options.add_option(std::string("-DDATA_TYPE_PROMOTED=" + get_cl_type_from_data_type(data_type)));
Michele Di Giorgio14cbfb22019-10-23 10:53:10 +0100522
523 if(is_data_type_quantized(data_type))
524 {
525 const UniformQuantizationInfo iqinfo = _input->info()->quantization_info().uniform();
526 const UniformQuantizationInfo wqinfo = _weights->info()->quantization_info().uniform();
527 const UniformQuantizationInfo oqinfo = _output->info()->quantization_info().uniform();
528
529 float multiplier = iqinfo.scale * wqinfo.scale / oqinfo.scale;
530 int output_multiplier = 0;
531 int output_shift = 0;
532 quantization::calculate_quantized_multiplier(multiplier, &output_multiplier, &output_shift);
533 build_options.add_option("-DOUTPUT_MULTIPLIER=" + support::cpp11::to_string(output_multiplier));
534 build_options.add_option("-DOUTPUT_SHIFT=" + support::cpp11::to_string(output_shift));
535 build_options.add_option("-DKERNEL_SIZE=" + support::cpp11::to_string(kernel_size));
536
537 // Create kernel
538 _kernel = static_cast<cl::Kernel>(CLKernelLibrary::get().create_kernel("direct_convolution_quantized", build_options.options()));
539
540 // Set static kernel arguments
541 unsigned int idx = 3 * num_arguments_per_3D_tensor() + ((_biases != nullptr) ? num_arguments_per_1D_tensor() : 0) + 1;
542 _kernel.setArg(idx++, -iqinfo.offset);
543 _kernel.setArg(idx++, -wqinfo.offset);
544 _kernel.setArg(idx++, oqinfo.offset);
545 }
546 else
547 {
548 // Create kernel
549 _kernel = static_cast<cl::Kernel>(CLKernelLibrary::get().create_kernel(kernel_name.str(), build_options.options()));
550 }
Giorgio Arena59486342017-12-01 10:42:47 +0000551 }
552
553 // Configure kernel window
554 auto win_config = validate_and_configure_window(input->info(), weights->info(), output->info(), conv_info, gpu_target);
555 ARM_COMPUTE_ERROR_THROW_ON(win_config.first);
Anthony Barbierb6eb3532018-08-08 13:20:04 +0100556 ICLKernel::configure_internal(win_config.second);
Giorgio Arena59486342017-12-01 10:42:47 +0000557
Giorgio Arena59486342017-12-01 10:42:47 +0000558 // Set config_id for enabling LWS tuning
559 _config_id = "direct_convolution_";
560 _config_id += lower_string(string_from_data_type(data_type));
561 _config_id += "_";
562 _config_id += support::cpp11::to_string(kernel_size);
563 _config_id += "_";
Georgios Pinitas15997872018-02-19 13:58:22 +0000564 _config_id += support::cpp11::to_string(border_size().left);
Giorgio Arena59486342017-12-01 10:42:47 +0000565 _config_id += "_";
Georgios Pinitas15997872018-02-19 13:58:22 +0000566 _config_id += support::cpp11::to_string(border_size().top);
Giorgio Arena59486342017-12-01 10:42:47 +0000567 _config_id += "_";
Georgios Pinitas15997872018-02-19 13:58:22 +0000568 _config_id += support::cpp11::to_string(border_size().right);
Giorgio Arena59486342017-12-01 10:42:47 +0000569 _config_id += "_";
Georgios Pinitas15997872018-02-19 13:58:22 +0000570 _config_id += support::cpp11::to_string(border_size().bottom);
Giorgio Arena59486342017-12-01 10:42:47 +0000571 _config_id += "_";
572 _config_id += support::cpp11::to_string(_conv_stride_x);
573 _config_id += "_";
574 _config_id += support::cpp11::to_string(_conv_stride_y);
575 _config_id += "_";
Pablo Tello3d319462018-06-21 15:13:17 +0100576 _config_id += support::cpp11::to_string(output->info()->dimension(width_idx));
Giorgio Arena59486342017-12-01 10:42:47 +0000577 _config_id += "_";
Pablo Tello3d319462018-06-21 15:13:17 +0100578 _config_id += support::cpp11::to_string(output->info()->dimension(height_idx));
579 _config_id += "_";
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000580 _config_id += lower_string(string_from_data_layout(_data_layout));
Giorgio Arena59486342017-12-01 10:42:47 +0000581}
582
Georgios Pinitas631c41a2017-12-06 11:53:03 +0000583Status CLDirectConvolutionLayerKernel::validate(const ITensorInfo *input, const ITensorInfo *weights, const ITensorInfo *biases, const ITensorInfo *output, const PadStrideInfo &conv_info,
584 const GPUTarget target)
Giorgio Arena59486342017-12-01 10:42:47 +0000585{
586 ARM_COMPUTE_RETURN_ON_ERROR(validate_arguments(input, weights, biases, output, conv_info));
587 ARM_COMPUTE_RETURN_ON_ERROR(validate_and_configure_window(input->clone().get(), weights->clone().get(), output->clone().get(), conv_info, target).first);
588
Georgios Pinitas631c41a2017-12-06 11:53:03 +0000589 return Status{};
Giorgio Arena59486342017-12-01 10:42:47 +0000590}
591
SiCong Lic51b72f2017-07-28 14:46:20 +0100592void CLDirectConvolutionLayerKernel::run(const Window &window, cl::CommandQueue &queue)
steniu0127b386c2017-07-18 17:37:43 +0100593{
594 ARM_COMPUTE_ERROR_ON_UNCONFIGURED_KERNEL(this);
595 ARM_COMPUTE_ERROR_ON_INVALID_SUBWINDOW(IKernel::window(), window);
596
597 // Get initial windows
598 Window slice = window.first_slice_window_3D();
599 Window win_in = window;
600
Jaroslaw Rzepecki2ecbada2017-11-29 13:51:34 +0000601 win_in.adjust(Window::DimX, -_border_size.left, true);
602 win_in.adjust(Window::DimY, -_border_size.top, true);
steniu0127b386c2017-07-18 17:37:43 +0100603
Georgios Pinitas7fdcfb12020-01-09 16:45:46 +0000604 const int width_idx = get_data_layout_dimension_index(_data_layout, DataLayoutDimension::WIDTH);
605 const int height_idx = get_data_layout_dimension_index(_data_layout, DataLayoutDimension::HEIGHT);
steniu0127b386c2017-07-18 17:37:43 +0100606
Pablo Tello3d319462018-06-21 15:13:17 +0100607 win_in.set_dimension_step(width_idx, window[width_idx].step() * _conv_stride_x);
608 win_in.set_dimension_step(height_idx, window[height_idx].step() * _conv_stride_y);
609
610 Window slice_in = win_in.first_slice_window_3D();
611 unsigned int idx1 = 2 * num_arguments_per_3D_tensor();
steniu0127b386c2017-07-18 17:37:43 +0100612 add_3D_tensor_argument(idx1, _weights, slice);
613
614 if(_biases != nullptr)
615 {
616 Window slice_biases;
SiCong Li86b53332017-08-23 11:02:43 +0100617 slice_biases.use_tensor_dimensions(_biases->info()->tensor_shape());
steniu0127b386c2017-07-18 17:37:43 +0100618 add_1D_tensor_argument(idx1, _biases, slice_biases);
619 }
620
Gian Marco Iodice5cb4d6a2017-08-08 10:53:00 +0100621 _kernel.setArg(idx1++, static_cast<unsigned int>(_weights->info()->strides_in_bytes()[3]));
622
steniu0127b386c2017-07-18 17:37:43 +0100623 do
624 {
625 unsigned int idx = 0;
626 add_3D_tensor_argument(idx, _input, slice_in);
627 add_3D_tensor_argument(idx, _output, slice);
Anthony Barbierb6eb3532018-08-08 13:20:04 +0100628 enqueue(queue, *this, slice, lws_hint());
steniu0127b386c2017-07-18 17:37:43 +0100629 }
630 while(window.slide_window_slice_3D(slice) && win_in.slide_window_slice_3D(slice_in));
631}
Michele Di Giorgio14cbfb22019-10-23 10:53:10 +0100632} // namespace arm_compute