#ifndef CAFFE2_OPERATORS_INT8_AVERAGE_POOL_OP_H_
#define CAFFE2_OPERATORS_INT8_AVERAGE_POOL_OP_H_
+#include <qnnpack.h>
+
#include "caffe2/core/context.h"
#include "caffe2/core/operator.h"
#include "caffe2/core/tensor_int8.h"
#include "caffe2/operators/conv_pool_op_base.h"
-#include "caffe2/operators/quantized/int8_simd.h"
#include "caffe2/operators/quantized/int8_utils.h"
namespace caffe2 {
namespace int8 {
-namespace {
-
-/*
- * Implementation based on TensorFlow Lite kernels:
- * - Repo: https://github.com/tensorflow/tensorflow
- * - Path: tensorflow/contrib/lite/kernels/internal/optimized/optimized_ops.h
- * - Hash: d4ad9c73969c45d1a224ebfc43eb645b9860216b
- */
-
-/* Copyright 2017 The TensorFlow Authors. All Rights Reserved.
-Licensed under the Apache License, Version 2.0 (the "License");
-you may not use this file except in compliance with the License.
-You may obtain a copy of the License at
- http://www.apache.org/licenses/LICENSE-2.0
-Unless required by applicable law or agreed to in writing, software
-distributed under the License is distributed on an "AS IS" BASIS,
-WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
-See the License for the specific language governing permissions and
-limitations under the License.
-==============================================================================*/
-
-void Int8AveragePool(
- const uint8_t* input_data,
- at::IntList input_dims,
- int stride_width,
- int stride_height,
- int pad_width,
- int pad_height,
- int filter_width,
- int filter_height,
- uint8_t* output_data,
- at::IntList output_dims,
- uint8_t output_activation_min,
- uint8_t output_activation_max) {
- DCHECK_LE(output_activation_min, output_activation_max);
- const int batches = input_dims[0];
- const int depth = input_dims[3];
- const int input_height = input_dims[1];
- const int input_width = input_dims[2];
- const int output_height = output_dims[1];
- const int output_width = output_dims[2];
- for (int batch = 0; batch < batches; ++batch) {
- for (int out_y = 0; out_y < output_height; ++out_y) {
- for (int out_x = 0; out_x < output_width; ++out_x) {
- const int in_x_origin = (out_x * stride_width) - pad_width;
- const int in_y_origin = (out_y * stride_height) - pad_height;
- const int filter_x_start = std::max(0, -in_x_origin);
- const int filter_x_end =
- std::min(filter_width, input_width - in_x_origin);
- const int filter_y_start = std::max(0, -in_y_origin);
- const int filter_y_end =
- std::min(filter_height, input_height - in_y_origin);
- const int filter_count =
- (filter_x_end - filter_x_start) * (filter_y_end - filter_y_start);
- // 1280 required by Inception v3
- static constexpr int kAccBufferMaxSize = 2048;
- DCHECK_LE(depth, kAccBufferMaxSize);
- uint16_t acc[kAccBufferMaxSize];
- memset(acc, 0, depth * sizeof(acc[0]));
- const uint8_t* input_ptr = input_data + depth * in_x_origin +
- depth * input_width * in_y_origin +
- depth * input_width * input_height * batch;
- for (int fy = filter_y_start; fy < filter_y_end; fy++) {
- const uint8_t* input_row_ptr =
- input_ptr + fy * input_width * depth + filter_x_start * depth;
- for (int fx = filter_x_start; fx < filter_x_end; fx++) {
- int channel = 0;
-#ifdef INT8_NEON_SIMD
- for (; channel <= depth - 16; channel += 16) {
- uint16x8_t acc_reg[2];
- for (int i = 0; i < 2; i++) {
- acc_reg[i] = vld1q_u16(acc + channel + 8 * i);
- }
- uint8x16_t input_reg = vld1q_u8(input_row_ptr);
- input_row_ptr += 16;
- acc_reg[0] = vaddw_u8(acc_reg[0], vget_low_u8(input_reg));
- acc_reg[1] = vaddw_u8(acc_reg[1], vget_high_u8(input_reg));
- for (int i = 0; i < 2; i++) {
- vst1q_u16(acc + channel + 8 * i, acc_reg[i]);
- }
- }
- for (; channel <= depth - 8; channel += 8) {
- uint16x8_t acc_reg = vld1q_u16(acc + channel);
- uint8x8_t input_reg = vld1_u8(input_row_ptr);
- input_row_ptr += 8;
- acc_reg = vaddw_u8(acc_reg, input_reg);
- vst1q_u16(acc + channel, acc_reg);
- }
-#endif
- for (; channel < depth; ++channel) {
- acc[channel] += *input_row_ptr++;
- }
- }
- }
- uint8_t* output_ptr = output_data + out_x * depth +
- out_y * depth * output_width +
- batch * depth * output_width * output_height;
- int channel = 0;
-#ifdef INT8_NEON_SIMD
-#define AVGPOOL_DIVIDING_BY(FILTER_COUNT) \
- if (filter_count == FILTER_COUNT) { \
- for (; channel <= depth - 8; channel += 8) { \
- uint16_t buf[8]; \
- for (int i = 0; i < 8; i++) { \
- buf[i] = (acc[channel + i] + FILTER_COUNT / 2) / FILTER_COUNT; \
- } \
- uint8x8_t buf8 = vqmovn_u16(vld1q_u16(buf)); \
- buf8 = vmin_u8(buf8, vdup_n_u8(output_activation_max)); \
- buf8 = vmax_u8(buf8, vdup_n_u8(output_activation_min)); \
- vst1_u8(output_ptr + channel, buf8); \
- } \
- }
- AVGPOOL_DIVIDING_BY(9)
- AVGPOOL_DIVIDING_BY(15)
-#undef AVGPOOL_DIVIDING_BY
- for (; channel <= depth - 8; channel += 8) {
- uint16_t buf[8];
- for (int i = 0; i < 8; i++) {
- buf[i] = (acc[channel + i] + filter_count / 2) / filter_count;
- }
- uint8x8_t buf8 = vqmovn_u16(vld1q_u16(buf));
- buf8 = vmin_u8(buf8, vdup_n_u8(output_activation_max));
- buf8 = vmax_u8(buf8, vdup_n_u8(output_activation_min));
- vst1_u8(output_ptr + channel, buf8);
- }
-#endif
- for (; channel < depth; ++channel) {
- uint16_t a = (acc[channel] + filter_count / 2) / filter_count;
- a = std::max<uint16_t>(a, output_activation_min);
- a = std::min<uint16_t>(a, output_activation_max);
- output_ptr[channel] = static_cast<uint8_t>(a);
- }
- }
- }
- }
-}
-
-} // namespace
-
template <Activation Ac>
class Int8AveragePoolOp final : public ConvPoolOpBase<CPUContext> {
public:
Int8AveragePoolOp(const OperatorDef& operator_def, Workspace* ws)
: ConvPoolOpBase<CPUContext>(operator_def, ws) {
OPERATOR_NEEDS_FEATURE(
- this->order_ == StorageOrder::NHWC, "Int8 only supports NCHW order.");
+ this->order_ == StorageOrder::NHWC, "Int8 only supports NHWC order.");
+ }
+
+ ~Int8AveragePoolOp() {
+ if (this->qnnpackOperator_ != nullptr) {
+ qnnp_delete_operator(this->qnnpackOperator_);
+ this->qnnpackOperator_ = nullptr;
+ }
+ if (this->qnnpackGlobalOperator_ != nullptr) {
+ qnnp_delete_operator(this->qnnpackGlobalOperator_);
+ this->qnnpackGlobalOperator_ = nullptr;
+ }
}
bool RunOnDeviceWithOrderNHWC() override {
const auto& X = Inputs()[0]->template Get<Int8TensorCPU>();
auto* Y = Outputs()[0]->template GetMutable<Int8TensorCPU>();
- Y->scale = X.scale;
- Y->zero_point = X.zero_point;
- int32_t Y_offset = this->template GetSingleArgument<int>("Y_zero_point", 0);
+ int32_t Y_zero_point =
+ this->template GetSingleArgument<int>("Y_zero_point", 0);
auto Y_scale = this->template GetSingleArgument<float>("Y_scale", 1);
- CHECK_EQ(Y_offset, X.zero_point);
- CHECK_EQ(Y_scale, X.scale);
+ Y->scale = Y_scale;
+ Y->zero_point = Y_zero_point;
CHECK_EQ(X.t.dim(), 4);
const int channels = X.t.dim32(3);
ConvPoolOpBase<CPUContext>::SetOutputSize(X.t, &(Y->t), channels);
- Int8AveragePool(
- X.t.template data<uint8_t>(),
- X.t.sizes(),
- stride_w(),
- stride_h(),
- pad_l(),
- pad_t(),
- kernel_w(),
- kernel_h(),
- Y->t.template mutable_data<uint8_t>(),
- Y->t.sizes(),
- activationLimits(Y->scale, Y->zero_point, Ac).first,
- activationLimits(Y->scale, Y->zero_point, Ac).second);
+ initQNNPACK();
+
+ const bool anyPadding =
+ pad_t() != 0 || pad_r() != 0 || pad_b() != 0 || pad_l() != 0;
+ const bool anyStride = stride_h() > 1 || stride_w() > 1;
+ const bool globalPooling = !anyPadding && !anyStride &&
+ (X.t.dim32(1) == kernel_h() && X.t.dim32(2) == kernel_w());
+ if (globalPooling) {
+ if (this->qnnpackGlobalOperator_ == nullptr) {
+ const qnnp_status createStatus =
+ qnnp_create_global_average_pooling_nwc_q8(
+ channels,
+ X.zero_point, X.scale,
+ Y->zero_point, Y->scale,
+ activationLimits(Y->scale, Y->zero_point, Ac).first,
+ activationLimits(Y->scale, Y->zero_point, Ac).second,
+ &this->qnnpackGlobalOperator_);
+ CAFFE_ENFORCE(
+ createStatus == qnnp_status_success,
+ "failed to create QNNPACK Global Average Pooling operator");
+ CAFFE_ENFORCE(this->qnnpackGlobalOperator_ != nullptr);
+ }
+
+ const qnnp_status setupStatus = qnnp_setup_global_average_pooling_nwc_q8(
+ this->qnnpackGlobalOperator_,
+ X.t.dim32(0), X.t.dim32(1) * X.t.dim32(2),
+ X.t.template data<uint8_t>(), channels,
+ Y->t.template mutable_data<uint8_t>(), channels);
+ CAFFE_ENFORCE(
+ setupStatus == qnnp_status_success,
+ "failed to setup QNNPACK Global Average Pooling operator");
+
+#ifdef FBCODE_CAFFE2
+ const qnnp_status runStatus =
+ qnnp_run_operator(this->qnnpackGlobalOperator_,
+ nullptr /* thread pool */);
+#else
+ pthreadpool_t threadpool =
+ reinterpret_cast<pthreadpool_t>(ws_->GetThreadPool());
+ const qnnp_status runStatus =
+ qnnp_run_operator(this->qnnpackGlobalOperator_, threadpool);
+#endif
+ CAFFE_ENFORCE(
+ runStatus == qnnp_status_success,
+ "failed to run QNNPACK Global Average Pooling operator");
+ } else {
+ if (this->qnnpackOperator_ == nullptr) {
+ const qnnp_status createStatus = qnnp_create_average_pooling2d_nhwc_q8(
+ pad_t(), pad_r(), pad_b(), pad_l(),
+ kernel_h(), kernel_w(),
+ stride_h(), stride_w(),
+ channels,
+ X.zero_point, X.scale,
+ Y->zero_point, Y->scale,
+ activationLimits(Y->scale, Y->zero_point, Ac).first,
+ activationLimits(Y->scale, Y->zero_point, Ac).second,
+ &this->qnnpackOperator_);
+ CAFFE_ENFORCE(
+ createStatus == qnnp_status_success,
+ "failed to create QNNPACK Average Pooling operator");
+ CAFFE_ENFORCE(this->qnnpackOperator_ != nullptr);
+ }
+
+ const qnnp_status setupStatus = qnnp_setup_average_pooling2d_nhwc_q8(
+ this->qnnpackOperator_,
+ X.t.dim32(0), X.t.dim32(1), X.t.dim32(2),
+ X.t.template data<uint8_t>(), channels,
+ Y->t.template mutable_data<uint8_t>(), channels,
+ nullptr /* thread pool */);
+ CAFFE_ENFORCE(
+ setupStatus == qnnp_status_success,
+ "failed to setup QNNPACK Average Pooling operator");
+
+#ifdef FBCODE_CAFFE2
+ const qnnp_status runStatus =
+ qnnp_run_operator(this->qnnpackOperator_, nullptr /* thread pool */);
+#else
+ pthreadpool_t threadpool =
+ reinterpret_cast<pthreadpool_t>(ws_->GetThreadPool());
+ const qnnp_status runStatus =
+ qnnp_run_operator(this->qnnpackOperator_, threadpool);
+#endif
+ CAFFE_ENFORCE(
+ runStatus == qnnp_status_success,
+ "failed to run QNNPACK Average Pooling operator");
+ }
+
return true;
}
+ private:
+ // QNNPACK Average Pooling operator
+ qnnp_operator_t qnnpackOperator_{nullptr};
+ // QNNPACK Global Average Pooling operator
+ qnnp_operator_t qnnpackGlobalOperator_{nullptr};
};
} // namespace int8