/*
 * Copyright (c) 2021-2024 Arm Limited.
 *
 * SPDX-License-Identifier: MIT
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to
 * deal in the Software without restriction, including without limitation the
 * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
 * sell copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in all
 * copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
 * SOFTWARE.
 */

#include "pooling.hpp"
#include <cstdint>
#include <cstddef>
#include <cstring>
#include <cmath>


#if defined(ARM_COMPUTE_ENABLE_SVE)

namespace arm_conv {
namespace pooling {

namespace {
  struct RescaleParams
  {
    int32_t multiplier, shift;
  };

  constexpr RescaleParams rescale_params[8] = {
    {0x40000000, -0},  // 1/2
    {0x55555556, -1},  // 1/3
    {0x40000000, -1},  // 1/4
    {0x66666666, -2},  // 1/5
    {0x55555556, -2},  // 1/6
    {0x49249249, -2},  // 1/7
    {0x40000000, -2},  // 1/8
    {0x71c71c72, -3},  // 1/9
  };
}

void sve_u8q_nhwc_avg_generic_depthfirst_impl(
  const uint64_t window_cells,
  const uint64_t n_valid_cells,
  uint64_t n_channels,
  const uint8_t *const *const inptrs,
  uint8_t *outptr,
  const Requantize32 &qp
)
{
  if (n_valid_cells == 1 && window_cells == 1)
  {
    // In this case, simply copy from the input to the output
    std::memcpy(outptr, *inptrs, n_channels);
    return;
  }

  // Compute (or look up) the rescale values
  int32_t shift_value = 0, rescale_value = 0;
  if (2 <= window_cells && window_cells <= 9)
  {
    auto &params = rescale_params[window_cells - 2];
    rescale_value = params.multiplier;
    shift_value = params.shift;
  }
  else
  {
    auto f_rescale_value = 1.0f / static_cast<float>(window_cells);

    shift_value = 0;
    while (f_rescale_value < 0.5f)
    {
      shift_value--;
      f_rescale_value *= 2.0f;
    }

    int64_t long_rescale_value = round(f_rescale_value * static_cast<float>(1ll << 31));
    if (long_rescale_value == (1ll << 31))
    {
      shift_value++;
      long_rescale_value >>= 1;
    }
    rescale_value = static_cast<int32_t>(long_rescale_value);
  }


  // Initialise the accumulators such that the offsets are subtracted for all
  // valid inputs.
  const int32_t accumulator_init = -qp.input_offset * n_valid_cells;

  // Combine together the rescale value for the requantization and the scaling
  // factor for the average pool.
  const int32_t shift = qp.per_layer_left_shift - qp.per_layer_right_shift + shift_value;
  const int32_t left_shift = shift > 0 ? shift : 0;
  const int32_t right_shift = shift <= 0 ? shift : 0;

  int32_t combined_rescale_value = 0;
  __asm__ __volatile__ (
      "mov v16.s[0], %w[per_layer_mul]\n"
      "mov v17.s[0], %w[rescale_value]\n"
      "sqrdmulh s18, s16, s17\n"
      "mov %w[combined_rescale_value], v18.s[0]\n"
    : [combined_rescale_value] "=r" (combined_rescale_value)
    : [per_layer_mul] "r" (qp.per_layer_mul), [rescale_value] "r" (rescale_value)
    : "v16", "v17", "v18"
  );

  __asm__ __volatile__(
    "mov x27, #0x0\n"
    "cntb x26\n"
    "cntb x25, ALL, MUL #2\n"
    "cntb x24, ALL, MUL #3\n"
    "ptrue p4.b\n"
    "whilelt p3.b, x27, %x[n_channels]\n"
    "whilelt p2.b, x26, %x[n_channels]\n"
    "whilelt p1.b, x25, %x[n_channels]\n"
    "whilelt p0.b, x24, %x[n_channels]\n"
    "b.none 7f\n"
    "1:"  // 4-vectors of channels
    "ld1rw { z15.s }, p4/Z, [%x[accumulator_init]]\n"
    "lsr x23, %x[n_valid_cells], #0x1\n"
    "mov x22, %x[inptrs]\n"
    "mov z14.d, z15.d\n"
    "mov z13.d, z15.d\n"
    "mov z12.d, z15.d\n"
    "mov z11.d, z15.d\n"
    "mov z10.d, z15.d\n"
    "mov z9.d, z15.d\n"
    "mov z8.d, z15.d\n"
    "mov z7.d, z15.d\n"
    "mov z6.d, z15.d\n"
    "mov z5.d, z15.d\n"
    "mov z4.d, z15.d\n"
    "mov z3.d, z15.d\n"
    "mov z2.d, z15.d\n"
    "mov z1.d, z15.d\n"
    "mov z0.d, z15.d\n"
    "cbz x23, 4f\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    "ld1b { z31.b }, p3/Z, [x21, x27]\n"
    "ld1b { z30.b }, p3/Z, [x20, x27]\n"
    "ld1b { z29.b }, p2/Z, [x21, x26]\n"
    "ld1b { z28.b }, p2/Z, [x20, x26]\n"
    "ld1b { z27.b }, p1/Z, [x21, x25]\n"
    "ld1b { z26.b }, p1/Z, [x20, x25]\n"
    "ld1b { z25.b }, p0/Z, [x21, x24]\n"
    "ld1b { z24.b }, p0/Z, [x20, x24]\n"
    "beq 3f\n"
    "2:"  // 4-vectors of channels: 2 inputs loop
    ".inst 0x455e0bf7  // uaddlb z23.h, z31.b, z30.b\n"
    ".inst 0x455e0ff6  // uaddlt z22.h, z31.b, z30.b\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    ".inst 0x455c0bb5  // uaddlb z21.h, z29.b, z28.b\n"
    ".inst 0x455c0fb4  // uaddlt z20.h, z29.b, z28.b\n"
    "add x22, x22, #0x10\n"
    ".inst 0x455a0b73  // uaddlb z19.h, z27.b, z26.b\n"
    ".inst 0x455a0f72  // uaddlt z18.h, z27.b, z26.b\n"
    ".inst 0x45580b31  // uaddlb z17.h, z25.b, z24.b\n"
    ".inst 0x45580f30  // uaddlt z16.h, z25.b, z24.b\n"
    "ld1b { z31.b }, p3/Z, [x21, x27]\n"
    "ld1b { z30.b }, p3/Z, [x20, x27]\n"
    ".inst 0x459749ef  // uaddwb z15.s, z15.s, z23.h\n"
    ".inst 0x45974dce  // uaddwt z14.s, z14.s, z23.h\n"
    "ld1b { z29.b }, p2/Z, [x21, x26]\n"
    "ld1b { z28.b }, p2/Z, [x20, x26]\n"
    ".inst 0x459649ad  // uaddwb z13.s, z13.s, z22.h\n"
    ".inst 0x45964d8c  // uaddwt z12.s, z12.s, z22.h\n"
    "ld1b { z27.b }, p1/Z, [x21, x25]\n"
    "ld1b { z26.b }, p1/Z, [x20, x25]\n"
    ".inst 0x4595496b  // uaddwb z11.s, z11.s, z21.h\n"
    ".inst 0x45954d4a  // uaddwt z10.s, z10.s, z21.h\n"
    "ld1b { z25.b }, p0/Z, [x21, x24]\n"
    "ld1b { z24.b }, p0/Z, [x20, x24]\n"
    ".inst 0x45944929  // uaddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944d08  // uaddwt z8.s, z8.s, z20.h\n"
    ".inst 0x459348e7  // uaddwb z7.s, z7.s, z19.h\n"
    ".inst 0x45934cc6  // uaddwt z6.s, z6.s, z19.h\n"
    ".inst 0x459248a5  // uaddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924c84  // uaddwt z4.s, z4.s, z18.h\n"
    ".inst 0x45914863  // uaddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914c42  // uaddwt z2.s, z2.s, z17.h\n"
    ".inst 0x45904821  // uaddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904c00  // uaddwt z0.s, z0.s, z16.h\n"
    "bgt 2b\n"
    "3:"  // 4-vectors of channels: 2 inputs tail
    ".inst 0x455e0bf7  // uaddlb z23.h, z31.b, z30.b\n"
    ".inst 0x455e0ff6  // uaddlt z22.h, z31.b, z30.b\n"
    ".inst 0x455c0bb5  // uaddlb z21.h, z29.b, z28.b\n"
    ".inst 0x455c0fb4  // uaddlt z20.h, z29.b, z28.b\n"
    ".inst 0x455a0b73  // uaddlb z19.h, z27.b, z26.b\n"
    ".inst 0x455a0f72  // uaddlt z18.h, z27.b, z26.b\n"
    ".inst 0x45580b31  // uaddlb z17.h, z25.b, z24.b\n"
    ".inst 0x45580f30  // uaddlt z16.h, z25.b, z24.b\n"
    ".inst 0x459749ef  // uaddwb z15.s, z15.s, z23.h\n"
    ".inst 0x45974dce  // uaddwt z14.s, z14.s, z23.h\n"
    ".inst 0x459649ad  // uaddwb z13.s, z13.s, z22.h\n"
    ".inst 0x45964d8c  // uaddwt z12.s, z12.s, z22.h\n"
    ".inst 0x4595496b  // uaddwb z11.s, z11.s, z21.h\n"
    ".inst 0x45954d4a  // uaddwt z10.s, z10.s, z21.h\n"
    ".inst 0x45944929  // uaddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944d08  // uaddwt z8.s, z8.s, z20.h\n"
    ".inst 0x459348e7  // uaddwb z7.s, z7.s, z19.h\n"
    ".inst 0x45934cc6  // uaddwt z6.s, z6.s, z19.h\n"
    ".inst 0x459248a5  // uaddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924c84  // uaddwt z4.s, z4.s, z18.h\n"
    ".inst 0x45914863  // uaddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914c42  // uaddwt z2.s, z2.s, z17.h\n"
    ".inst 0x45904821  // uaddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904c00  // uaddwt z0.s, z0.s, z16.h\n"
    "4:"  // 4-vectors of channels: After loop
    "ands x21, %x[n_valid_cells], #0x1\n"
    "beq 6f\n"
    "5:"  // 4-vectors of channels: Single input loop
    "ldr x20, [x22], #0x8\n"
    "subs x21, x21, #0x1\n"
    "ld1b { z19.b }, p3/Z, [x20, x27]\n"
    "ld1b { z18.b }, p2/Z, [x20, x26]\n"
    "ld1b { z17.b }, p1/Z, [x20, x25]\n"
    "ld1b { z16.b }, p0/Z, [x20, x24]\n"
    ".inst 0x4508aa77  // ushllb z23.h, z19.b, #0x0\n"
    ".inst 0x4508ae76  // ushllt z22.h, z19.b, #0x0\n"
    ".inst 0x4508aa55  // ushllb z21.h, z18.b, #0x0\n"
    ".inst 0x4508ae54  // ushllt z20.h, z18.b, #0x0\n"
    ".inst 0x4508aa33  // ushllb z19.h, z17.b, #0x0\n"
    ".inst 0x4508ae32  // ushllt z18.h, z17.b, #0x0\n"
    ".inst 0x4508aa11  // ushllb z17.h, z16.b, #0x0\n"
    ".inst 0x4508ae10  // ushllt z16.h, z16.b, #0x0\n"
    ".inst 0x459749ef  // uaddwb z15.s, z15.s, z23.h\n"
    ".inst 0x45974dce  // uaddwt z14.s, z14.s, z23.h\n"
    ".inst 0x459649ad  // uaddwb z13.s, z13.s, z22.h\n"
    ".inst 0x45964d8c  // uaddwt z12.s, z12.s, z22.h\n"
    ".inst 0x4595496b  // uaddwb z11.s, z11.s, z21.h\n"
    ".inst 0x45954d4a  // uaddwt z10.s, z10.s, z21.h\n"
    ".inst 0x45944929  // uaddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944d08  // uaddwt z8.s, z8.s, z20.h\n"
    ".inst 0x459348e7  // uaddwb z7.s, z7.s, z19.h\n"
    ".inst 0x45934cc6  // uaddwt z6.s, z6.s, z19.h\n"
    ".inst 0x459248a5  // uaddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924c84  // uaddwt z4.s, z4.s, z18.h\n"
    ".inst 0x45914863  // uaddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914c42  // uaddwt z2.s, z2.s, z17.h\n"
    ".inst 0x45904821  // uaddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904c00  // uaddwt z0.s, z0.s, z16.h\n"
    "bgt 5b\n"
    "6:"  // 4-vectors of channels: Single input loop: End
    "ld1rw { z21.s }, p4/Z, [%x[left_shift]]\n"
    "ld1rw { z19.s }, p4/Z, [%x[combined_rescale_value]]\n"
    "add x20, %x[quant_params], %[offsetof_qp_output_offset]\n"
    "mov z18.s, #0x0\n"
    "ld1rw { z17.s }, p4/Z, [%x[right_shift]]\n"
    "ld1rw { z16.s }, p4/Z, [x20]\n"
    "mov z20.s, #0xff\n"
    ".inst 0x448292af  // srshl z15.s, p4/M, z15.s, z21.s\n"
    ".inst 0x448292ae  // srshl z14.s, p4/M, z14.s, z21.s\n"
    ".inst 0x448292ad  // srshl z13.s, p4/M, z13.s, z21.s\n"
    ".inst 0x448292ac  // srshl z12.s, p4/M, z12.s, z21.s\n"
    ".inst 0x448292ab  // srshl z11.s, p4/M, z11.s, z21.s\n"
    ".inst 0x448292aa  // srshl z10.s, p4/M, z10.s, z21.s\n"
    ".inst 0x448292a9  // srshl z9.s, p4/M, z9.s, z21.s\n"
    ".inst 0x04b375ef  // sqrdmulh z15.s, z15.s, z19.s\n"
    ".inst 0x448292a8  // srshl z8.s, p4/M, z8.s, z21.s\n"
    ".inst 0x448292a7  // srshl z7.s, p4/M, z7.s, z21.s\n"
    ".inst 0x04b375ce  // sqrdmulh z14.s, z14.s, z19.s\n"
    ".inst 0x04b375ad  // sqrdmulh z13.s, z13.s, z19.s\n"
    ".inst 0x448292a6  // srshl z6.s, p4/M, z6.s, z21.s\n"
    ".inst 0x448292a5  // srshl z5.s, p4/M, z5.s, z21.s\n"
    ".inst 0x04b3758c  // sqrdmulh z12.s, z12.s, z19.s\n"
    ".inst 0x04b3756b  // sqrdmulh z11.s, z11.s, z19.s\n"
    ".inst 0x448292a4  // srshl z4.s, p4/M, z4.s, z21.s\n"
    ".inst 0x448292a3  // srshl z3.s, p4/M, z3.s, z21.s\n"
    ".inst 0x04b3754a  // sqrdmulh z10.s, z10.s, z19.s\n"
    ".inst 0x04b37529  // sqrdmulh z9.s, z9.s, z19.s\n"
    ".inst 0x448292a2  // srshl z2.s, p4/M, z2.s, z21.s\n"
    ".inst 0x448292a1  // srshl z1.s, p4/M, z1.s, z21.s\n"
    ".inst 0x04b37508  // sqrdmulh z8.s, z8.s, z19.s\n"
    ".inst 0x04b374e7  // sqrdmulh z7.s, z7.s, z19.s\n"
    ".inst 0x448292a0  // srshl z0.s, p4/M, z0.s, z21.s\n"
    ".inst 0x04b374c6  // sqrdmulh z6.s, z6.s, z19.s\n"
    ".inst 0x04b374a5  // sqrdmulh z5.s, z5.s, z19.s\n"
    ".inst 0x4482922f  // srshl z15.s, p4/M, z15.s, z17.s\n"
    ".inst 0x04b37484  // sqrdmulh z4.s, z4.s, z19.s\n"
    ".inst 0x04b37463  // sqrdmulh z3.s, z3.s, z19.s\n"
    ".inst 0x4482922e  // srshl z14.s, p4/M, z14.s, z17.s\n"
    ".inst 0x4482922d  // srshl z13.s, p4/M, z13.s, z17.s\n"
    ".inst 0x04b37442  // sqrdmulh z2.s, z2.s, z19.s\n"
    ".inst 0x04b37421  // sqrdmulh z1.s, z1.s, z19.s\n"
    ".inst 0x4482922c  // srshl z12.s, p4/M, z12.s, z17.s\n"
    ".inst 0x4482922b  // srshl z11.s, p4/M, z11.s, z17.s\n"
    ".inst 0x04b37400  // sqrdmulh z0.s, z0.s, z19.s\n"
    ".inst 0x4482922a  // srshl z10.s, p4/M, z10.s, z17.s\n"
    ".inst 0x44829229  // srshl z9.s, p4/M, z9.s, z17.s\n"
    "add z15.s, z15.s, z16.s\n"
    ".inst 0x44829228  // srshl z8.s, p4/M, z8.s, z17.s\n"
    ".inst 0x44829227  // srshl z7.s, p4/M, z7.s, z17.s\n"
    "add z14.s, z14.s, z16.s\n"
    "add z13.s, z13.s, z16.s\n"
    ".inst 0x44829226  // srshl z6.s, p4/M, z6.s, z17.s\n"
    ".inst 0x44829225  // srshl z5.s, p4/M, z5.s, z17.s\n"
    "add z12.s, z12.s, z16.s\n"
    "add z11.s, z11.s, z16.s\n"
    ".inst 0x44829224  // srshl z4.s, p4/M, z4.s, z17.s\n"
    ".inst 0x44829223  // srshl z3.s, p4/M, z3.s, z17.s\n"
    "add z10.s, z10.s, z16.s\n"
    "add z9.s, z9.s, z16.s\n"
    ".inst 0x44829222  // srshl z2.s, p4/M, z2.s, z17.s\n"
    ".inst 0x44829221  // srshl z1.s, p4/M, z1.s, z17.s\n"
    "add z8.s, z8.s, z16.s\n"
    "add z7.s, z7.s, z16.s\n"
    ".inst 0x44829220  // srshl z0.s, p4/M, z0.s, z17.s\n"
    "add z6.s, z6.s, z16.s\n"
    "add z5.s, z5.s, z16.s\n"
    "smax z15.s, p4/M, z15.s, z18.s\n"
    "add z4.s, z4.s, z16.s\n"
    "add z3.s, z3.s, z16.s\n"
    "smax z14.s, p4/M, z14.s, z18.s\n"
    "smax z13.s, p4/M, z13.s, z18.s\n"
    "add z2.s, z2.s, z16.s\n"
    "add z1.s, z1.s, z16.s\n"
    "smax z12.s, p4/M, z12.s, z18.s\n"
    "smax z11.s, p4/M, z11.s, z18.s\n"
    "add z0.s, z0.s, z16.s\n"
    "smax z10.s, p4/M, z10.s, z18.s\n"
    "smax z9.s, p4/M, z9.s, z18.s\n"
    "smax z8.s, p4/M, z8.s, z18.s\n"
    "smax z7.s, p4/M, z7.s, z18.s\n"
    "smax z6.s, p4/M, z6.s, z18.s\n"
    "smax z5.s, p4/M, z5.s, z18.s\n"
    "smax z4.s, p4/M, z4.s, z18.s\n"
    "smax z3.s, p4/M, z3.s, z18.s\n"
    "smax z2.s, p4/M, z2.s, z18.s\n"
    "smax z1.s, p4/M, z1.s, z18.s\n"
    "smax z0.s, p4/M, z0.s, z18.s\n"
    "smin z15.s, p4/M, z15.s, z20.s\n"
    "smin z14.s, p4/M, z14.s, z20.s\n"
    "smin z13.s, p4/M, z13.s, z20.s\n"
    "smin z12.s, p4/M, z12.s, z20.s\n"
    "smin z11.s, p4/M, z11.s, z20.s\n"
    "smin z10.s, p4/M, z10.s, z20.s\n"
    "smin z9.s, p4/M, z9.s, z20.s\n"
    "smin z8.s, p4/M, z8.s, z20.s\n"
    "smin z7.s, p4/M, z7.s, z20.s\n"
    "trn1 z19.h, z15.h, z14.h\n"
    "smin z6.s, p4/M, z6.s, z20.s\n"
    "smin z5.s, p4/M, z5.s, z20.s\n"
    "trn1 z17.h, z13.h, z12.h\n"
    "smin z4.s, p4/M, z4.s, z20.s\n"
    "smin z3.s, p4/M, z3.s, z20.s\n"
    "trn1 z18.h, z11.h, z10.h\n"
    "smin z2.s, p4/M, z2.s, z20.s\n"
    "smin z1.s, p4/M, z1.s, z20.s\n"
    "trn1 z16.h, z9.h, z8.h\n"
    "smin z0.s, p4/M, z0.s, z20.s\n"
    "trn1 z21.h, z7.h, z6.h\n"
    "trn1 z20.b, z19.b, z17.b\n"
    "trn1 z17.h, z5.h, z4.h\n"
    "trn1 z19.h, z3.h, z2.h\n"
    "trn1 z18.b, z18.b, z16.b\n"
    "trn1 z16.h, z1.h, z0.h\n"
    "st1b { z20.b }, p3, [%x[outptr], x27]\n"
    "incb x27, ALL, MUL #4\n"
    "trn1 z17.b, z21.b, z17.b\n"
    "trn1 z16.b, z19.b, z16.b\n"
    "st1b { z18.b }, p2, [%x[outptr], x26]\n"
    "incb x26, ALL, MUL #4\n"
    "st1b { z17.b }, p1, [%x[outptr], x25]\n"
    "incb x25, ALL, MUL #4\n"
    "st1b { z16.b }, p0, [%x[outptr], x24]\n"
    "incb x24, ALL, MUL #4\n"
    "whilelt p0.b, x24, %x[n_channels]\n"
    "b.any 1b\n"
    "7:"  // Single vector of channels
    "whilelt p3.b, x27, %x[n_channels]\n"
    "b.none 14f\n"
    "8:"  // Single vector of channels: Loop
    "ld1rw { z15.s }, p4/Z, [%x[accumulator_init]]\n"
    "lsr x23, %x[n_valid_cells], #0x1\n"
    "mov x22, %x[inptrs]\n"
    "mov z14.d, z15.d\n"
    "mov z13.d, z15.d\n"
    "mov z12.d, z15.d\n"
    "cbz x23, 11f\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    "ld1b { z31.b }, p3/Z, [x21, x27]\n"
    "ld1b { z30.b }, p3/Z, [x20, x27]\n"
    "beq 10f\n"
    "9:"  // Single vector of channels: Loop: 2 inputs loop
    ".inst 0x455e0bf1  // uaddlb z17.h, z31.b, z30.b\n"
    ".inst 0x455e0ff0  // uaddlt z16.h, z31.b, z30.b\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    ".inst 0x459149ef  // uaddwb z15.s, z15.s, z17.h\n"
    ".inst 0x45914dce  // uaddwt z14.s, z14.s, z17.h\n"
    ".inst 0x459049ad  // uaddwb z13.s, z13.s, z16.h\n"
    ".inst 0x45904d8c  // uaddwt z12.s, z12.s, z16.h\n"
    "ld1b { z31.b }, p3/Z, [x21, x27]\n"
    "ld1b { z30.b }, p3/Z, [x20, x27]\n"
    "bgt 9b\n"
    "10:"  // Single vector of channels: Loop: 2 inputs tail
    ".inst 0x455e0bf1  // uaddlb z17.h, z31.b, z30.b\n"
    ".inst 0x455e0ff0  // uaddlt z16.h, z31.b, z30.b\n"
    ".inst 0x459149ef  // uaddwb z15.s, z15.s, z17.h\n"
    ".inst 0x45914dce  // uaddwt z14.s, z14.s, z17.h\n"
    ".inst 0x459049ad  // uaddwb z13.s, z13.s, z16.h\n"
    ".inst 0x45904d8c  // uaddwt z12.s, z12.s, z16.h\n"
    "11:"  // Single vector of channels: Loop: After loop
    "ands x21, %x[n_valid_cells], #0x1\n"
    "beq 13f\n"
    "12:"  // Single vector of channels: Loop: Single input loop
    "ldr x20, [x22], #0x8\n"
    "subs x21, x21, #0x1\n"
    "ld1b { z16.b }, p3/Z, [x20, x27]\n"
    ".inst 0x4508aa11  // ushllb z17.h, z16.b, #0x0\n"
    ".inst 0x4508ae10  // ushllt z16.h, z16.b, #0x0\n"
    ".inst 0x459149ef  // uaddwb z15.s, z15.s, z17.h\n"
    ".inst 0x45914dce  // uaddwt z14.s, z14.s, z17.h\n"
    ".inst 0x459049ad  // uaddwb z13.s, z13.s, z16.h\n"
    ".inst 0x45904d8c  // uaddwt z12.s, z12.s, z16.h\n"
    "bgt 12b\n"
    "13:"  // Single vector of channels: Loop: Single input loop: End
    "ld1rw { z21.s }, p4/Z, [%x[left_shift]]\n"
    "ld1rw { z20.s }, p4/Z, [%x[combined_rescale_value]]\n"
    "add x20, %x[quant_params], %[offsetof_qp_output_offset]\n"
    "mov z19.s, #0x0\n"
    "ld1rw { z18.s }, p4/Z, [%x[right_shift]]\n"
    "ld1rw { z17.s }, p4/Z, [x20]\n"
    "mov z16.s, #0xff\n"
    ".inst 0x448292af  // srshl z15.s, p4/M, z15.s, z21.s\n"
    ".inst 0x448292ae  // srshl z14.s, p4/M, z14.s, z21.s\n"
    ".inst 0x448292ad  // srshl z13.s, p4/M, z13.s, z21.s\n"
    ".inst 0x448292ac  // srshl z12.s, p4/M, z12.s, z21.s\n"
    ".inst 0x04b475ef  // sqrdmulh z15.s, z15.s, z20.s\n"
    ".inst 0x04b475ce  // sqrdmulh z14.s, z14.s, z20.s\n"
    ".inst 0x04b475ad  // sqrdmulh z13.s, z13.s, z20.s\n"
    ".inst 0x04b4758c  // sqrdmulh z12.s, z12.s, z20.s\n"
    ".inst 0x4482924f  // srshl z15.s, p4/M, z15.s, z18.s\n"
    ".inst 0x4482924e  // srshl z14.s, p4/M, z14.s, z18.s\n"
    ".inst 0x4482924d  // srshl z13.s, p4/M, z13.s, z18.s\n"
    ".inst 0x4482924c  // srshl z12.s, p4/M, z12.s, z18.s\n"
    "add z15.s, z15.s, z17.s\n"
    "add z14.s, z14.s, z17.s\n"
    "add z13.s, z13.s, z17.s\n"
    "add z12.s, z12.s, z17.s\n"
    "smax z15.s, p4/M, z15.s, z19.s\n"
    "smax z14.s, p4/M, z14.s, z19.s\n"
    "smax z13.s, p4/M, z13.s, z19.s\n"
    "smax z12.s, p4/M, z12.s, z19.s\n"
    "smin z15.s, p4/M, z15.s, z16.s\n"
    "smin z14.s, p4/M, z14.s, z16.s\n"
    "smin z13.s, p4/M, z13.s, z16.s\n"
    "smin z12.s, p4/M, z12.s, z16.s\n"
    "trn1 z17.h, z15.h, z14.h\n"
    "trn1 z16.h, z13.h, z12.h\n"
    "trn1 z16.b, z17.b, z16.b\n"
    "st1b { z16.b }, p3, [%x[outptr], x27]\n"
    "incb x27\n"
    "whilelt p3.b, x27, %x[n_channels]\n"
    "b.any 8b\n"
    "14:"  // End
    :
    : [accumulator_init] "r" (&accumulator_init), [combined_rescale_value] "r" (&combined_rescale_value), [inptrs] "r" (inptrs), [left_shift] "r" (&left_shift), [n_channels] "r" (n_channels), [n_valid_cells] "r" (n_valid_cells), [offsetof_qp_output_offset] "I" (offsetof(Requantize32, output_offset)), [outptr] "r" (outptr), [quant_params] "r" (&qp), [right_shift] "r" (&right_shift)
    : "cc", "memory", "p0", "p1", "p2", "p3", "p4", "x20", "x21", "x22", "x23", "x24", "x25", "x26", "x27", "z0", "z1", "z2", "z3", "z4", "z5", "z6", "z7", "z8", "z9", "z10", "z11", "z12", "z13", "z14", "z15", "z16", "z17", "z18", "z19", "z20", "z21", "z22", "z23", "z24", "z25", "z26", "z27", "z28", "z29", "z30", "z31"
  );
}

}  // namespace pooling
}  // namespace arm_conv

#endif  // defined(ARM_COMPUTE_ENABLE_SVE)
