blob: 31726444884f20bcc2da7c049b1a29831d0ffd0d [file] [log] [blame]
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +01001/*
Milos Puzovic13b623e2022-07-27 17:53:21 +00002 * Copyright (c) 2021-2022 Arm Limited.
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +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 */
Georgios Pinitas7891a732021-08-20 21:39:25 +010024#include "src/cpu/operators/CpuFullyConnected.h"
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +010025
26#include "arm_compute/core/Helpers.h"
27#include "arm_compute/core/ITensorPack.h"
28#include "arm_compute/core/Validate.h"
29#include "arm_compute/core/utils/misc/ShapeCalculator.h"
30#include "arm_compute/core/utils/quantization/AsymmHelpers.h"
31#include "arm_compute/runtime/NEON/NEScheduler.h"
ramelg013ae3d882021-09-12 23:07:47 +010032#include "src/common/utils/Log.h"
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +010033#include "src/core/helpers/AutoConfiguration.h"
34#include "src/core/helpers/MemoryHelpers.h"
Georgios Pinitas7891a732021-08-20 21:39:25 +010035#include "src/cpu/kernels/CpuTransposeKernel.h"
36#include "src/cpu/operators/CpuConvertFullyConnectedWeights.h"
37#include "src/cpu/operators/CpuFlatten.h"
38#include "src/cpu/operators/CpuGemm.h"
39#include "src/cpu/operators/CpuGemmLowpMatrixMultiplyCore.h"
40#include "src/cpu/utils/CpuAuxTensorHandler.h"
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +010041
42namespace arm_compute
43{
44namespace cpu
45{
46using namespace arm_compute::experimental;
47using namespace arm_compute::misc::shape_calculator;
48
49namespace
50{
51// Get min, max bound of a quantized asymmetric dst tensor, with the effect of fused activation
52std::pair<PixelValue, PixelValue> get_quantized_asymmetric_output_min_max(const QuantizationInfo &q_info, const ActivationLayerInfo &act_info, DataType data_type)
53{
54 PixelValue type_min{};
55 PixelValue type_max{};
Milos Puzovic13b623e2022-07-27 17:53:21 +000056 std::tie(type_min, type_max) = get_min_max(data_type);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +010057 const UniformQuantizationInfo q_unif = q_info.uniform();
58
59 if(act_info.enabled())
60 {
61 switch(act_info.activation())
62 {
63 case ActivationLayerInfo::ActivationFunction::RELU:
64 type_min = PixelValue(q_unif.offset);
65 break;
66 case ActivationLayerInfo::ActivationFunction::BOUNDED_RELU:
67 type_min = PixelValue(q_unif.offset);
68 type_max = PixelValue(act_info.a(), data_type, q_info);
69 break;
70 case ActivationLayerInfo::ActivationFunction::LU_BOUNDED_RELU:
71 type_min = PixelValue(act_info.b(), data_type, q_info);
72 type_max = PixelValue(act_info.a(), data_type, q_info);
73 break;
74 default:
75 ARM_COMPUTE_ERROR("Activation function not supported.");
76 break;
77 }
78 }
79
80 return std::make_pair(type_min, type_max);
81}
82
83Status get_gemmlowp_output_stage_info(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *dst, const ActivationLayerInfo &act,
84 GEMMLowpOutputStageInfo &gemmlowp_output_stage_info)
85{
86 const auto data_type = src->data_type();
87 const QuantizationInfo oq_info = dst->quantization_info();
88 const UniformQuantizationInfo iq_unif = src->quantization_info().uniform();
89 const UniformQuantizationInfo wq_unif = weights->quantization_info().uniform();
90 const UniformQuantizationInfo oq_unif = oq_info.uniform();
91
92 float multiplier = (iq_unif.scale * wq_unif.scale) / oq_unif.scale;
93 int32_t output_multiplier;
94 int32_t output_shift;
95
96 ARM_COMPUTE_RETURN_ON_ERROR(quantization::calculate_quantized_multiplier(multiplier, &output_multiplier, &output_shift));
97
98 PixelValue type_min{};
99 PixelValue type_max{};
100 std::tie(type_min, type_max) = get_quantized_asymmetric_output_min_max(oq_info, act, data_type);
101
102 gemmlowp_output_stage_info.gemmlowp_multiplier = output_multiplier;
103 gemmlowp_output_stage_info.gemmlowp_shift = output_shift;
104 gemmlowp_output_stage_info.gemmlowp_offset = oq_unif.offset;
105 gemmlowp_output_stage_info.type = GEMMLowpOutputStageType::QUANTIZE_DOWN_FIXEDPOINT;
106 gemmlowp_output_stage_info.gemmlowp_min_bound = type_min.get<int32_t>();
107 gemmlowp_output_stage_info.gemmlowp_max_bound = type_max.get<int32_t>();
108
109 return Status{};
110}
111
cfRodf2c022e2021-11-05 11:29:53 +0000112Status validate_mm(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, const ITensorInfo *dst, const ActivationLayerInfo &act, bool enable_fast_math)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100113{
114 if(is_data_type_quantized_asymmetric(src->data_type()))
115 {
116 // Since we need negative offsets for computing convolution, we need to change QuantizationInfo()
117 // Extract and negate src and weights offset
118 const QuantizationInfo src_quantization_info(src->quantization_info().uniform().scale, -src->quantization_info().uniform().offset);
119 const QuantizationInfo weights_quantization_info(weights->quantization_info().uniform().scale, -weights->quantization_info().uniform().offset);
120
121 GEMMLowpOutputStageInfo gemmlowp_output_stage_info;
122 ARM_COMPUTE_RETURN_ON_ERROR(get_gemmlowp_output_stage_info(src, weights, dst, act, gemmlowp_output_stage_info));
123
124 GEMMInfo gemm_info;
125 gemm_info.set_gemmlowp_output_stage(gemmlowp_output_stage_info);
cfRodf2c022e2021-11-05 11:29:53 +0000126 gemm_info.set_fast_math(enable_fast_math);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100127
128 // Validate gemmlowp function
129 TensorInfo src_info = src->clone()->set_quantization_info(src_quantization_info);
130 TensorInfo weights_info = weights->clone()->set_quantization_info(weights_quantization_info);
131 ARM_COMPUTE_RETURN_ON_ERROR(CpuGemmLowpMatrixMultiplyCore::validate(&src_info,
132 &weights_info,
133 biases,
134 dst,
135 gemm_info));
136 }
137 else
138 {
cfRodf2c022e2021-11-05 11:29:53 +0000139 GEMMInfo gemm_info(false, false, true /* Reshape weights only for the first run */);
140 gemm_info.set_fast_math(enable_fast_math);
141 ARM_COMPUTE_RETURN_ON_ERROR(CpuGemm::validate(src, weights, biases, dst, 1.f, 1.0f, gemm_info));
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100142 }
143
144 return Status{};
145}
146} // namespace
147
148CpuFullyConnected::CpuFullyConnected()
149 : _flatten(nullptr),
150 _convert_weights(nullptr),
151 _transpose_weights(nullptr),
152 _mm_gemm(nullptr),
153 _mm_gemmlowp(nullptr),
154 _flattened_src(),
155 _converted_weights(),
156 _reshaped_weights(),
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100157 _trans_weights(),
158 _trans_weights_idx(AuxTensorIdx::Count),
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100159 _aux_mem(Count),
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100160 _needs_weights_conversion(false),
161 _needs_weights_reshape(false),
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100162 _is_fc_after_conv(false),
163 _is_quantized_asymmetric(false),
cfRodf2c022e2021-11-05 11:29:53 +0000164 _is_prepared(false),
Milos Puzovic13b623e2022-07-27 17:53:21 +0000165 _enable_fast_math(false),
166 _fixed_format(false),
167 _weight_format(arm_compute::WeightFormat::UNSPECIFIED)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100168{
169}
170
171CpuFullyConnected::~CpuFullyConnected() = default;
172
173void CpuFullyConnected::configure_mm(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, ITensorInfo *dst, const ActivationLayerInfo &act)
174{
175 if(_is_quantized_asymmetric)
176 {
177 // Since we need negative offsets for computing convolution, we need to change QuantizationInfo()
178 // Extract and negate src and weights offset
179 const QuantizationInfo src_quantization_info(src->quantization_info().uniform().scale, -src->quantization_info().uniform().offset);
180 const QuantizationInfo weights_quantization_info(weights->quantization_info().uniform().scale, -weights->quantization_info().uniform().offset);
181
182 TensorInfo src_info = src->clone()->set_quantization_info(src_quantization_info);
183 TensorInfo weights_info = weights->clone()->set_quantization_info(weights_quantization_info);
184
185 // Configure gemmlowp function and output stage for asymmetric quantized types
186 GEMMLowpOutputStageInfo gemmlowp_output_stage_info;
187 const Status status = get_gemmlowp_output_stage_info(&src_info, &weights_info, dst, act, gemmlowp_output_stage_info);
188 ARM_COMPUTE_ERROR_ON(status.error_code() != ErrorCode::OK);
189
190 GEMMInfo gemm_info;
191 gemm_info.set_gemmlowp_output_stage(gemmlowp_output_stage_info);
192 gemm_info.set_activation_info(act);
cfRodf2c022e2021-11-05 11:29:53 +0000193 gemm_info.set_fast_math(_enable_fast_math);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100194 _mm_gemmlowp = std::make_unique<CpuGemmLowpMatrixMultiplyCore>();
195 _mm_gemmlowp->configure(&src_info, &weights_info, biases, dst, gemm_info);
196 }
197 else
198 {
199 // Configure matrix multiply kernel
200 GEMMInfo gemm_info(false, false, true /* Reshape weights only for the first run */);
201 gemm_info.set_activation_info(act);
cfRodf2c022e2021-11-05 11:29:53 +0000202 gemm_info.set_fast_math(_enable_fast_math);
Milos Puzovic13b623e2022-07-27 17:53:21 +0000203 gemm_info.set_fixed_format(_fixed_format);
204 gemm_info.set_weight_format(_weight_format);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100205 _mm_gemm = std::make_unique<CpuGemm>();
206 _mm_gemm->configure(src, weights, biases, dst, 1.f, 1.0f, gemm_info);
207 }
208}
209
210void CpuFullyConnected::configure_conv_fc(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, ITensorInfo *dst, const ActivationLayerInfo &act)
211{
212 ARM_COMPUTE_ERROR_ON((weights->dimension(1) != (src->dimension(0) * src->dimension(1) * src->dimension(2))));
213
214 // If the fully connected layer is called after a convolution layer, the src tensor must be linearized
215
216 // Initialize output tensor for flatten
217 auto_init_if_empty(_flattened_src, src->clone()->set_tensor_shape(compute_flatten_shape(src)));
218
219 _flatten = std::make_unique<CpuFlatten>();
220 _flatten->configure(src, &_flattened_src);
221
222 // Configure matrix multiply kernel
223 configure_mm(&_flattened_src, weights, biases, dst, act);
224}
225
226void CpuFullyConnected::configure_fc_fc(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, ITensorInfo *dst, const ActivationLayerInfo &act)
227{
228 ARM_COMPUTE_ERROR_ON(src->dimension(0) != weights->dimension(1));
229
230 // Configure matrix multiply kernel
231 configure_mm(src, weights, biases, dst, act);
232}
233
234void CpuFullyConnected::configure(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, ITensorInfo *dst,
Milos Puzovic13b623e2022-07-27 17:53:21 +0000235 FullyConnectedLayerInfo fc_info, const WeightsInfo &weights_info)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100236{
237 // Perform validate step
238 ARM_COMPUTE_ERROR_ON_NULLPTR(src, weights, dst);
239 ARM_COMPUTE_ERROR_THROW_ON(CpuFullyConnected::validate(src,
240 weights,
241 biases != nullptr ? biases : nullptr,
242 dst,
243 fc_info));
ramelg013ae3d882021-09-12 23:07:47 +0100244 ARM_COMPUTE_LOG_PARAMS(src, weights, biases, dst, fc_info);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100245
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100246 _needs_weights_conversion = false;
247 _needs_weights_reshape = fc_info.transpose_weights ? !fc_info.are_weights_reshaped : false;
248 _needs_weights_reshape = _needs_weights_reshape && !fc_info.retain_internal_weights;
249 _is_fc_after_conv = true;
250 _is_quantized_asymmetric = is_data_type_quantized_asymmetric(src->data_type());
251 _is_prepared = false;
252 _trans_weights_idx = AuxTensorIdx::Count;
cfRodf2c022e2021-11-05 11:29:53 +0000253 _enable_fast_math = fc_info.enable_fast_math;
Milos Puzovic13b623e2022-07-27 17:53:21 +0000254 _fixed_format = weights_info.weight_format() != WeightFormat::UNSPECIFIED;
255 _weight_format = weights_info.weight_format();
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100256
257 // With the Fully Connected layer we can have 4 different cases:
258 // 1) Convolution layer -> Fully Connected layer without batches
259 // 2) Fully Connected layer -> Fully Connected layer without batches
260 // 3) Convolution layer -> Fully Connected layer with batches
261 // 4) Fully Connected layer -> Fully Connected layer with batches
262
263 const ITensorInfo *weights_to_use = weights;
264
265 // Check if we have a fully connected layer with batches
266 const bool is_batched_fc_layer = dst->dimension(1) > 1;
267 if(is_batched_fc_layer)
268 {
Milos Puzovic13b623e2022-07-27 17:53:21 +0000269 _is_fc_after_conv = (TensorShape::num_max_dimensions >= 4) && (std::equal(src->tensor_shape().cbegin() + 3, src->tensor_shape().cend(), dst->tensor_shape().cbegin() + 1));
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100270 }
271 else
272 {
273 _is_fc_after_conv = src->num_dimensions() > 1;
274 }
275
276 // Reshape weights if needed
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100277 if(_needs_weights_reshape)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100278 {
279 // Reshape the weights
280 _transpose_weights = std::make_unique<kernels::CpuTransposeKernel>();
281 _transpose_weights->configure(weights, &_reshaped_weights);
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100282 weights_to_use = &_reshaped_weights;
283 _trans_weights_idx = AuxTensorIdx::TransposedWeights;
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100284 }
285
286 // Convert weights if needed
287 if(_is_fc_after_conv && (src->data_layout() != fc_info.weights_trained_layout))
288 {
289 // Convert weights
290 _convert_weights = std::make_unique<CpuConvertFullyConnectedWeights>();
291 _convert_weights->configure(weights_to_use,
292 &_converted_weights,
293 src->tensor_shape(),
294 fc_info.weights_trained_layout);
295
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100296 weights_to_use = &_converted_weights;
297 _needs_weights_conversion = true;
298 _trans_weights_idx = AuxTensorIdx::ConvertedWeights;
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100299 }
300
301 if(_is_fc_after_conv)
302 {
303 // Fully Connected layer after a Convolution Layer without batches
304 configure_conv_fc(src, weights_to_use, biases, dst, fc_info.activation_info);
305 }
306 else
307 {
308 // Fully Connected layer after a Fully Connected Layer without batches
309 configure_fc_fc(src, weights_to_use, biases, dst, fc_info.activation_info);
310 }
311
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100312 // Retain the tensorinfo with the weights to use
313 if(_needs_weights_reshape || _needs_weights_conversion)
314 {
315 _trans_weights = *weights_to_use;
316 }
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100317
318 // Set auxiliary memory requirements
319 auto gemm_mem_req = (_is_quantized_asymmetric) ? _mm_gemmlowp->workspace() : _mm_gemm->workspace();
320 for(unsigned int i = 0; i < gemm_mem_req.size(); ++i)
321 {
322 _aux_mem[i] = gemm_mem_req[i];
323 }
324
325 if(_aux_mem[Pretranspose].size > 0)
326 {
Giorgio Arena63e0beb2021-09-24 14:04:27 +0100327 // Release permuted weights at the end of prepare as they are further transposed by the assembly dispatch
328 // Do not release them if biases are dynamic and data type is quantized, since the weights tensor will be used for biases offset calculation
Milos Puzovic13b623e2022-07-27 17:53:21 +0000329 _aux_mem[TransposedWeights] = MemoryInfo(offset_int_vec(TransposedWeights), (_is_quantized_asymmetric && biases
330 && !(biases->are_values_constant())) ? MemoryLifetime::Persistent : MemoryLifetime::Prepare,
Giorgio Arena63e0beb2021-09-24 14:04:27 +0100331 _reshaped_weights.total_size());
Milos Puzovic13b623e2022-07-27 17:53:21 +0000332 _aux_mem[ConvertedWeights] = MemoryInfo(offset_int_vec(ConvertedWeights), MemoryLifetime::Prepare, _converted_weights.total_size());
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100333 }
334 else
335 {
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100336 _aux_mem[TransposedWeights] = MemoryInfo(offset_int_vec(TransposedWeights), _needs_weights_conversion ? MemoryLifetime::Prepare : MemoryLifetime::Persistent, _reshaped_weights.total_size());
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100337 _aux_mem[ConvertedWeights] = MemoryInfo(offset_int_vec(ConvertedWeights), MemoryLifetime::Persistent, _converted_weights.total_size());
338 }
339 _aux_mem[FlattenedSrc] = MemoryInfo(offset_int_vec(FlattenedSrc), MemoryLifetime::Temporary, _flattened_src.total_size());
340}
341
Milos Puzovic13b623e2022-07-27 17:53:21 +0000342Status CpuFullyConnected::has_opt_impl(arm_compute::WeightFormat &expected_weight_format, const ITensorInfo *src, const ITensorInfo *weights,
343 const ITensorInfo *biases, const ITensorInfo *dst, FullyConnectedLayerInfo fc_info, WeightsInfo weights_info)
344{
345 GEMMInfo gemm_info(false, false, true /* Reshape weights only for the first run */);
346 gemm_info.set_activation_info(fc_info.activation_info);
347 gemm_info.set_fast_math(fc_info.enable_fast_math);
348 gemm_info.set_fixed_format(weights_info.weight_format() != WeightFormat::UNSPECIFIED);
349 gemm_info.set_weight_format(weights_info.weight_format());
350
351 return CpuGemm::has_opt_impl(expected_weight_format, src, weights, biases, dst, gemm_info);
352}
353
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100354Status CpuFullyConnected::validate(const ITensorInfo *src, const ITensorInfo *weights, const ITensorInfo *biases, const ITensorInfo *dst,
355 FullyConnectedLayerInfo fc_info)
356{
357 ARM_COMPUTE_UNUSED(fc_info.retain_internal_weights);
358 ARM_COMPUTE_RETURN_ERROR_ON_NULLPTR(src, weights, dst);
359 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(src, 1, DataType::QASYMM8, DataType::QASYMM8_SIGNED, DataType::F16, DataType::F32);
360 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(src, weights, dst);
361 ARM_COMPUTE_RETURN_ERROR_ON(weights->num_dimensions() > 2);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100362 ARM_COMPUTE_RETURN_ERROR_ON(fc_info.activation_info.enabled() && is_data_type_quantized(src->data_type()) && fc_info.activation_info.activation() != ActivationLayerInfo::ActivationFunction::RELU
363 && fc_info.activation_info.activation() != ActivationLayerInfo::ActivationFunction::BOUNDED_RELU && fc_info.activation_info.activation() != ActivationLayerInfo::ActivationFunction::LU_BOUNDED_RELU);
Giorgio Arena63e0beb2021-09-24 14:04:27 +0100364 ARM_COMPUTE_RETURN_ERROR_ON(!weights->are_values_constant() && (!fc_info.are_weights_reshaped || fc_info.transpose_weights));
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100365
366 bool weights_reshaped = fc_info.transpose_weights ? fc_info.are_weights_reshaped : true;
367 bool is_fc_after_conv = true;
368
369 const ITensorInfo &flatten_src = TensorInfo(src->clone()->set_is_resizable(true).reset_padding().set_tensor_shape(compute_flatten_shape(src)));
370 const ITensorInfo &reshaped_weights = TensorInfo(weights->clone()->set_is_resizable(true).reset_padding().set_tensor_shape(compute_transposed_shape(*weights)));
371 const ITensorInfo &converted_weights = weights_reshaped ? TensorInfo(weights->clone()->set_is_resizable(true).reset_padding()) : TensorInfo(*reshaped_weights.clone());
372
373 // With the Fully Connected layer we can have 4 different cases:
374 // 1) Convolution layer -> Fully Connected layer without batches
375 // 2) Fully Connected layer -> Fully Connected layer without batches
376 // 3) Convolution layer -> Fully Connected layer with batches
377 // 4) Fully Connected layer -> Fully Connected layer with batches
378
379 const ITensorInfo *src_to_use = src;
380 const ITensorInfo *weights_to_use = weights;
381
382 // Check if we have a fully connected layer with batches
383 const bool is_batched_fc_layer = dst->dimension(1) > 1;
384
Giorgio Arena63e0beb2021-09-24 14:04:27 +0100385 if(biases != nullptr)
386 {
387 ARM_COMPUTE_RETURN_ERROR_ON(biases->num_dimensions() > 1);
388 if(is_data_type_quantized(src->data_type()))
389 {
390 ARM_COMPUTE_RETURN_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(biases, 1, DataType::S32);
391 }
392 else
393 {
394 ARM_COMPUTE_RETURN_ERROR_ON_MISMATCHING_DATA_TYPES(src, biases);
395 }
396 }
397
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100398 if(is_batched_fc_layer)
399 {
Milos Puzovic13b623e2022-07-27 17:53:21 +0000400 is_fc_after_conv = (TensorShape::num_max_dimensions >= 4) && (std::equal(src->tensor_shape().cbegin() + 3, src->tensor_shape().cend(), dst->tensor_shape().cbegin() + 1));
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100401 }
402 else
403 {
404 is_fc_after_conv = src->num_dimensions() > 1;
405 }
406
407 if(!weights_reshaped)
408 {
409 // Validate reshape weights kernel
410 ARM_COMPUTE_RETURN_ON_ERROR(kernels::CpuTransposeKernel::validate(weights, &reshaped_weights));
411 weights_to_use = &reshaped_weights;
412 }
413
414 if(is_fc_after_conv && (src->data_layout() != fc_info.weights_trained_layout))
415 {
416 // Validate convert weights kernel
417 ARM_COMPUTE_RETURN_ON_ERROR(CpuConvertFullyConnectedWeights::validate(weights_to_use,
418 &converted_weights,
419 src->tensor_shape(),
420 fc_info.weights_trained_layout));
421 weights_to_use = &converted_weights;
422 }
423
424 if(is_fc_after_conv)
425 {
426 // Fully Connected layer after a Convolution Layer without batches
427 ARM_COMPUTE_RETURN_ERROR_ON((weights_to_use->dimension(1) != (src->dimension(0) * src->dimension(1) * src->dimension(2))));
428
429 // Validate flatten kernel
430 ARM_COMPUTE_RETURN_ON_ERROR(CpuFlatten::validate(src, &flatten_src));
431 src_to_use = &flatten_src;
432 }
433 else
434 {
435 // Fully Connected layer after a Fully Connected Layer without batches
436 ARM_COMPUTE_RETURN_ERROR_ON(src->dimension(0) != weights_to_use->dimension(1));
437 }
438 // Validate matrix multiply kernel
cfRodf2c022e2021-11-05 11:29:53 +0000439 ARM_COMPUTE_RETURN_ON_ERROR(validate_mm(src_to_use, weights_to_use, biases, dst, fc_info.activation_info, fc_info.enable_fast_math));
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100440
441 return Status{};
442}
443
444void CpuFullyConnected::run(ITensorPack &tensors)
445{
446 prepare(tensors);
447
448 auto src = tensors.get_const_tensor(ACL_SRC_0);
449
450 CpuAuxTensorHandler flattened_src(offset_int_vec(FlattenedSrc), _flattened_src, tensors, false);
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100451 CpuAuxTensorHandler transformed_wei(offset_int_vec(_trans_weights_idx), _trans_weights, tensors, false);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100452
453 // Linearize src if it comes from a convolutional layer
454 if(_is_fc_after_conv)
455 {
456 ITensorPack flatten_pack{ { ACL_SRC, src }, { ACL_DST, flattened_src.get() } };
457 _flatten->run(flatten_pack);
458 }
459
460 ITensorPack gemm_pack = tensors;
461 gemm_pack.add_const_tensor(ACL_SRC_0, (_is_fc_after_conv) ? flattened_src.get() : src);
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100462 if(_needs_weights_reshape || _needs_weights_conversion)
463 {
464 gemm_pack.add_const_tensor(ACL_SRC_1, transformed_wei.get());
465 }
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100466
467 // Run matrix multiply
468 if(_is_quantized_asymmetric)
469 {
470 _mm_gemmlowp->run(gemm_pack);
471 }
472 else
473 {
474 _mm_gemm->run(gemm_pack);
475 }
476}
477
478void CpuFullyConnected::prepare(ITensorPack &tensors)
479{
480 if(!_is_prepared)
481 {
482 auto weights = tensors.get_const_tensor(ACL_SRC_1);
483
484 CpuAuxTensorHandler reshaped_weights(offset_int_vec(TransposedWeights), _reshaped_weights, tensors, false);
485 CpuAuxTensorHandler converted_weights(offset_int_vec(ConvertedWeights), _converted_weights, tensors, false);
486
487 // Pointer to current weights
488 const ITensor *cur_weights = weights;
489
490 // Reshape of the weights (happens only once)
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100491 if(_needs_weights_reshape)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100492 {
493 // Run reshape weights kernel and mark weights as unused
494 ITensorPack transpose_pack{ { ACL_SRC, weights }, { ACL_DST, reshaped_weights.get() } };
495 NEScheduler::get().schedule_op(_transpose_weights.get(), Window::DimY, _transpose_weights->window(), transpose_pack);
496
497 cur_weights->mark_as_unused();
498 cur_weights = reshaped_weights.get();
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100499 }
500
501 // Convert weights if needed (happens only once)
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100502 if(_needs_weights_conversion)
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100503 {
504 ITensorPack convert_pack{ { ACL_SRC, cur_weights }, { ACL_DST, converted_weights.get() } };
505 _convert_weights->run(convert_pack);
506
507 cur_weights->mark_as_unused();
508 cur_weights = converted_weights.get();
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100509 }
510
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100511 ITensorPack gemm_pack = tensors;
512 gemm_pack.add_const_tensor(ACL_SRC_1, cur_weights);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100513
514 // Prepare GEMM prepare and release unused weights
515 if(!_is_quantized_asymmetric)
516 {
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100517 _mm_gemm->prepare(gemm_pack);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100518 }
519 else
520 {
Georgios Pinitasfa1db172021-08-12 06:28:09 +0100521 _mm_gemmlowp->prepare(gemm_pack);
Michele Di Giorgiod9cdf142021-07-02 15:17:08 +0100522 }
523
524 _is_prepared = true;
525 }
526}
527
528experimental::MemoryRequirements CpuFullyConnected::workspace() const
529{
530 return _aux_mem;
531}
532} // namespace cpu
533} // namespace arm_compute