blob: 20e1369e730d7ec48775e8e2f8ad3da2259532f8 [file] [log] [blame]
Anthony Barbier2a07e182017-08-04 18:20:27 +01001/*
Giorgio Arena33b103b2021-01-08 10:37:15 +00002 * Copyright (c) 2017-2021 Arm Limited.
Anthony Barbier2a07e182017-08-04 18:20:27 +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
25#include "utils/GraphUtils.h"
hakanardof36ac352018-02-16 10:06:34 +010026
Georgios Pinitascac13b12018-04-27 19:07:19 +010027#include "arm_compute/core/Helpers.h"
28#include "arm_compute/core/Types.h"
Georgios Pinitas12be7ab2018-07-03 12:06:23 +010029#include "arm_compute/graph/Logger.h"
hakanardof36ac352018-02-16 10:06:34 +010030#include "arm_compute/runtime/SubTensor.h"
Michalis Spyrou6bff1952019-10-02 17:22:11 +010031
32#pragma GCC diagnostic push
33#pragma GCC diagnostic ignored "-Wunused-parameter"
Georgios Pinitas7c3b9242018-06-21 19:01:25 +010034#include "utils/ImageLoader.h"
Michalis Spyrou6bff1952019-10-02 17:22:11 +010035#pragma GCC diagnostic pop
Anthony Barbier2a07e182017-08-04 18:20:27 +010036#include "utils/Utils.h"
37
Michalis Spyrou7c60c992019-10-10 14:33:47 +010038#include <inttypes.h>
Gian Marco44ec2e72017-10-19 14:13:38 +010039#include <iomanip>
Georgios Pinitas12be7ab2018-07-03 12:06:23 +010040#include <limits>
Anthony Barbier2a07e182017-08-04 18:20:27 +010041
42using namespace arm_compute::graph_utils;
43
Georgios Pinitascac13b12018-04-27 19:07:19 +010044namespace
45{
Anthony Barbier4ead11a2018-08-06 09:25:36 +010046std::pair<arm_compute::TensorShape, arm_compute::PermutationVector> compute_permutation_parameters(const arm_compute::TensorShape &shape,
Georgios Pinitascac13b12018-04-27 19:07:19 +010047 arm_compute::DataLayout data_layout)
48{
49 // Set permutation parameters if needed
50 arm_compute::TensorShape permuted_shape = shape;
51 arm_compute::PermutationVector perm;
52 // Permute only if num_dimensions greater than 2
53 if(shape.num_dimensions() > 2)
54 {
55 perm = (data_layout == arm_compute::DataLayout::NHWC) ? arm_compute::PermutationVector(2U, 0U, 1U) : arm_compute::PermutationVector(1U, 2U, 0U);
56
57 arm_compute::PermutationVector perm_shape = (data_layout == arm_compute::DataLayout::NCHW) ? arm_compute::PermutationVector(2U, 0U, 1U) : arm_compute::PermutationVector(1U, 2U, 0U);
58 arm_compute::permute(permuted_shape, perm_shape);
59 }
60
61 return std::make_pair(permuted_shape, perm);
62}
63} // namespace
64
Georgios Pinitasbe2772a2018-08-17 15:33:39 +010065TFPreproccessor::TFPreproccessor(float min_range, float max_range)
66 : _min_range(min_range), _max_range(max_range)
67{
68}
Georgios Pinitas140fdc72018-02-16 11:42:38 +000069void TFPreproccessor::preprocess(ITensor &tensor)
70{
giuros01351bd132019-08-23 14:27:30 +010071 if(tensor.info()->data_type() == DataType::F32)
72 {
73 preprocess_typed<float>(tensor);
74 }
75 else if(tensor.info()->data_type() == DataType::F16)
76 {
77 preprocess_typed<half>(tensor);
78 }
79 else
80 {
81 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
82 }
83}
84
85template <typename T>
86void TFPreproccessor::preprocess_typed(ITensor &tensor)
87{
Georgios Pinitas140fdc72018-02-16 11:42:38 +000088 Window window;
89 window.use_tensor_dimensions(tensor.info()->tensor_shape());
90
Georgios Pinitasbe2772a2018-08-17 15:33:39 +010091 const float range = _max_range - _min_range;
Georgios Pinitas140fdc72018-02-16 11:42:38 +000092 execute_window_loop(window, [&](const Coordinates & id)
93 {
giuros01351bd132019-08-23 14:27:30 +010094 const T value = *reinterpret_cast<T *>(tensor.ptr_to_element(id));
95 float res = value / 255.f; // Normalize to [0, 1]
96 res = res * range + _min_range; // Map to [min_range, max_range]
97 *reinterpret_cast<T *>(tensor.ptr_to_element(id)) = res;
Georgios Pinitas140fdc72018-02-16 11:42:38 +000098 });
99}
100
Georgios Pinitasb54c6442019-04-03 13:18:14 +0100101CaffePreproccessor::CaffePreproccessor(std::array<float, 3> mean, bool bgr, float scale)
102 : _mean(mean), _bgr(bgr), _scale(scale)
Georgios Pinitas140fdc72018-02-16 11:42:38 +0000103{
104 if(_bgr)
105 {
106 std::swap(_mean[0], _mean[2]);
107 }
108}
109
110void CaffePreproccessor::preprocess(ITensor &tensor)
111{
giuros01351bd132019-08-23 14:27:30 +0100112 if(tensor.info()->data_type() == DataType::F32)
113 {
114 preprocess_typed<float>(tensor);
115 }
116 else if(tensor.info()->data_type() == DataType::F16)
117 {
118 preprocess_typed<half>(tensor);
119 }
120 else
121 {
122 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
123 }
124}
125
126template <typename T>
127void CaffePreproccessor::preprocess_typed(ITensor &tensor)
128{
Georgios Pinitas140fdc72018-02-16 11:42:38 +0000129 Window window;
130 window.use_tensor_dimensions(tensor.info()->tensor_shape());
Georgios Pinitas7d66a8e2018-07-17 12:28:42 +0100131 const int channel_idx = get_data_layout_dimension_index(tensor.info()->data_layout(), DataLayoutDimension::CHANNEL);
132
Georgios Pinitas140fdc72018-02-16 11:42:38 +0000133 execute_window_loop(window, [&](const Coordinates & id)
134 {
giuros01351bd132019-08-23 14:27:30 +0100135 const T value = *reinterpret_cast<T *>(tensor.ptr_to_element(id)) - T(_mean[id[channel_idx]]);
136 *reinterpret_cast<T *>(tensor.ptr_to_element(id)) = value * T(_scale);
Georgios Pinitas140fdc72018-02-16 11:42:38 +0000137 });
138}
139
Anthony Barbier2a07e182017-08-04 18:20:27 +0100140PPMWriter::PPMWriter(std::string name, unsigned int maximum)
141 : _name(std::move(name)), _iterator(0), _maximum(maximum)
142{
143}
144
145bool PPMWriter::access_tensor(ITensor &tensor)
146{
147 std::stringstream ss;
148 ss << _name << _iterator << ".ppm";
Gian Marco44ec2e72017-10-19 14:13:38 +0100149
150 arm_compute::utils::save_to_ppm(tensor, ss.str());
Anthony Barbier2a07e182017-08-04 18:20:27 +0100151
152 _iterator++;
153 if(_maximum == 0)
154 {
155 return true;
156 }
157 return _iterator < _maximum;
158}
159
160DummyAccessor::DummyAccessor(unsigned int maximum)
161 : _iterator(0), _maximum(maximum)
162{
163}
164
165bool DummyAccessor::access_tensor(ITensor &tensor)
166{
167 ARM_COMPUTE_UNUSED(tensor);
Anthony Barbier8a042112018-08-21 18:16:53 +0100168 bool ret = _maximum == 0 || _iterator < _maximum;
Anthony Barbier2a07e182017-08-04 18:20:27 +0100169 if(_iterator == _maximum)
170 {
171 _iterator = 0;
172 }
173 else
174 {
175 _iterator++;
176 }
177 return ret;
178}
179
Isabella Gottardia7acb3c2019-01-08 13:48:44 +0000180NumPyAccessor::NumPyAccessor(std::string npy_path, TensorShape shape, DataType data_type, DataLayout data_layout, std::ostream &output_stream)
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100181 : _npy_tensor(), _filename(std::move(npy_path)), _output_stream(output_stream)
182{
Isabella Gottardia7acb3c2019-01-08 13:48:44 +0000183 NumPyBinLoader loader(_filename, data_layout);
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100184
185 TensorInfo info(shape, 1, data_type);
Isabella Gottardia7acb3c2019-01-08 13:48:44 +0000186 info.set_data_layout(data_layout);
187
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100188 _npy_tensor.allocator()->init(info);
189 _npy_tensor.allocator()->allocate();
190
191 loader.access_tensor(_npy_tensor);
192}
193
194template <typename T>
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000195void NumPyAccessor::access_numpy_tensor(ITensor &tensor, T tolerance)
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100196{
Gian Marco Iodicead486e22018-08-07 17:17:06 +0100197 const int num_elements = tensor.info()->tensor_shape().total_size();
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000198 int num_mismatches = utils::compare_tensor<T>(tensor, _npy_tensor, tolerance);
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100199 float percentage_mismatches = static_cast<float>(num_mismatches) / num_elements;
200
201 _output_stream << "Results: " << 100.f - (percentage_mismatches * 100) << " % matches with the provided output[" << _filename << "]." << std::endl;
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000202 _output_stream << " " << num_elements - num_mismatches << " out of " << num_elements << " matches with the provided output[" << _filename << "]." << std::endl
203 << std::endl;
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100204}
205
206bool NumPyAccessor::access_tensor(ITensor &tensor)
207{
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000208 ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&tensor, 1, DataType::F32, DataType::QASYMM8);
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100209 ARM_COMPUTE_ERROR_ON(_npy_tensor.info()->dimension(0) != tensor.info()->dimension(0));
210
211 switch(tensor.info()->data_type())
212 {
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000213 case DataType::QASYMM8:
214 access_numpy_tensor<qasymm8_t>(tensor, 0);
215 break;
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100216 case DataType::F32:
Isabella Gottardi0ae5de92019-03-14 10:32:11 +0000217 access_numpy_tensor<float>(tensor, 0.0001f);
Isabella Gottardi88d5b222018-04-06 12:24:55 +0100218 break;
219 default:
220 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
221 }
222
223 return false;
224}
225
Isabella Gottardicd4e9ab2019-11-05 17:50:27 +0000226#ifdef ARM_COMPUTE_ASSERTS_ENABLED
227PrintAccessor::PrintAccessor(std::ostream &output_stream, IOFormatInfo io_fmt)
228 : _output_stream(output_stream), _io_fmt(io_fmt)
229{
230}
231
232bool PrintAccessor::access_tensor(ITensor &tensor)
233{
234 tensor.print(_output_stream, _io_fmt);
235 return false;
236}
237#endif /* ARM_COMPUTE_ASSERTS_ENABLED */
238
Isabella Gottardi2ea37612019-07-16 11:48:51 +0100239SaveNumPyAccessor::SaveNumPyAccessor(std::string npy_name, const bool is_fortran)
240 : _npy_name(std::move(npy_name)), _is_fortran(is_fortran)
241{
242}
243
244bool SaveNumPyAccessor::access_tensor(ITensor &tensor)
245{
246 ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&tensor, 1, DataType::F32);
247
248 utils::save_to_npy(tensor, _npy_name, _is_fortran);
249
250 return false;
251}
252
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100253ImageAccessor::ImageAccessor(std::string filename, bool bgr, std::unique_ptr<IPreprocessor> preprocessor)
Anthony Barbier8a042112018-08-21 18:16:53 +0100254 : _already_loaded(false), _filename(std::move(filename)), _bgr(bgr), _preprocessor(std::move(preprocessor))
Gian Marco44ec2e72017-10-19 14:13:38 +0100255{
256}
257
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100258bool ImageAccessor::access_tensor(ITensor &tensor)
Gian Marco44ec2e72017-10-19 14:13:38 +0100259{
Anthony Barbier8a042112018-08-21 18:16:53 +0100260 if(!_already_loaded)
Georgios Pinitascac13b12018-04-27 19:07:19 +0100261 {
Anthony Barbier8a042112018-08-21 18:16:53 +0100262 auto image_loader = utils::ImageLoaderFactory::create(_filename);
263 ARM_COMPUTE_EXIT_ON_MSG(image_loader == nullptr, "Unsupported image type");
Isabella Gottardia4c61882017-11-03 12:11:55 +0000264
Anthony Barbier8a042112018-08-21 18:16:53 +0100265 // Open image file
266 image_loader->open(_filename);
Gian Marco44ec2e72017-10-19 14:13:38 +0100267
Anthony Barbier8a042112018-08-21 18:16:53 +0100268 // Get permutated shape and permutation parameters
269 TensorShape permuted_shape = tensor.info()->tensor_shape();
270 arm_compute::PermutationVector perm;
271 if(tensor.info()->data_layout() != DataLayout::NCHW)
272 {
273 std::tie(permuted_shape, perm) = compute_permutation_parameters(tensor.info()->tensor_shape(), tensor.info()->data_layout());
274 }
Sheri Zhangd80792a2020-11-05 10:43:37 +0000275
276#ifdef __arm__
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100277 ARM_COMPUTE_EXIT_ON_MSG_VAR(image_loader->width() != permuted_shape.x() || image_loader->height() != permuted_shape.y(),
278 "Failed to load image file: dimensions [%d,%d] not correct, expected [%" PRIu32 ",%" PRIu32 "].",
279 image_loader->width(), image_loader->height(), permuted_shape.x(), permuted_shape.y());
Sheri Zhangd80792a2020-11-05 10:43:37 +0000280#else // __arm__
281 ARM_COMPUTE_EXIT_ON_MSG_VAR(image_loader->width() != permuted_shape.x() || image_loader->height() != permuted_shape.y(),
282 "Failed to load image file: dimensions [%d,%d] not correct, expected [%" PRIu64 ",%" PRIu64 "].",
Georgios Pinitas45514032020-12-30 00:03:09 +0000283 image_loader->width(), image_loader->height(),
284 static_cast<uint64_t>(permuted_shape.x()), static_cast<uint64_t>(permuted_shape.y()));
Sheri Zhangd80792a2020-11-05 10:43:37 +0000285#endif // __arm__
Anthony Barbier8a042112018-08-21 18:16:53 +0100286
287 // Fill the tensor with the PPM content (BGR)
288 image_loader->fill_planar_tensor(tensor, _bgr);
289
290 // Preprocess tensor
291 if(_preprocessor)
292 {
293 _preprocessor->preprocess(tensor);
294 }
Georgios Pinitas140fdc72018-02-16 11:42:38 +0000295 }
Gian Marco44ec2e72017-10-19 14:13:38 +0100296
Anthony Barbier8a042112018-08-21 18:16:53 +0100297 _already_loaded = !_already_loaded;
298 return _already_loaded;
Gian Marco44ec2e72017-10-19 14:13:38 +0100299}
300
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100301ValidationInputAccessor::ValidationInputAccessor(const std::string &image_list,
302 std::string images_path,
303 std::unique_ptr<IPreprocessor> preprocessor,
304 bool bgr,
305 unsigned int start,
Anthony Barbier40606df2018-07-23 14:41:59 +0100306 unsigned int end,
307 std::ostream &output_stream)
308 : _path(std::move(images_path)), _images(), _preprocessor(std::move(preprocessor)), _bgr(bgr), _offset(0), _output_stream(output_stream)
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100309{
Anthony Barbier40606df2018-07-23 14:41:59 +0100310 ARM_COMPUTE_EXIT_ON_MSG(start > end, "Invalid validation range!");
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100311
312 std::ifstream ifs;
313 try
314 {
315 ifs.exceptions(std::ifstream::badbit);
316 ifs.open(image_list, std::ios::in | std::ios::binary);
317
318 // Parse image names
319 unsigned int counter = 0;
320 for(std::string line; !std::getline(ifs, line).fail() && counter <= end; ++counter)
321 {
322 // Add image to process if withing range
323 if(counter >= start)
324 {
325 std::stringstream linestream(line);
326 std::string image_name;
327
328 linestream >> image_name;
329 _images.emplace_back(std::move(image_name));
330 }
331 }
332 }
333 catch(const std::ifstream::failure &e)
334 {
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100335 ARM_COMPUTE_ERROR_VAR("Accessing %s: %s", image_list.c_str(), e.what());
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100336 }
337}
338
339bool ValidationInputAccessor::access_tensor(arm_compute::ITensor &tensor)
340{
341 bool ret = _offset < _images.size();
342 if(ret)
343 {
344 utils::JPEGLoader jpeg;
345
346 // Open JPEG file
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100347 std::string image_name = _path + _images[_offset++];
348 jpeg.open(image_name);
Anthony Barbier40606df2018-07-23 14:41:59 +0100349 _output_stream << "[" << _offset << "/" << _images.size() << "] Validating " << image_name << std::endl;
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100350
351 // Get permutated shape and permutation parameters
352 TensorShape permuted_shape = tensor.info()->tensor_shape();
353 arm_compute::PermutationVector perm;
354 if(tensor.info()->data_layout() != DataLayout::NCHW)
355 {
Anthony Barbier4ead11a2018-08-06 09:25:36 +0100356 std::tie(permuted_shape, perm) = compute_permutation_parameters(tensor.info()->tensor_shape(),
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100357 tensor.info()->data_layout());
358 }
Sheri Zhangd80792a2020-11-05 10:43:37 +0000359
360#ifdef __arm__
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100361 ARM_COMPUTE_EXIT_ON_MSG_VAR(jpeg.width() != permuted_shape.x() || jpeg.height() != permuted_shape.y(),
362 "Failed to load image file: dimensions [%d,%d] not correct, expected [%" PRIu32 ",%" PRIu32 "].",
363 jpeg.width(), jpeg.height(), permuted_shape.x(), permuted_shape.y());
Sheri Zhangd80792a2020-11-05 10:43:37 +0000364#else // __arm__
365 ARM_COMPUTE_EXIT_ON_MSG_VAR(jpeg.width() != permuted_shape.x() || jpeg.height() != permuted_shape.y(),
366 "Failed to load image file: dimensions [%d,%d] not correct, expected [%" PRIu64 ",%" PRIu64 "].",
Georgios Pinitas45514032020-12-30 00:03:09 +0000367 jpeg.width(), jpeg.height(),
368 static_cast<uint64_t>(permuted_shape.x()), static_cast<uint64_t>(permuted_shape.y()));
Sheri Zhangd80792a2020-11-05 10:43:37 +0000369#endif // __arm__
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100370
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100371 // Fill the tensor with the JPEG content (BGR)
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100372 jpeg.fill_planar_tensor(tensor, _bgr);
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100373
374 // Preprocess tensor
375 if(_preprocessor)
376 {
377 _preprocessor->preprocess(tensor);
378 }
Georgios Pinitas7c3b9242018-06-21 19:01:25 +0100379 }
380
381 return ret;
382}
383
Georgios Pinitas7908de72018-06-27 12:34:20 +0100384ValidationOutputAccessor::ValidationOutputAccessor(const std::string &image_list,
Georgios Pinitas7908de72018-06-27 12:34:20 +0100385 std::ostream &output_stream,
386 unsigned int start,
387 unsigned int end)
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100388 : _results(), _output_stream(output_stream), _offset(0), _positive_samples_top1(0), _positive_samples_top5(0)
Georgios Pinitas7908de72018-06-27 12:34:20 +0100389{
Anthony Barbier40606df2018-07-23 14:41:59 +0100390 ARM_COMPUTE_EXIT_ON_MSG(start > end, "Invalid validation range!");
Georgios Pinitas7908de72018-06-27 12:34:20 +0100391
392 std::ifstream ifs;
393 try
394 {
395 ifs.exceptions(std::ifstream::badbit);
396 ifs.open(image_list, std::ios::in | std::ios::binary);
397
398 // Parse image correctly classified labels
399 unsigned int counter = 0;
400 for(std::string line; !std::getline(ifs, line).fail() && counter <= end; ++counter)
401 {
402 // Add label if within range
403 if(counter >= start)
404 {
405 std::stringstream linestream(line);
406 std::string image_name;
407 int result;
408
409 linestream >> image_name >> result;
410 _results.emplace_back(result);
411 }
412 }
413 }
414 catch(const std::ifstream::failure &e)
415 {
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100416 ARM_COMPUTE_ERROR_VAR("Accessing %s: %s", image_list.c_str(), e.what());
Georgios Pinitas7908de72018-06-27 12:34:20 +0100417 }
418}
419
420void ValidationOutputAccessor::reset()
421{
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100422 _offset = 0;
423 _positive_samples_top1 = 0;
424 _positive_samples_top5 = 0;
Georgios Pinitas7908de72018-06-27 12:34:20 +0100425}
426
427bool ValidationOutputAccessor::access_tensor(arm_compute::ITensor &tensor)
428{
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100429 bool ret = _offset < _results.size();
430 if(ret)
Georgios Pinitas7908de72018-06-27 12:34:20 +0100431 {
432 // Get results
433 std::vector<size_t> tensor_results;
434 switch(tensor.info()->data_type())
435 {
436 case DataType::QASYMM8:
437 tensor_results = access_predictions_tensor<uint8_t>(tensor);
438 break;
giuros01351bd132019-08-23 14:27:30 +0100439 case DataType::F16:
440 tensor_results = access_predictions_tensor<half>(tensor);
441 break;
Georgios Pinitas7908de72018-06-27 12:34:20 +0100442 case DataType::F32:
443 tensor_results = access_predictions_tensor<float>(tensor);
444 break;
445 default:
446 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
447 }
448
449 // Check if tensor results are within top-n accuracy
450 size_t correct_label = _results[_offset++];
Georgios Pinitas7908de72018-06-27 12:34:20 +0100451
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100452 aggregate_sample(tensor_results, _positive_samples_top1, 1, correct_label);
453 aggregate_sample(tensor_results, _positive_samples_top5, 5, correct_label);
Georgios Pinitas7908de72018-06-27 12:34:20 +0100454 }
455
456 // Report top_n accuracy
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100457 if(_offset >= _results.size())
Georgios Pinitas7908de72018-06-27 12:34:20 +0100458 {
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100459 report_top_n(1, _results.size(), _positive_samples_top1);
460 report_top_n(5, _results.size(), _positive_samples_top5);
Georgios Pinitas7908de72018-06-27 12:34:20 +0100461 }
462
463 return ret;
464}
465
466template <typename T>
467std::vector<size_t> ValidationOutputAccessor::access_predictions_tensor(arm_compute::ITensor &tensor)
468{
469 // Get the predicted class
470 std::vector<size_t> index;
471
472 const auto output_net = reinterpret_cast<T *>(tensor.buffer() + tensor.info()->offset_first_element_in_bytes());
473 const size_t num_classes = tensor.info()->dimension(0);
474
475 index.resize(num_classes);
476
477 // Sort results
478 std::iota(std::begin(index), std::end(index), static_cast<size_t>(0));
479 std::sort(std::begin(index), std::end(index),
480 [&](size_t a, size_t b)
481 {
482 return output_net[a] > output_net[b];
483 });
484
485 return index;
486}
487
Georgios Pinitas12be7ab2018-07-03 12:06:23 +0100488void ValidationOutputAccessor::aggregate_sample(const std::vector<size_t> &res, size_t &positive_samples, size_t top_n, size_t correct_label)
489{
490 auto is_valid_label = [correct_label](size_t label)
491 {
492 return label == correct_label;
493 };
494
495 if(std::any_of(std::begin(res), std::begin(res) + top_n, is_valid_label))
496 {
497 ++positive_samples;
498 }
499}
500
501void ValidationOutputAccessor::report_top_n(size_t top_n, size_t total_samples, size_t positive_samples)
502{
503 size_t negative_samples = total_samples - positive_samples;
504 float accuracy = positive_samples / static_cast<float>(total_samples);
505
506 _output_stream << "----------Top " << top_n << " accuracy ----------" << std::endl
507 << std::endl;
508 _output_stream << "Positive samples : " << positive_samples << std::endl;
509 _output_stream << "Negative samples : " << negative_samples << std::endl;
510 _output_stream << "Accuracy : " << accuracy << std::endl;
511}
512
Isabella Gottardi7234ed82018-11-27 08:51:10 +0000513DetectionOutputAccessor::DetectionOutputAccessor(const std::string &labels_path, std::vector<TensorShape> &imgs_tensor_shapes, std::ostream &output_stream)
514 : _labels(), _tensor_shapes(std::move(imgs_tensor_shapes)), _output_stream(output_stream)
515{
516 _labels.clear();
517
518 std::ifstream ifs;
519
520 try
521 {
522 ifs.exceptions(std::ifstream::badbit);
523 ifs.open(labels_path, std::ios::in | std::ios::binary);
524
525 for(std::string line; !std::getline(ifs, line).fail();)
526 {
527 _labels.emplace_back(line);
528 }
529 }
530 catch(const std::ifstream::failure &e)
531 {
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100532 ARM_COMPUTE_ERROR_VAR("Accessing %s: %s", labels_path.c_str(), e.what());
Isabella Gottardi7234ed82018-11-27 08:51:10 +0000533 }
534}
535
536template <typename T>
537void DetectionOutputAccessor::access_predictions_tensor(ITensor &tensor)
538{
539 const size_t num_detection = tensor.info()->valid_region().shape.y();
540 const auto output_prt = reinterpret_cast<T *>(tensor.buffer() + tensor.info()->offset_first_element_in_bytes());
541
542 if(num_detection > 0)
543 {
544 _output_stream << "---------------------- Detections ----------------------" << std::endl
545 << std::endl;
546
547 _output_stream << std::left << std::setprecision(4) << std::setw(8) << "Image | " << std::setw(8) << "Label | " << std::setw(12) << "Confidence | "
548 << "[ xmin, ymin, xmax, ymax ]" << std::endl;
549
550 for(size_t i = 0; i < num_detection; ++i)
551 {
552 auto im = static_cast<const int>(output_prt[i * 7]);
553 _output_stream << std::setw(8) << im << std::setw(8)
554 << _labels[output_prt[i * 7 + 1]] << std::setw(12) << output_prt[i * 7 + 2]
555 << " [" << (output_prt[i * 7 + 3] * _tensor_shapes[im].x())
556 << ", " << (output_prt[i * 7 + 4] * _tensor_shapes[im].y())
557 << ", " << (output_prt[i * 7 + 5] * _tensor_shapes[im].x())
558 << ", " << (output_prt[i * 7 + 6] * _tensor_shapes[im].y())
559 << "]" << std::endl;
560 }
561 }
562 else
563 {
564 _output_stream << "No detection found." << std::endl;
565 }
566}
567
568bool DetectionOutputAccessor::access_tensor(ITensor &tensor)
569{
570 ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&tensor, 1, DataType::F32);
571
572 switch(tensor.info()->data_type())
573 {
574 case DataType::F32:
575 access_predictions_tensor<float>(tensor);
576 break;
577 default:
578 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
579 }
580
581 return false;
582}
583
Gian Marco44ec2e72017-10-19 14:13:38 +0100584TopNPredictionsAccessor::TopNPredictionsAccessor(const std::string &labels_path, size_t top_n, std::ostream &output_stream)
585 : _labels(), _output_stream(output_stream), _top_n(top_n)
586{
587 _labels.clear();
588
589 std::ifstream ifs;
590
591 try
592 {
593 ifs.exceptions(std::ifstream::badbit);
594 ifs.open(labels_path, std::ios::in | std::ios::binary);
595
596 for(std::string line; !std::getline(ifs, line).fail();)
597 {
598 _labels.emplace_back(line);
599 }
600 }
601 catch(const std::ifstream::failure &e)
602 {
Michalis Spyrou7c60c992019-10-10 14:33:47 +0100603 ARM_COMPUTE_ERROR_VAR("Accessing %s: %s", labels_path.c_str(), e.what());
Gian Marco44ec2e72017-10-19 14:13:38 +0100604 }
605}
606
Giorgio Arenaa66eaa22017-12-21 19:50:06 +0000607template <typename T>
608void TopNPredictionsAccessor::access_predictions_tensor(ITensor &tensor)
Gian Marco44ec2e72017-10-19 14:13:38 +0100609{
Gian Marco44ec2e72017-10-19 14:13:38 +0100610 // Get the predicted class
Giorgio Arenaa66eaa22017-12-21 19:50:06 +0000611 std::vector<T> classes_prob;
Gian Marco44ec2e72017-10-19 14:13:38 +0100612 std::vector<size_t> index;
613
Giorgio Arenaa66eaa22017-12-21 19:50:06 +0000614 const auto output_net = reinterpret_cast<T *>(tensor.buffer() + tensor.info()->offset_first_element_in_bytes());
Gian Marco44ec2e72017-10-19 14:13:38 +0100615 const size_t num_classes = tensor.info()->dimension(0);
616
617 classes_prob.resize(num_classes);
618 index.resize(num_classes);
619
620 std::copy(output_net, output_net + num_classes, classes_prob.begin());
621
622 // Sort results
623 std::iota(std::begin(index), std::end(index), static_cast<size_t>(0));
624 std::sort(std::begin(index), std::end(index),
625 [&](size_t a, size_t b)
626 {
627 return classes_prob[a] > classes_prob[b];
628 });
629
630 _output_stream << "---------- Top " << _top_n << " predictions ----------" << std::endl
631 << std::endl;
632 for(size_t i = 0; i < _top_n; ++i)
633 {
634 _output_stream << std::fixed << std::setprecision(4)
Giorgio Arenaa66eaa22017-12-21 19:50:06 +0000635 << +classes_prob[index.at(i)]
Gian Marco44ec2e72017-10-19 14:13:38 +0100636 << " - [id = " << index.at(i) << "]"
637 << ", " << _labels[index.at(i)] << std::endl;
638 }
Giorgio Arenaa66eaa22017-12-21 19:50:06 +0000639}
640
641bool TopNPredictionsAccessor::access_tensor(ITensor &tensor)
642{
643 ARM_COMPUTE_ERROR_ON_DATA_TYPE_CHANNEL_NOT_IN(&tensor, 1, DataType::F32, DataType::QASYMM8);
644 ARM_COMPUTE_ERROR_ON(_labels.size() != tensor.info()->dimension(0));
645
646 switch(tensor.info()->data_type())
647 {
648 case DataType::QASYMM8:
649 access_predictions_tensor<uint8_t>(tensor);
650 break;
651 case DataType::F32:
652 access_predictions_tensor<float>(tensor);
653 break;
654 default:
655 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
656 }
Gian Marco44ec2e72017-10-19 14:13:38 +0100657
658 return false;
659}
660
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100661RandomAccessor::RandomAccessor(PixelValue lower, PixelValue upper, std::random_device::result_type seed)
662 : _lower(lower), _upper(upper), _seed(seed)
663{
664}
665
666template <typename T, typename D>
667void RandomAccessor::fill(ITensor &tensor, D &&distribution)
668{
669 std::mt19937 gen(_seed);
670
hakanardof36ac352018-02-16 10:06:34 +0100671 if(tensor.info()->padding().empty() && (dynamic_cast<SubTensor *>(&tensor) == nullptr))
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100672 {
673 for(size_t offset = 0; offset < tensor.info()->total_size(); offset += tensor.info()->element_size())
674 {
Michele Di Giorgio88731f02018-09-25 16:49:27 +0100675 const auto value = static_cast<T>(distribution(gen));
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100676 *reinterpret_cast<T *>(tensor.buffer() + offset) = value;
677 }
678 }
679 else
680 {
681 // If tensor has padding accessing tensor elements through execution window.
682 Window window;
683 window.use_tensor_dimensions(tensor.info()->tensor_shape());
684
685 execute_window_loop(window, [&](const Coordinates & id)
686 {
Michele Di Giorgio88731f02018-09-25 16:49:27 +0100687 const auto value = static_cast<T>(distribution(gen));
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100688 *reinterpret_cast<T *>(tensor.ptr_to_element(id)) = value;
689 });
690 }
691}
692
693bool RandomAccessor::access_tensor(ITensor &tensor)
694{
695 switch(tensor.info()->data_type())
696 {
John Kesapidesfb68ca12019-01-21 14:13:27 +0000697 case DataType::QASYMM8:
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100698 case DataType::U8:
699 {
700 std::uniform_int_distribution<uint8_t> distribution_u8(_lower.get<uint8_t>(), _upper.get<uint8_t>());
701 fill<uint8_t>(tensor, distribution_u8);
702 break;
703 }
704 case DataType::S8:
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100705 {
706 std::uniform_int_distribution<int8_t> distribution_s8(_lower.get<int8_t>(), _upper.get<int8_t>());
707 fill<int8_t>(tensor, distribution_s8);
708 break;
709 }
710 case DataType::U16:
711 {
712 std::uniform_int_distribution<uint16_t> distribution_u16(_lower.get<uint16_t>(), _upper.get<uint16_t>());
713 fill<uint16_t>(tensor, distribution_u16);
714 break;
715 }
716 case DataType::S16:
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100717 {
718 std::uniform_int_distribution<int16_t> distribution_s16(_lower.get<int16_t>(), _upper.get<int16_t>());
719 fill<int16_t>(tensor, distribution_s16);
720 break;
721 }
722 case DataType::U32:
723 {
724 std::uniform_int_distribution<uint32_t> distribution_u32(_lower.get<uint32_t>(), _upper.get<uint32_t>());
725 fill<uint32_t>(tensor, distribution_u32);
726 break;
727 }
728 case DataType::S32:
729 {
730 std::uniform_int_distribution<int32_t> distribution_s32(_lower.get<int32_t>(), _upper.get<int32_t>());
731 fill<int32_t>(tensor, distribution_s32);
732 break;
733 }
734 case DataType::U64:
735 {
736 std::uniform_int_distribution<uint64_t> distribution_u64(_lower.get<uint64_t>(), _upper.get<uint64_t>());
737 fill<uint64_t>(tensor, distribution_u64);
738 break;
739 }
740 case DataType::S64:
741 {
742 std::uniform_int_distribution<int64_t> distribution_s64(_lower.get<int64_t>(), _upper.get<int64_t>());
743 fill<int64_t>(tensor, distribution_s64);
744 break;
745 }
746 case DataType::F16:
747 {
Giorgio Arena33b103b2021-01-08 10:37:15 +0000748 arm_compute::utils::uniform_real_distribution_16bit<half> distribution_f16(_lower.get<float>(), _upper.get<float>());
Michele Di Giorgio88731f02018-09-25 16:49:27 +0100749 fill<half>(tensor, distribution_f16);
Michalis Spyrou53b405f2017-09-27 15:55:31 +0100750 break;
751 }
752 case DataType::F32:
753 {
754 std::uniform_real_distribution<float> distribution_f32(_lower.get<float>(), _upper.get<float>());
755 fill<float>(tensor, distribution_f32);
756 break;
757 }
758 case DataType::F64:
759 {
760 std::uniform_real_distribution<double> distribution_f64(_lower.get<double>(), _upper.get<double>());
761 fill<double>(tensor, distribution_f64);
762 break;
763 }
764 default:
765 ARM_COMPUTE_ERROR("NOT SUPPORTED!");
766 }
767 return true;
768}
769
Georgios Pinitascac13b12018-04-27 19:07:19 +0100770NumPyBinLoader::NumPyBinLoader(std::string filename, DataLayout file_layout)
Anthony Barbier8a042112018-08-21 18:16:53 +0100771 : _already_loaded(false), _filename(std::move(filename)), _file_layout(file_layout)
Anthony Barbier2a07e182017-08-04 18:20:27 +0100772{
773}
774
775bool NumPyBinLoader::access_tensor(ITensor &tensor)
776{
Anthony Barbier8a042112018-08-21 18:16:53 +0100777 if(!_already_loaded)
778 {
779 utils::NPYLoader loader;
780 loader.open(_filename, _file_layout);
781 loader.fill_tensor(tensor);
782 }
Anthony Barbier2a07e182017-08-04 18:20:27 +0100783
Anthony Barbier8a042112018-08-21 18:16:53 +0100784 _already_loaded = !_already_loaded;
785 return _already_loaded;
Anthony Barbier87f21cd2017-11-10 16:27:32 +0000786}