/* * Copyright (c) 2022 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. */ #ifdef __ARM_FEATURE_SVE #ifdef ARM_COMPUTE_ENABLE_SME2 #include "arm_gemm.hpp" #include #include "../../asmlib.hpp" #include "../../utils.hpp" namespace arm_gemm { void sme2_interleaved_nomerge_s8s32_mopa_2VLx2VL(const int8_t *const A, const int8_t *const B, int32_t *const C, int ldc, const int M, const int N, const int K, const int32_t *const bias, const Activation act, bool accumulate, int32_t *const accumulator_buffer) { ARM_COMPUTE_UNUSED(act); struct KernelArgs { KernelArgs( const int8_t *const A, const int8_t *const B, int32_t *const C, const int ldc, const int M, const int N, const int K, const int32_t *const bias, bool accumulate, int32_t *const accumulator_buffer ) : A(A), B(B), kstride_bytes(roundup(K, 4) * sizeof(int8_t)), C(C), ldcb(ldc * sizeof(int32_t)), M(M), N(N), K(K), n_loops(((K / 4) - 1) / 2), n_tail_iters(((K / 4) - 1) % 2), bias(bias), accumulator_buffer(accumulator_buffer), flags(0x0) { if (accumulate) { flags |= 1 << 0; // FILL_ACCUMULATORS_FROM_BUFFER } if (C == nullptr) { flags |= 1 << 1; // STORE_ACCUMULATORS_TO_BUFFER } } const int8_t *const A; const int8_t *const B; const long kstride_bytes; int32_t *const C; const long ldcb; const long M, N, K, n_loops, n_tail_iters; const int32_t *const bias; int32_t *const accumulator_buffer; uint64_t flags; }; // Construct arguments for this kernel KernelArgs args(A, B, C, ldc, M, N, K, bias, accumulate, accumulator_buffer); __asm__ __volatile__( "ldr x15, [%x[args], %[offsetof_flags]]\n" ".inst 0xd503477f // SMSTART ZA\n" "ptrue p0.b\n" ".inst 0x25207811 // ptrue pn9.b\n" "ldr x14, [%x[args], %[offsetof_accumulator_buffer]]\n" "ldr x13, [%x[args], %[offsetof_accumulator_buffer]]\n" "tbz x15, #0, 2f\n" "mov x12, #0x0\n" "cntw x19\n" "1:" // Initial accumulator load from buffer: Loop ".inst 0xa040c5c8 // ld1w { z8.s-z11.s }, pn9.b/Z, [x14]\n" ".inst 0xc0840500 // mova za0h.s[x12], { z8.s-z11.s }\n" ".inst 0xa041c5d8 // ld1w { z24.s-z27.s }, pn9.b/Z, [x14, #0x4, MUL VL]\n" ".inst 0xc0840701 // mova za1h.s[x12], { z24.s-z27.s }\n" ".inst 0xa042c5dc // ld1w { z28.s-z31.s }, pn9.b/Z, [x14, #0x8, MUL VL]\n" ".inst 0xc0840782 // mova za2h.s[x12], { z28.s-z31.s }\n" ".inst 0xa043c5d0 // ld1w { z16.s-z19.s }, pn9.b/Z, [x14, #0xc, MUL VL]\n" ".inst 0xc0840603 // mova za3h.s[x12], { z16.s-z19.s }\n" "add x12, x12, #0x4\n" "cmp x12, x19\n" "addvl x14, x14, #16\n" "blt 1b\n" "2:" // Initial accumulator load from buffer: End "ldr w11, [%x[args], %[offsetof_M]]\n" "mov x10, #0x0\n" "mov x9, #0x0\n" "ldr w28, [%x[args], %[offsetof_N]]\n" "ldr x27, [%x[args], %[offsetof_A]]\n" "3:" // M and N loop "mov x26, x27\n" ".inst 0x25bc4530 // whilelt pn8.s, x9, x28, VLx2\n" "tbnz x15, #0, 4f\n" "ldr x19, [%x[args], %[offsetof_bias]]\n" ".inst 0xc00800ff // zero { zad0, zad1, zad2, zad3, zad4, zad5, zad6, zad7 }\n" "cbz x19, 5f\n" ".inst 0xa109427c // ldnt1w { z20.s, z28.s }, p8/Z, [x19, x9, LSL #2]\n" ".inst 0xc0900280 // addha za0.s, p0/M, p0/M, z20.s\n" ".inst 0xc0900381 // addha za1.s, p0/M, p0/M, z28.s\n" ".inst 0xc0900282 // addha za2.s, p0/M, p0/M, z20.s\n" ".inst 0xc0900383 // addha za3.s, p0/M, p0/M, z28.s\n" "4:" // Prepare accumulators: Test for last block "mov x19, x9\n" "mov x20, x10\n" "incw x19, ALL, MUL #2\n" "incw x20, ALL, MUL #2\n" "cmp x19, x28\n" "csel x20, x10, x20, LT\n" "mov x19, x15\n" "bfm x15, XZR, #0x0, #0x0 // bfc x15, #0x0, #0x1\n" "cmp x20, x11\n" "csel x15, x19, x15, LT\n" "5:" // Prepare accumulators: End "ldr x19, [%x[args], %[offsetof_K]]\n" "add x19, x19, #0x3\n" "lsr x19, x19, #0x2\n" "ldr x22, [%x[args], %[offsetof_B]]\n" "lsr x21, x19, #0x2\n" "and x20, x19, #0x3\n" "ldr x19, [%x[args], %[offsetof_kstride_bytes]]\n" "madd x22, x9, x19, x22\n" // bptr = B + n * kstride_bytes "cbz x21, 8f\n" "subs x21, x21, #0x1\n" ".inst 0xa1400756 // ld1b { z22.b, z30.b }, pn9.b/Z, [x26]\n" ".inst 0xa14006d9 // ldnt1b { z17.b, z25.b }, pn9.b/Z, [x22]\n" ".inst 0xa1410750 // ld1b { z16.b, z24.b }, pn9.b/Z, [x26, #0x2, MUL VL]\n" ".inst 0xa14106cb // ldnt1b { z3.b, z11.b }, pn9.b/Z, [x22, #0x2, MUL VL]\n" ".inst 0xa0420748 // ld1b { z8.b-z9.b }, pn9.b/Z, [x26, #0x4, MUL VL]\n" ".inst 0xa04206d3 // ldnt1b { z18.b-z19.b }, pn9.b/Z, [x22, #0x4, MUL VL]\n" ".inst 0xa0430744 // ld1b { z4.b-z5.b }, pn9.b/Z, [x26, #0x6, MUL VL]\n" "addvl x26, x26, #8\n" ".inst 0xa14306dd // ldnt1b { z21.b, z29.b }, pn9.b/Z, [x22, #0x6, MUL VL]\n" "addvl x22, x22, #8\n" "ble 7f\n" "6:" // K loop ".inst 0xa09102c0 // smopa za0.s, p0/M, p0/M, z22.b, z17.b\n" "subs x21, x21, #0x1\n" ".inst 0xa09902c1 // smopa za1.s, p0/M, p0/M, z22.b, z25.b\n" ".inst 0xa09103c2 // smopa za2.s, p0/M, p0/M, z30.b, z17.b\n" ".inst 0xa09903c3 // smopa za3.s, p0/M, p0/M, z30.b, z25.b\n" ".inst 0xa1400756 // ld1b { z22.b, z30.b }, pn9.b/Z, [x26]\n" ".inst 0xa0830200 // smopa za0.s, p0/M, p0/M, z16.b, z3.b\n" ".inst 0xa14006d9 // ldnt1b { z17.b, z25.b }, pn9.b/Z, [x22]\n" ".inst 0xa08b0201 // smopa za1.s, p0/M, p0/M, z16.b, z11.b\n" ".inst 0xa0830302 // smopa za2.s, p0/M, p0/M, z24.b, z3.b\n" ".inst 0xa08b0303 // smopa za3.s, p0/M, p0/M, z24.b, z11.b\n" ".inst 0xa1410750 // ld1b { z16.b, z24.b }, pn9.b/Z, [x26, #0x2, MUL VL]\n" ".inst 0xa0920100 // smopa za0.s, p0/M, p0/M, z8.b, z18.b\n" ".inst 0xa14106cb // ldnt1b { z3.b, z11.b }, pn9.b/Z, [x22, #0x2, MUL VL]\n" ".inst 0xa0930101 // smopa za1.s, p0/M, p0/M, z8.b, z19.b\n" ".inst 0xa0920122 // smopa za2.s, p0/M, p0/M, z9.b, z18.b\n" ".inst 0xa0930123 // smopa za3.s, p0/M, p0/M, z9.b, z19.b\n" ".inst 0xa0420748 // ld1b { z8.b-z9.b }, pn9.b/Z, [x26, #0x4, MUL VL]\n" ".inst 0xa04206d3 // ldnt1b { z18.b-z19.b }, pn9.b/Z, [x22, #0x4, MUL VL]\n" ".inst 0xa0950080 // smopa za0.s, p0/M, p0/M, z4.b, z21.b\n" ".inst 0xa09d0081 // smopa za1.s, p0/M, p0/M, z4.b, z29.b\n" ".inst 0xa09500a2 // smopa za2.s, p0/M, p0/M, z5.b, z21.b\n" ".inst 0xa09d00a3 // smopa za3.s, p0/M, p0/M, z5.b, z29.b\n" ".inst 0xa0430744 // ld1b { z4.b-z5.b }, pn9.b/Z, [x26, #0x6, MUL VL]\n" "addvl x26, x26, #8\n" ".inst 0xa14306dd // ldnt1b { z21.b, z29.b }, pn9.b/Z, [x22, #0x6, MUL VL]\n" "addvl x22, x22, #8\n" "bgt 6b\n" "7:" // K loop tail ".inst 0xa09102c0 // smopa za0.s, p0/M, p0/M, z22.b, z17.b\n" ".inst 0xa09902c1 // smopa za1.s, p0/M, p0/M, z22.b, z25.b\n" ".inst 0xa09103c2 // smopa za2.s, p0/M, p0/M, z30.b, z17.b\n" ".inst 0xa09903c3 // smopa za3.s, p0/M, p0/M, z30.b, z25.b\n" ".inst 0xa0830200 // smopa za0.s, p0/M, p0/M, z16.b, z3.b\n" ".inst 0xa08b0201 // smopa za1.s, p0/M, p0/M, z16.b, z11.b\n" ".inst 0xa0830302 // smopa za2.s, p0/M, p0/M, z24.b, z3.b\n" ".inst 0xa08b0303 // smopa za3.s, p0/M, p0/M, z24.b, z11.b\n" ".inst 0xa0920100 // smopa za0.s, p0/M, p0/M, z8.b, z18.b\n" ".inst 0xa0930101 // smopa za1.s, p0/M, p0/M, z8.b, z19.b\n" ".inst 0xa0920122 // smopa za2.s, p0/M, p0/M, z9.b, z18.b\n" ".inst 0xa0930123 // smopa za3.s, p0/M, p0/M, z9.b, z19.b\n" ".inst 0xa0950080 // smopa za0.s, p0/M, p0/M, z4.b, z21.b\n" ".inst 0xa09d0081 // smopa za1.s, p0/M, p0/M, z4.b, z29.b\n" ".inst 0xa09500a2 // smopa za2.s, p0/M, p0/M, z5.b, z21.b\n" ".inst 0xa09d00a3 // smopa za3.s, p0/M, p0/M, z5.b, z29.b\n" "8:" // K oddments "cbz x20, 10f\n" "9:" // K oddments: Loop ".inst 0xa1400756 // ld1b { z22.b, z30.b }, pn9.b/Z, [x26]\n" "subs x20, x20, #0x1\n" "addvl x26, x26, #2\n" ".inst 0xa14006d1 // ld1b { z17.b, z25.b }, pn9.b/Z, [x22]\n" "addvl x22, x22, #2\n" ".inst 0xa09102c0 // smopa za0.s, p0/M, p0/M, z22.b, z17.b\n" ".inst 0xa09902c1 // smopa za1.s, p0/M, p0/M, z22.b, z25.b\n" ".inst 0xa09103c2 // smopa za2.s, p0/M, p0/M, z30.b, z17.b\n" ".inst 0xa09903c3 // smopa za3.s, p0/M, p0/M, z30.b, z25.b\n" "bgt 9b\n" "10:" // K oddments: End "tbz x15, #1, 14f\n" "tbz x15, #0, 12f\n" "mov x12, #0x0\n" "cntw x19\n" "11:" // Store to partial result buffer: Store and refill: Loop ".inst 0xa040c5cc // ld1w { z12.s-z15.s }, pn9.b/Z, [x14]\n" ".inst 0xc0860418 // mova { z24.s-z27.s }, za0h.s[x12]\n" ".inst 0xc0840580 // mova za0h.s[x12], { z12.s-z15.s }\n" ".inst 0xc0860420 // mova { z0.s-z3.s }, za1h.s[x12]\n" ".inst 0xa041c5cc // ld1w { z12.s-z15.s }, pn9.b/Z, [x14, #0x4, MUL VL]\n" ".inst 0xc0840581 // mova za1h.s[x12], { z12.s-z15.s }\n" ".inst 0xc086045c // mova { z28.s-z31.s }, za2h.s[x12]\n" ".inst 0xc0860468 // mova { z8.s-z11.s }, za3h.s[x12]\n" ".inst 0xa042c5cc // ld1w { z12.s-z15.s }, pn9.b/Z, [x14, #0x8, MUL VL]\n" ".inst 0xc0840582 // mova za2h.s[x12], { z12.s-z15.s }\n" ".inst 0xa043c5c4 // ld1w { z4.s-z7.s }, pn9.b/Z, [x14, #0xc, MUL VL]\n" ".inst 0xc0840483 // mova za3h.s[x12], { z4.s-z7.s }\n" "add x12, x12, #0x4\n" "cmp x12, x19\n" ".inst 0xa060c5b8 // st1w { z24.s-z27.s }, pn9.b, [x13]\n" "addvl x14, x14, #16\n" ".inst 0xa061c5a0 // st1w { z0.s-z3.s }, pn9.b, [x13, #0x4, MUL VL]\n" ".inst 0xa062c5bc // st1w { z28.s-z31.s }, pn9.b, [x13, #0x8, MUL VL]\n" ".inst 0xa063c5a8 // st1w { z8.s-z11.s }, pn9.b, [x13, #0xc, MUL VL]\n" "addvl x13, x13, #16\n" "blt 11b\n" "b 23f\n" "12:" // Store to partial result buffer: Store only "mov x12, #0x0\n" "cntw x19\n" "13:" // Store to partial result buffer: Store only: Loop ".inst 0xc0860400 // mova { z0.s-z3.s }, za0h.s[x12]\n" ".inst 0xc086042c // mova { z12.s-z15.s }, za1h.s[x12]\n" ".inst 0xa060c5a0 // st1w { z0.s-z3.s }, pn9.b, [x13]\n" ".inst 0xc086045c // mova { z28.s-z31.s }, za2h.s[x12]\n" ".inst 0xc0860464 // mova { z4.s-z7.s }, za3h.s[x12]\n" ".inst 0xa061c5ac // st1w { z12.s-z15.s }, pn9.b, [x13, #0x4, MUL VL]\n" "add x12, x12, #0x4\n" "cmp x12, x19\n" ".inst 0xa062c5bc // st1w { z28.s-z31.s }, pn9.b, [x13, #0x8, MUL VL]\n" ".inst 0xa063c5a4 // st1w { z4.s-z7.s }, pn9.b, [x13, #0xc, MUL VL]\n" "addvl x13, x13, #16\n" "blt 13b\n" "b 23f\n" "14:" // Store to output array "ldr x25, [%x[args], %[offsetof_C]]\n" "sub x24, x11, x10\n" "cntw x23\n" "ldr x22, [%x[args], %[offsetof_ldcb]]\n" "cmp x24, x23\n" "csel x21, x24, x23, LT\n" "add x25, x25, x9, LSL #2\n" // C += n "lsr x20, x21, #0x2\n" "madd x25, x10, x22, x25\n" // C += m * ldc "mov x12, #0x0\n" "and x19, x21, #0x3\n" "cbz x20, 16f\n" "15:" // Store to output array: Accumulator row 0 loop ".inst 0xc0860410 // mova { z16.s-z19.s }, za0h.s[x12]\n" ".inst 0xc0860438 // mova { z24.s-z27.s }, za1h.s[x12]\n" ".inst 0xa1604330 // st1w { z16.s, z24.s }, p8, [x25]\n" "add x25, x25, x22\n" ".inst 0xa1604331 // st1w { z17.s, z25.s }, p8, [x25]\n" "add x25, x25, x22\n" "add x12, x12, #0x4\n" ".inst 0xa1604332 // st1w { z18.s, z26.s }, p8, [x25]\n" "add x25, x25, x22\n" "cmp x12, x20, LSL #2\n" ".inst 0xa1604333 // st1w { z19.s, z27.s }, p8, [x25]\n" "add x25, x25, x22\n" "blt 15b\n" "16:" // Store to output array: Accumulator row 0 oddments "cbz x19, 17f\n" "subs x19, x19, #0x1\n" ".inst 0xc0860400 // mova { z0.s-z3.s }, za0h.s[x12]\n" ".inst 0xc0860428 // mova { z8.s-z11.s }, za1h.s[x12]\n" ".inst 0xa1604320 // st1w { z0.s, z8.s }, p8, [x25]\n" "add x25, x25, x22\n" "beq 17f\n" "subs x19, x19, #0x1\n" ".inst 0xa1604321 // st1w { z1.s, z9.s }, p8, [x25]\n" "add x25, x25, x22\n" "beq 17f\n" ".inst 0xa1604322 // st1w { z2.s, z10.s }, p8, [x25]\n" "add x25, x25, x22\n" "17:" // Store to output array: Accumulator row 0 oddments: End "subs x24, x24, x21\n" "beq 21f\n" "cmp x24, x23\n" "csel x19, x24, x23, LT\n" "lsr x20, x19, #0x2\n" "mov x12, #0x0\n" "and x19, x19, #0x3\n" "cbz x20, 19f\n" "18:" // Store to output array: Accumulator row 1 loop ".inst 0xc0860450 // mova { z16.s-z19.s }, za2h.s[x12]\n" ".inst 0xc0860478 // mova { z24.s-z27.s }, za3h.s[x12]\n" ".inst 0xa1604330 // st1w { z16.s, z24.s }, p8, [x25]\n" "add x25, x25, x22\n" ".inst 0xa1604331 // st1w { z17.s, z25.s }, p8, [x25]\n" "add x25, x25, x22\n" "add x12, x12, #0x4\n" ".inst 0xa1604332 // st1w { z18.s, z26.s }, p8, [x25]\n" "add x25, x25, x22\n" "cmp x12, x20, LSL #2\n" ".inst 0xa1604333 // st1w { z19.s, z27.s }, p8, [x25]\n" "add x25, x25, x22\n" "blt 18b\n" "19:" // Store to output array: Accumulator row 1 oddments "cbz x19, 20f\n" "subs x19, x19, #0x1\n" ".inst 0xc0860444 // mova { z4.s-z7.s }, za2h.s[x12]\n" ".inst 0xc086046c // mova { z12.s-z15.s }, za3h.s[x12]\n" ".inst 0xa1604324 // st1w { z4.s, z12.s }, p8, [x25]\n" "add x25, x25, x22\n" "beq 20f\n" "subs x19, x19, #0x1\n" ".inst 0xa1604325 // st1w { z5.s, z13.s }, p8, [x25]\n" "add x25, x25, x22\n" "beq 20f\n" ".inst 0xa1604326 // st1w { z6.s, z14.s }, p8, [x25]\n" "20:" // Store to output array: Accumulator row 1 oddments: End "21:" // Store to output array: End "tbz x15, #0, 23f\n" "mov x12, #0x0\n" "cntw x19\n" "22:" // Store to output array: Refill accumulators: Loop ".inst 0xa040c5d0 // ld1w { z16.s-z19.s }, pn9.b/Z, [x14]\n" ".inst 0xc0840600 // mova za0h.s[x12], { z16.s-z19.s }\n" ".inst 0xa041c5d0 // ld1w { z16.s-z19.s }, pn9.b/Z, [x14, #0x4, MUL VL]\n" ".inst 0xc0840601 // mova za1h.s[x12], { z16.s-z19.s }\n" ".inst 0xa042c5c0 // ld1w { z0.s-z3.s }, pn9.b/Z, [x14, #0x8, MUL VL]\n" ".inst 0xc0840402 // mova za2h.s[x12], { z0.s-z3.s }\n" ".inst 0xa043c5cc // ld1w { z12.s-z15.s }, pn9.b/Z, [x14, #0xc, MUL VL]\n" ".inst 0xc0840583 // mova za3h.s[x12], { z12.s-z15.s }\n" "add x12, x12, #0x4\n" "cmp x12, x19\n" "addvl x14, x14, #16\n" "blt 22b\n" "23:" // End block "incw x9, ALL, MUL #2\n" "cmp x9, x28\n" "blt 3b\n" "incw x10, ALL, MUL #2\n" "cmp x10, x11\n" "mov x9, #0x0\n" "mov x27, x26\n" "blt 3b\n" ".inst 0xd503467f // SMSTOP\n" : : [args] "r" (&args), [offsetof_A] "I" (offsetof(KernelArgs, A)), [offsetof_B] "I" (offsetof(KernelArgs, B)), [offsetof_C] "I" (offsetof(KernelArgs, C)), [offsetof_K] "I" (offsetof(KernelArgs, K)), [offsetof_M] "I" (offsetof(KernelArgs, M)), [offsetof_N] "I" (offsetof(KernelArgs, N)), [offsetof_accumulator_buffer] "I" (offsetof(KernelArgs, accumulator_buffer)), [offsetof_bias] "I" (offsetof(KernelArgs, bias)), [offsetof_flags] "I" (offsetof(KernelArgs, flags)), [offsetof_kstride_bytes] "I" (offsetof(KernelArgs, kstride_bytes)), [offsetof_ldcb] "I" (offsetof(KernelArgs, ldcb)) : "cc", "memory", "p0", "p1", "p2", "p3", "p4", "p5", "p6", "p7", "p8", "p9", "p10", "p11", "p12", "p13", "p14", "p15", "x9", "x10", "x11", "x12", "x13", "x14", "x15", "x19", "x20", "x21", "x22", "x23", "x24", "x25", "x26", "x27", "x28", "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 arm_gemm #endif // ARM_COMPUTE_ENABLE_SME2 #endif // __ARM_FEATURE_SVE