/* * Copyright (c) 2021 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 #include #include #if defined(__ARM_FEATURE_SVE) && defined(SVE2) 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_u8_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 ) { 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 ¶ms = rescale_params[window_cells - 2]; rescale_value = params.multiplier; shift_value = params.shift; } else { auto f_rescale_value = 1.0f / static_cast(window_cells); shift_value = 0; while (f_rescale_value < 0.5f) { shift_value--; f_rescale_value *= 2.0f; } rescale_value = static_cast(round(f_rescale_value * static_cast(1ll << 31))); if (static_cast(rescale_value) == (1ll << 31)) { shift_value++; rescale_value >>= 1; } } __asm__ __volatile__( "ptrue p4.b\n" "mov x26, #0x0\n" "cntb x25\n" "cntb x24, ALL, MUL #2\n" "cntb x23, ALL, MUL #3\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" "whilelt p0.b, x23, %x[n_channels]\n" "b.none 7f\n" "1:" // 4-vectors of channels "mov z15.s, #0x0\n" "mov x19, %x[inptrs]\n" "mov z14.s, #0x0\n" "lsr x22, %x[n_valid_cells], #0x1\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 x22, 4f\n" "ldp x21, x20, [x19, #0x0]\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" "add x19, x19, #0x10\n" "ld1b { z30.b }, p3/Z, [x20, x26]\n" "subs x22, x22, #0x1\n" "ld1b { z29.b }, p2/Z, [x21, x25]\n" "ld1b { z28.b }, p2/Z, [x20, x25]\n" "ld1b { z27.b }, p1/Z, [x21, x24]\n" "ld1b { z26.b }, p1/Z, [x20, x24]\n" "ld1b { z25.b }, p0/Z, [x21, x23]\n" "ld1b { z24.b }, p0/Z, [x20, x23]\n" "beq 3f\n" "2:" // 4-vectors of channels: 2 inputs loop ".inst 0x455e0bf7 // uaddlb z23.h, z31.b, z30.b\n" "ldp x21, x20, [x19, #0x0]\n" "add x19, x19, #0x10\n" ".inst 0x455e0ff6 // uaddlt z22.h, z31.b, z30.b\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" ".inst 0x455c0bb5 // uaddlb z21.h, z29.b, z28.b\n" "subs x22, x22, #0x1\n" ".inst 0x455c0fb4 // uaddlt z20.h, z29.b, z28.b\n" "ld1b { z30.b }, p3/Z, [x20, x26]\n" ".inst 0x455a0b73 // uaddlb z19.h, z27.b, z26.b\n" "ld1b { z29.b }, p2/Z, [x21, x25]\n" ".inst 0x455a0f72 // uaddlt z18.h, z27.b, z26.b\n" "ld1b { z28.b }, p2/Z, [x20, x25]\n" ".inst 0x45580b31 // uaddlb z17.h, z25.b, z24.b\n" "ld1b { z27.b }, p1/Z, [x21, x24]\n" ".inst 0x45580f30 // uaddlt z16.h, z25.b, z24.b\n" "ld1b { z26.b }, p1/Z, [x20, x24]\n" ".inst 0x459749ef // uaddwb z15.s, z15.s, z23.h\n" "ld1b { z25.b }, p0/Z, [x21, x23]\n" ".inst 0x45974dce // uaddwt z14.s, z14.s, z23.h\n" "ld1b { z24.b }, p0/Z, [x20, x23]\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 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 x20, %x[n_valid_cells], #0x1\n" "beq 6f\n" "5:" // 4-vectors of channels: Single input loop "ldr x21, [x19], #0x8\n" "subs x20, x20, #0x1\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" ".inst 0x4508abf1 // ushllb z17.h, z31.b, #0x0\n" "ld1b { z29.b }, p2/Z, [x21, x25]\n" ".inst 0x4508aff0 // ushllt z16.h, z31.b, #0x0\n" "ld1b { z27.b }, p1/Z, [x21, x24]\n" ".inst 0x459149ef // uaddwb z15.s, z15.s, z17.h\n" "ld1b { z25.b }, p0/Z, [x21, x23]\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" ".inst 0x4508abb0 // ushllb z16.h, z29.b, #0x0\n" ".inst 0x4590496b // uaddwb z11.s, z11.s, z16.h\n" ".inst 0x45904d4a // uaddwt z10.s, z10.s, z16.h\n" ".inst 0x4508afb0 // ushllt z16.h, z29.b, #0x0\n" ".inst 0x45904929 // uaddwb z9.s, z9.s, z16.h\n" ".inst 0x45904d08 // uaddwt z8.s, z8.s, z16.h\n" ".inst 0x4508ab70 // ushllb z16.h, z27.b, #0x0\n" ".inst 0x459048e7 // uaddwb z7.s, z7.s, z16.h\n" ".inst 0x45904cc6 // uaddwt z6.s, z6.s, z16.h\n" ".inst 0x4508af70 // ushllt z16.h, z27.b, #0x0\n" ".inst 0x459048a5 // uaddwb z5.s, z5.s, z16.h\n" ".inst 0x45904c84 // uaddwt z4.s, z4.s, z16.h\n" ".inst 0x4508ab30 // ushllb z16.h, z25.b, #0x0\n" ".inst 0x45904863 // uaddwb z3.s, z3.s, z16.h\n" ".inst 0x45904c42 // uaddwt z2.s, z2.s, z16.h\n" ".inst 0x4508af30 // ushllt z16.h, z25.b, #0x0\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 "mov z20.s, #0x0\n" "ld1rw { z17.s }, p4/Z, [%x[rescale_ptr]]\n" "mov z19.s, #0xff\n" "ld1rw { z16.s }, p4/Z, [%x[shift_ptr]]\n" ".inst 0x04b175ef // sqdmulh z15.s, z15.s, z17.s\n" ".inst 0x04b175ce // sqdmulh z14.s, z14.s, z17.s\n" ".inst 0x04b175ad // sqdmulh z13.s, z13.s, z17.s\n" ".inst 0x04b1758c // sqdmulh z12.s, z12.s, z17.s\n" ".inst 0x04b1756b // sqdmulh z11.s, z11.s, z17.s\n" ".inst 0x04b1754a // sqdmulh z10.s, z10.s, z17.s\n" ".inst 0x04b17529 // sqdmulh z9.s, z9.s, z17.s\n" ".inst 0x04b17508 // sqdmulh z8.s, z8.s, z17.s\n" ".inst 0x04b174e7 // sqdmulh z7.s, z7.s, z17.s\n" ".inst 0x04b174c6 // sqdmulh z6.s, z6.s, z17.s\n" ".inst 0x04b174a5 // sqdmulh z5.s, z5.s, z17.s\n" ".inst 0x04b17484 // sqdmulh z4.s, z4.s, z17.s\n" ".inst 0x04b17463 // sqdmulh z3.s, z3.s, z17.s\n" ".inst 0x04b17442 // sqdmulh z2.s, z2.s, z17.s\n" ".inst 0x04b17421 // sqdmulh z1.s, z1.s, z17.s\n" ".inst 0x04b17400 // sqdmulh z0.s, z0.s, z17.s\n" ".inst 0x4482920f // srshl z15.s, p4/M, z15.s, z16.s\n" ".inst 0x4482920e // srshl z14.s, p4/M, z14.s, z16.s\n" ".inst 0x4482920d // srshl z13.s, p4/M, z13.s, z16.s\n" ".inst 0x4482920c // srshl z12.s, p4/M, z12.s, z16.s\n" ".inst 0x4482920b // srshl z11.s, p4/M, z11.s, z16.s\n" ".inst 0x4482920a // srshl z10.s, p4/M, z10.s, z16.s\n" ".inst 0x44829209 // srshl z9.s, p4/M, z9.s, z16.s\n" ".inst 0x44829208 // srshl z8.s, p4/M, z8.s, z16.s\n" ".inst 0x44829207 // srshl z7.s, p4/M, z7.s, z16.s\n" ".inst 0x44829206 // srshl z6.s, p4/M, z6.s, z16.s\n" ".inst 0x44829205 // srshl z5.s, p4/M, z5.s, z16.s\n" ".inst 0x44829204 // srshl z4.s, p4/M, z4.s, z16.s\n" ".inst 0x44829203 // srshl z3.s, p4/M, z3.s, z16.s\n" ".inst 0x44829202 // srshl z2.s, p4/M, z2.s, z16.s\n" ".inst 0x44829201 // srshl z1.s, p4/M, z1.s, z16.s\n" ".inst 0x44829200 // srshl z0.s, p4/M, z0.s, z16.s\n" "smax z15.s, p4/M, z15.s, z20.s\n" "smax z14.s, p4/M, z14.s, z20.s\n" "smax z13.s, p4/M, z13.s, z20.s\n" "smax z12.s, p4/M, z12.s, z20.s\n" "smin z15.s, p4/M, z15.s, z19.s\n" "smin z14.s, p4/M, z14.s, z19.s\n" "smin z13.s, p4/M, z13.s, z19.s\n" "smin z12.s, p4/M, z12.s, z19.s\n" "smax z11.s, p4/M, z11.s, z20.s\n" "trn1 z17.h, z15.h, z14.h\n" "smax z10.s, p4/M, z10.s, z20.s\n" "trn1 z16.h, z13.h, z12.h\n" "smin z11.s, p4/M, z11.s, z19.s\n" "trn1 z16.b, z17.b, z16.b\n" "st1b { z16.b }, p3, [%x[outptr], x26]\n" "smin z10.s, p4/M, z10.s, z19.s\n" "incb x26, ALL, MUL #4\n" "smax z9.s, p4/M, z9.s, z20.s\n" "smax z8.s, p4/M, z8.s, z20.s\n" "smax z7.s, p4/M, z7.s, z20.s\n" "smax z6.s, p4/M, z6.s, z20.s\n" "trn1 z18.h, z11.h, z10.h\n" "smin z9.s, p4/M, z9.s, z19.s\n" "smin z8.s, p4/M, z8.s, z19.s\n" "smin z7.s, p4/M, z7.s, z19.s\n" "smin z6.s, p4/M, z6.s, z19.s\n" "smax z5.s, p4/M, z5.s, z20.s\n" "trn1 z16.h, z9.h, z8.h\n" "smax z4.s, p4/M, z4.s, z20.s\n" "trn1 z17.h, z7.h, z6.h\n" "trn1 z16.b, z18.b, z16.b\n" "st1b { z16.b }, p2, [%x[outptr], x25]\n" "smin z5.s, p4/M, z5.s, z19.s\n" "incb x25, ALL, MUL #4\n" "smin z4.s, p4/M, z4.s, z19.s\n" "smax z3.s, p4/M, z3.s, z20.s\n" "smax z2.s, p4/M, z2.s, z20.s\n" "smax z1.s, p4/M, z1.s, z20.s\n" "smax z0.s, p4/M, z0.s, z20.s\n" "trn1 z16.h, z5.h, z4.h\n" "smin z3.s, p4/M, z3.s, z19.s\n" "trn1 z16.b, z17.b, z16.b\n" "st1b { z16.b }, p1, [%x[outptr], x24]\n" "smin z2.s, p4/M, z2.s, z19.s\n" "incb x24, ALL, MUL #4\n" "smin z1.s, p4/M, z1.s, z19.s\n" "smin z0.s, p4/M, z0.s, z19.s\n" "trn1 z17.h, z3.h, z2.h\n" "trn1 z16.h, z1.h, z0.h\n" "trn1 z16.b, z17.b, z16.b\n" "st1b { z16.b }, p0, [%x[outptr], x23]\n" "incb x23, ALL, MUL #4\n" "whilelt p0.b, x23, %x[n_channels]\n" "b.any 1b\n" "7:" // Single vector of channels "whilelt p3.b, x26, %x[n_channels]\n" "b.none 14f\n" "8:" // Single vector of channels: Loop "mov z15.s, #0x0\n" "mov x19, %x[inptrs]\n" "mov z14.s, #0x0\n" "lsr x22, %x[n_valid_cells], #0x1\n" "mov z13.s, #0x0\n" "mov z12.s, #0x0\n" "cbz x22, 11f\n" "ldp x21, x20, [x19, #0x0]\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" "add x19, x19, #0x10\n" "ld1b { z30.b }, p3/Z, [x20, x26]\n" "subs x22, x22, #0x1\n" "beq 10f\n" "9:" // Single vector of channels: Loop: 2 inputs loop ".inst 0x455e0bf7 // uaddlb z23.h, z31.b, z30.b\n" "ldp x21, x20, [x19, #0x0]\n" "add x19, x19, #0x10\n" ".inst 0x455e0ff6 // uaddlt z22.h, z31.b, z30.b\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" ".inst 0x459749ef // uaddwb z15.s, z15.s, z23.h\n" "subs x22, x22, #0x1\n" ".inst 0x45974dce // uaddwt z14.s, z14.s, z23.h\n" "ld1b { z30.b }, p3/Z, [x20, x26]\n" ".inst 0x459649ad // uaddwb z13.s, z13.s, z22.h\n" ".inst 0x45964d8c // uaddwt z12.s, z12.s, z22.h\n" "bgt 9b\n" "10:" // Single vector of channels: Loop: 2 inputs tail ".inst 0x455e0bf7 // uaddlb z23.h, z31.b, z30.b\n" ".inst 0x455e0ff6 // uaddlt z22.h, z31.b, z30.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" "11:" // Single vector of channels: Loop: After loop "ands x20, %x[n_valid_cells], #0x1\n" "beq 13f\n" "12:" // Single vector of channels: Loop: Single input loop "ldr x21, [x19], #0x8\n" "subs x20, x20, #0x1\n" "ld1b { z31.b }, p3/Z, [x21, x26]\n" ".inst 0x4508abf1 // ushllb z17.h, z31.b, #0x0\n" ".inst 0x4508aff0 // ushllt z16.h, z31.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 "mov z20.s, #0x0\n" "ld1rw { z17.s }, p4/Z, [%x[rescale_ptr]]\n" "mov z19.s, #0xff\n" "ld1rw { z16.s }, p4/Z, [%x[shift_ptr]]\n" ".inst 0x04b175ef // sqdmulh z15.s, z15.s, z17.s\n" ".inst 0x04b175ce // sqdmulh z14.s, z14.s, z17.s\n" ".inst 0x04b175ad // sqdmulh z13.s, z13.s, z17.s\n" ".inst 0x04b1758c // sqdmulh z12.s, z12.s, z17.s\n" ".inst 0x4482920f // srshl z15.s, p4/M, z15.s, z16.s\n" ".inst 0x4482920e // srshl z14.s, p4/M, z14.s, z16.s\n" ".inst 0x4482920d // srshl z13.s, p4/M, z13.s, z16.s\n" ".inst 0x4482920c // srshl z12.s, p4/M, z12.s, z16.s\n" "smax z15.s, p4/M, z15.s, z20.s\n" "smax z14.s, p4/M, z14.s, z20.s\n" "smax z13.s, p4/M, z13.s, z20.s\n" "smax z12.s, p4/M, z12.s, z20.s\n" "smin z15.s, p4/M, z15.s, z19.s\n" "smin z14.s, p4/M, z14.s, z19.s\n" "smin z13.s, p4/M, z13.s, z19.s\n" "smin z12.s, p4/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 }, p3, [%x[outptr], x26]\n" "incb x26\n" "whilelt p3.b, x26, %x[n_channels]\n" "b.any 8b\n" "14:" // End : : [inptrs] "r" (inptrs), [n_channels] "r" (n_channels), [n_valid_cells] "r" (n_valid_cells), [outptr] "r" (outptr), [rescale_ptr] "r" (&rescale_value), [shift_ptr] "r" (&shift_value) : "cc", "memory", "p0", "p1", "p2", "p3", "p4", "x19", "x20", "x21", "x22", "x23", "x24", "x25", "x26", "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_FEATURE_SVE) && defined(SVE2)