/*
 * Copyright (c) 2022-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_SME)

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 sme_s8q_nhwc_avg_generic_depthfirst_impl(
  const uint64_t window_cells,
  const uint64_t n_valid_cells,
  uint64_t n_channels,
  const int8_t *const *const inptrs,
  int8_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);
  }

  // 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__(
    ".inst 0xd503477f  // SMSTART ZA\n"
    "mov x27, #0x0\n"
    "cntb x26\n"
    "cntb x25, ALL, MUL #2\n"
    "cntb x24, ALL, MUL #3\n"
    "whilelt p4.b, x27, %x[n_channels]\n"
    "whilelt p3.b, x26, %x[n_channels]\n"
    "whilelt p2.b, x25, %x[n_channels]\n"
    "whilelt p1.b, x24, %x[n_channels]\n"
    "ptrue p0.b\n"
    "b.none 7f\n"
    "1:"  // 4-vectors of channels
    "lsr x23, %x[n_valid_cells], #0x1\n"
    "mov z15.s, #0x0\n"
    "mov z14.s, #0x0\n"
    "mov x22, %x[inptrs]\n"
    "mov z13.s, #0x0\n"
    "mov z12.s, #0x0\n"
    "mov z11.s, #0x0\n"
    "mov z10.s, #0x0\n"
    "mov z9.s, #0x0\n"
    "mov z8.s, #0x0\n"
    "mov z7.s, #0x0\n"
    "mov z6.s, #0x0\n"
    "mov z5.s, #0x0\n"
    "mov z4.s, #0x0\n"
    "mov z3.s, #0x0\n"
    "mov z2.s, #0x0\n"
    "mov z1.s, #0x0\n"
    "mov z0.s, #0x0\n"
    "cbz x23, 4f\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    "ld1b { z31.b }, p4/Z, [x21, x27]\n"
    "ld1b { z30.b }, p4/Z, [x20, x27]\n"
    "ld1b { z29.b }, p3/Z, [x21, x26]\n"
    "ld1b { z28.b }, p3/Z, [x20, x26]\n"
    "ld1b { z27.b }, p2/Z, [x21, x25]\n"
    "ld1b { z26.b }, p2/Z, [x20, x25]\n"
    "ld1b { z25.b }, p1/Z, [x21, x24]\n"
    "ld1b { z24.b }, p1/Z, [x20, x24]\n"
    "beq 3f\n"
    "2:"  // 4-vectors of channels: 2 inputs loop
    ".inst 0x455e03f7  // saddlb z23.h, z31.b, z30.b\n"
    ".inst 0x455e07f6  // saddlt z22.h, z31.b, z30.b\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    ".inst 0x455c03b5  // saddlb z21.h, z29.b, z28.b\n"
    ".inst 0x455c07b4  // saddlt z20.h, z29.b, z28.b\n"
    "add x22, x22, #0x10\n"
    ".inst 0x455a0373  // saddlb z19.h, z27.b, z26.b\n"
    ".inst 0x455a0772  // saddlt z18.h, z27.b, z26.b\n"
    ".inst 0x45580331  // saddlb z17.h, z25.b, z24.b\n"
    ".inst 0x45580730  // saddlt z16.h, z25.b, z24.b\n"
    "ld1b { z31.b }, p4/Z, [x21, x27]\n"
    ".inst 0x459741ef  // saddwb z15.s, z15.s, z23.h\n"
    ".inst 0x459745ce  // saddwt z14.s, z14.s, z23.h\n"
    "ld1b { z30.b }, p4/Z, [x20, x27]\n"
    ".inst 0x459641ad  // saddwb z13.s, z13.s, z22.h\n"
    ".inst 0x4596458c  // saddwt z12.s, z12.s, z22.h\n"
    "ld1b { z29.b }, p3/Z, [x21, x26]\n"
    ".inst 0x4595416b  // saddwb z11.s, z11.s, z21.h\n"
    ".inst 0x4595454a  // saddwt z10.s, z10.s, z21.h\n"
    "ld1b { z28.b }, p3/Z, [x20, x26]\n"
    ".inst 0x45944129  // saddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944508  // saddwt z8.s, z8.s, z20.h\n"
    "ld1b { z27.b }, p2/Z, [x21, x25]\n"
    ".inst 0x459340e7  // saddwb z7.s, z7.s, z19.h\n"
    ".inst 0x459344c6  // saddwt z6.s, z6.s, z19.h\n"
    "ld1b { z26.b }, p2/Z, [x20, x25]\n"
    ".inst 0x459240a5  // saddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924484  // saddwt z4.s, z4.s, z18.h\n"
    "ld1b { z25.b }, p1/Z, [x21, x24]\n"
    ".inst 0x45914063  // saddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914442  // saddwt z2.s, z2.s, z17.h\n"
    "ld1b { z24.b }, p1/Z, [x20, x24]\n"
    ".inst 0x45904021  // saddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904400  // saddwt z0.s, z0.s, z16.h\n"
    "bgt 2b\n"
    "3:"  // 4-vectors of channels: 2 inputs tail
    ".inst 0x455e03f7  // saddlb z23.h, z31.b, z30.b\n"
    ".inst 0x455e07f6  // saddlt z22.h, z31.b, z30.b\n"
    ".inst 0x455c03b5  // saddlb z21.h, z29.b, z28.b\n"
    ".inst 0x455c07b4  // saddlt z20.h, z29.b, z28.b\n"
    ".inst 0x455a0373  // saddlb z19.h, z27.b, z26.b\n"
    ".inst 0x455a0772  // saddlt z18.h, z27.b, z26.b\n"
    ".inst 0x45580331  // saddlb z17.h, z25.b, z24.b\n"
    ".inst 0x45580730  // saddlt z16.h, z25.b, z24.b\n"
    ".inst 0x459741ef  // saddwb z15.s, z15.s, z23.h\n"
    ".inst 0x459745ce  // saddwt z14.s, z14.s, z23.h\n"
    ".inst 0x459641ad  // saddwb z13.s, z13.s, z22.h\n"
    ".inst 0x4596458c  // saddwt z12.s, z12.s, z22.h\n"
    ".inst 0x4595416b  // saddwb z11.s, z11.s, z21.h\n"
    ".inst 0x4595454a  // saddwt z10.s, z10.s, z21.h\n"
    ".inst 0x45944129  // saddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944508  // saddwt z8.s, z8.s, z20.h\n"
    ".inst 0x459340e7  // saddwb z7.s, z7.s, z19.h\n"
    ".inst 0x459344c6  // saddwt z6.s, z6.s, z19.h\n"
    ".inst 0x459240a5  // saddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924484  // saddwt z4.s, z4.s, z18.h\n"
    ".inst 0x45914063  // saddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914442  // saddwt z2.s, z2.s, z17.h\n"
    ".inst 0x45904021  // saddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904400  // saddwt 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 }, p4/Z, [x20, x27]\n"
    "ld1b { z18.b }, p3/Z, [x20, x26]\n"
    "ld1b { z17.b }, p2/Z, [x20, x25]\n"
    "ld1b { z16.b }, p1/Z, [x20, x24]\n"
    ".inst 0x4508a277  // sshllb z23.h, z19.b, #0x0\n"
    ".inst 0x4508a676  // sshllt z22.h, z19.b, #0x0\n"
    ".inst 0x4508a255  // sshllb z21.h, z18.b, #0x0\n"
    ".inst 0x4508a654  // sshllt z20.h, z18.b, #0x0\n"
    ".inst 0x4508a233  // sshllb z19.h, z17.b, #0x0\n"
    ".inst 0x4508a632  // sshllt z18.h, z17.b, #0x0\n"
    ".inst 0x4508a211  // sshllb z17.h, z16.b, #0x0\n"
    ".inst 0x4508a610  // sshllt z16.h, z16.b, #0x0\n"
    ".inst 0x459741ef  // saddwb z15.s, z15.s, z23.h\n"
    ".inst 0x459745ce  // saddwt z14.s, z14.s, z23.h\n"
    ".inst 0x459641ad  // saddwb z13.s, z13.s, z22.h\n"
    ".inst 0x4596458c  // saddwt z12.s, z12.s, z22.h\n"
    ".inst 0x4595416b  // saddwb z11.s, z11.s, z21.h\n"
    ".inst 0x4595454a  // saddwt z10.s, z10.s, z21.h\n"
    ".inst 0x45944129  // saddwb z9.s, z9.s, z20.h\n"
    ".inst 0x45944508  // saddwt z8.s, z8.s, z20.h\n"
    ".inst 0x459340e7  // saddwb z7.s, z7.s, z19.h\n"
    ".inst 0x459344c6  // saddwt z6.s, z6.s, z19.h\n"
    ".inst 0x459240a5  // saddwb z5.s, z5.s, z18.h\n"
    ".inst 0x45924484  // saddwt z4.s, z4.s, z18.h\n"
    ".inst 0x45914063  // saddwb z3.s, z3.s, z17.h\n"
    ".inst 0x45914442  // saddwt z2.s, z2.s, z17.h\n"
    ".inst 0x45904021  // saddwb z1.s, z1.s, z16.h\n"
    ".inst 0x45904400  // saddwt z0.s, z0.s, z16.h\n"
    "bgt 5b\n"
    "6:"  // 4-vectors of channels: Single input loop: End
    "ld1rw { z19.s }, p0/Z, [%x[left_shift]]\n"
    "mov z18.s, #0x7f\n"
    "ld1rw { z17.s }, p0/Z, [%x[combined_rescale_value]]\n"
    "ld1rw { z16.s }, p0/Z, [%x[right_shift]]\n"
    ".inst 0x4482826f  // srshl z15.s, p0/M, z15.s, z19.s\n"
    ".inst 0x4482826e  // srshl z14.s, p0/M, z14.s, z19.s\n"
    ".inst 0x4482826d  // srshl z13.s, p0/M, z13.s, z19.s\n"
    ".inst 0x4482826c  // srshl z12.s, p0/M, z12.s, z19.s\n"
    ".inst 0x4482826b  // srshl z11.s, p0/M, z11.s, z19.s\n"
    ".inst 0x4482826a  // srshl z10.s, p0/M, z10.s, z19.s\n"
    ".inst 0x44828269  // srshl z9.s, p0/M, z9.s, z19.s\n"
    ".inst 0x44828268  // srshl z8.s, p0/M, z8.s, z19.s\n"
    ".inst 0x44828267  // srshl z7.s, p0/M, z7.s, z19.s\n"
    ".inst 0x44828266  // srshl z6.s, p0/M, z6.s, z19.s\n"
    ".inst 0x44828265  // srshl z5.s, p0/M, z5.s, z19.s\n"
    ".inst 0x44828264  // srshl z4.s, p0/M, z4.s, z19.s\n"
    ".inst 0x44828263  // srshl z3.s, p0/M, z3.s, z19.s\n"
    ".inst 0x44828262  // srshl z2.s, p0/M, z2.s, z19.s\n"
    ".inst 0x44828261  // srshl z1.s, p0/M, z1.s, z19.s\n"
    ".inst 0x44828260  // srshl z0.s, p0/M, z0.s, z19.s\n"
    ".inst 0x04b175ef  // sqrdmulh z15.s, z15.s, z17.s\n"
    ".inst 0x04b175ce  // sqrdmulh z14.s, z14.s, z17.s\n"
    ".inst 0x04b175ad  // sqrdmulh z13.s, z13.s, z17.s\n"
    ".inst 0x04b1758c  // sqrdmulh z12.s, z12.s, z17.s\n"
    ".inst 0x04b1756b  // sqrdmulh z11.s, z11.s, z17.s\n"
    ".inst 0x04b1754a  // sqrdmulh z10.s, z10.s, z17.s\n"
    ".inst 0x04b17529  // sqrdmulh z9.s, z9.s, z17.s\n"
    ".inst 0x04b17508  // sqrdmulh z8.s, z8.s, z17.s\n"
    ".inst 0x04b174e7  // sqrdmulh z7.s, z7.s, z17.s\n"
    ".inst 0x04b174c6  // sqrdmulh z6.s, z6.s, z17.s\n"
    ".inst 0x04b174a5  // sqrdmulh z5.s, z5.s, z17.s\n"
    ".inst 0x04b17484  // sqrdmulh z4.s, z4.s, z17.s\n"
    ".inst 0x04b17463  // sqrdmulh z3.s, z3.s, z17.s\n"
    ".inst 0x04b17442  // sqrdmulh z2.s, z2.s, z17.s\n"
    ".inst 0x04b17421  // sqrdmulh z1.s, z1.s, z17.s\n"
    ".inst 0x04b17400  // sqrdmulh z0.s, z0.s, z17.s\n"
    ".inst 0x4482820f  // srshl z15.s, p0/M, z15.s, z16.s\n"
    ".inst 0x4482820e  // srshl z14.s, p0/M, z14.s, z16.s\n"
    ".inst 0x4482820d  // srshl z13.s, p0/M, z13.s, z16.s\n"
    ".inst 0x4482820c  // srshl z12.s, p0/M, z12.s, z16.s\n"
    ".inst 0x4482820b  // srshl z11.s, p0/M, z11.s, z16.s\n"
    ".inst 0x4482820a  // srshl z10.s, p0/M, z10.s, z16.s\n"
    ".inst 0x44828209  // srshl z9.s, p0/M, z9.s, z16.s\n"
    ".inst 0x44828208  // srshl z8.s, p0/M, z8.s, z16.s\n"
    ".inst 0x44828207  // srshl z7.s, p0/M, z7.s, z16.s\n"
    ".inst 0x44828206  // srshl z6.s, p0/M, z6.s, z16.s\n"
    ".inst 0x44828205  // srshl z5.s, p0/M, z5.s, z16.s\n"
    ".inst 0x44828204  // srshl z4.s, p0/M, z4.s, z16.s\n"
    ".inst 0x44828203  // srshl z3.s, p0/M, z3.s, z16.s\n"
    ".inst 0x44828202  // srshl z2.s, p0/M, z2.s, z16.s\n"
    ".inst 0x44828201  // srshl z1.s, p0/M, z1.s, z16.s\n"
    ".inst 0x44828200  // srshl z0.s, p0/M, z0.s, z16.s\n"
    "not z16.s, p0/M, z18.s\n"
    "smax z15.s, p0/M, z15.s, z16.s\n"
    "smax z14.s, p0/M, z14.s, z16.s\n"
    "smax z13.s, p0/M, z13.s, z16.s\n"
    "smax z12.s, p0/M, z12.s, z16.s\n"
    "smax z11.s, p0/M, z11.s, z16.s\n"
    "smax z10.s, p0/M, z10.s, z16.s\n"
    "smax z9.s, p0/M, z9.s, z16.s\n"
    "smax z8.s, p0/M, z8.s, z16.s\n"
    "smax z7.s, p0/M, z7.s, z16.s\n"
    "smax z6.s, p0/M, z6.s, z16.s\n"
    "smax z5.s, p0/M, z5.s, z16.s\n"
    "smax z4.s, p0/M, z4.s, z16.s\n"
    "smax z3.s, p0/M, z3.s, z16.s\n"
    "smax z2.s, p0/M, z2.s, z16.s\n"
    "smax z1.s, p0/M, z1.s, z16.s\n"
    "smax z0.s, p0/M, z0.s, z16.s\n"
    "smin z15.s, p0/M, z15.s, z18.s\n"
    "smin z14.s, p0/M, z14.s, z18.s\n"
    "smin z13.s, p0/M, z13.s, z18.s\n"
    "smin z12.s, p0/M, z12.s, z18.s\n"
    "smin z11.s, p0/M, z11.s, z18.s\n"
    "smin z10.s, p0/M, z10.s, z18.s\n"
    "smin z9.s, p0/M, z9.s, z18.s\n"
    "smin z8.s, p0/M, z8.s, z18.s\n"
    "smin z7.s, p0/M, z7.s, z18.s\n"
    "smin z6.s, p0/M, z6.s, z18.s\n"
    "trn1 z19.h, z15.h, z14.h\n"
    "smin z5.s, p0/M, z5.s, z18.s\n"
    "smin z4.s, p0/M, z4.s, z18.s\n"
    "trn1 z16.h, z13.h, z12.h\n"
    "smin z3.s, p0/M, z3.s, z18.s\n"
    "smin z2.s, p0/M, z2.s, z18.s\n"
    "trn1 z22.h, z11.h, z10.h\n"
    "smin z1.s, p0/M, z1.s, z18.s\n"
    "smin z0.s, p0/M, z0.s, z18.s\n"
    "trn1 z18.h, z9.h, z8.h\n"
    "trn1 z21.h, z7.h, z6.h\n"
    "trn1 z17.h, z5.h, z4.h\n"
    "trn1 z20.b, z19.b, z16.b\n"
    "trn1 z19.h, z3.h, z2.h\n"
    "trn1 z16.h, z1.h, z0.h\n"
    "trn1 z18.b, z22.b, z18.b\n"
    "trn1 z17.b, z21.b, z17.b\n"
    "st1b { z20.b }, p4, [%x[outptr], x27]\n"
    "incb x27, ALL, MUL #4\n"
    "trn1 z16.b, z19.b, z16.b\n"
    "st1b { z18.b }, p3, [%x[outptr], x26]\n"
    "incb x26, ALL, MUL #4\n"
    "st1b { z17.b }, p2, [%x[outptr], x25]\n"
    "incb x25, ALL, MUL #4\n"
    "st1b { z16.b }, p1, [%x[outptr], x24]\n"
    "incb x24, ALL, MUL #4\n"
    "whilelt p1.b, x24, %x[n_channels]\n"
    "b.any 1b\n"
    "7:"  // Single vector of channels
    "whilelt p4.b, x27, %x[n_channels]\n"
    "b.none 14f\n"
    "8:"  // Single vector of channels: Loop
    "lsr x23, %x[n_valid_cells], #0x1\n"
    "mov z15.s, #0x0\n"
    "mov z14.s, #0x0\n"
    "mov x22, %x[inptrs]\n"
    "mov z13.s, #0x0\n"
    "mov z12.s, #0x0\n"
    "cbz x23, 11f\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    "ld1b { z31.b }, p4/Z, [x21, x27]\n"
    "ld1b { z30.b }, p4/Z, [x20, x27]\n"
    "beq 10f\n"
    "9:"  // Single vector of channels: Loop: 2 inputs loop
    ".inst 0x455e03f1  // saddlb z17.h, z31.b, z30.b\n"
    ".inst 0x455e07f0  // saddlt z16.h, z31.b, z30.b\n"
    "ldp x21, x20, [x22, #0x0]\n"
    "subs x23, x23, #0x1\n"
    "add x22, x22, #0x10\n"
    "ld1b { z31.b }, p4/Z, [x21, x27]\n"
    ".inst 0x459141ef  // saddwb z15.s, z15.s, z17.h\n"
    ".inst 0x459145ce  // saddwt z14.s, z14.s, z17.h\n"
    "ld1b { z30.b }, p4/Z, [x20, x27]\n"
    ".inst 0x459041ad  // saddwb z13.s, z13.s, z16.h\n"
    ".inst 0x4590458c  // saddwt z12.s, z12.s, z16.h\n"
    "bgt 9b\n"
    "10:"  // Single vector of channels: Loop: 2 inputs tail
    ".inst 0x455e03f1  // saddlb z17.h, z31.b, z30.b\n"
    ".inst 0x455e07f0  // saddlt z16.h, z31.b, z30.b\n"
    ".inst 0x459141ef  // saddwb z15.s, z15.s, z17.h\n"
    ".inst 0x459145ce  // saddwt z14.s, z14.s, z17.h\n"
    ".inst 0x459041ad  // saddwb z13.s, z13.s, z16.h\n"
    ".inst 0x4590458c  // saddwt 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 }, p4/Z, [x20, x27]\n"
    ".inst 0x4508a211  // sshllb z17.h, z16.b, #0x0\n"
    ".inst 0x4508a610  // sshllt z16.h, z16.b, #0x0\n"
    ".inst 0x459141ef  // saddwb z15.s, z15.s, z17.h\n"
    ".inst 0x459145ce  // saddwt z14.s, z14.s, z17.h\n"
    ".inst 0x459041ad  // saddwb z13.s, z13.s, z16.h\n"
    ".inst 0x4590458c  // saddwt z12.s, z12.s, z16.h\n"
    "bgt 12b\n"
    "13:"  // Single vector of channels: Loop: Single input loop: End
    "ld1rw { z16.s }, p0/Z, [%x[left_shift]]\n"
    "mov z19.s, #0x7f\n"
    "ld1rw { z18.s }, p0/Z, [%x[combined_rescale_value]]\n"
    "ld1rw { z17.s }, p0/Z, [%x[right_shift]]\n"
    ".inst 0x4482820f  // srshl z15.s, p0/M, z15.s, z16.s\n"
    ".inst 0x4482820e  // srshl z14.s, p0/M, z14.s, z16.s\n"
    ".inst 0x4482820d  // srshl z13.s, p0/M, z13.s, z16.s\n"
    ".inst 0x4482820c  // srshl z12.s, p0/M, z12.s, z16.s\n"
    "not z16.s, p0/M, z19.s\n"
    ".inst 0x04b275ef  // sqrdmulh z15.s, z15.s, z18.s\n"
    ".inst 0x04b275ce  // sqrdmulh z14.s, z14.s, z18.s\n"
    ".inst 0x04b275ad  // sqrdmulh z13.s, z13.s, z18.s\n"
    ".inst 0x04b2758c  // sqrdmulh z12.s, z12.s, z18.s\n"
    ".inst 0x4482822f  // srshl z15.s, p0/M, z15.s, z17.s\n"
    ".inst 0x4482822e  // srshl z14.s, p0/M, z14.s, z17.s\n"
    ".inst 0x4482822d  // srshl z13.s, p0/M, z13.s, z17.s\n"
    ".inst 0x4482822c  // srshl z12.s, p0/M, z12.s, z17.s\n"
    "smax z15.s, p0/M, z15.s, z16.s\n"
    "smax z14.s, p0/M, z14.s, z16.s\n"
    "smax z13.s, p0/M, z13.s, z16.s\n"
    "smax z12.s, p0/M, z12.s, z16.s\n"
    "smin z15.s, p0/M, z15.s, z19.s\n"
    "smin z14.s, p0/M, z14.s, z19.s\n"
    "smin z13.s, p0/M, z13.s, z19.s\n"
    "smin z12.s, p0/M, z12.s, z19.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 }, p4, [%x[outptr], x27]\n"
    "incb x27\n"
    "whilelt p4.b, x27, %x[n_channels]\n"
    "b.any 8b\n"
    "14:"  // End
    ".inst 0xd503467f  // SMSTOP\n"
    :
    : [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), [outptr] "r" (outptr), [right_shift] "r" (&right_shift)
    : "cc", "memory", "p0", "p1", "p2", "p3", "p4", "p5", "p6", "p7", "p8", "p9", "p10", "p11", "p12", "p13", "p14", "p15", "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_SME)
