// Copyright 2019 Google LLC // // This source code is licensed under the BSD-style license found in the // LICENSE file in the root directory of this source tree. $assert DATATYPE in ["F32", "QC8", "QC4"] #include "xnnpack/assembly.h" $DATATYPE_SPEC = {"F32": "f32", "QC8": "f32_qc8w", "QC4": "f32_qc4w"}[DATATYPE] # void xnn_${DATATYPE_SPEC}_gemm${"_goi" if GOI else ""}_minmax_ukernel_1x8__asm_aarch64_neonfma_ld128${"_prfm" if PREFETCH else ""}( # size_t mr, (x0) - unused. mr = 1 # size_t nc, x1 # size_t kc, x2 / x0 # const float* a, x3 # size_t a_stride, (x4) - unused # const void* w, x5 # float* c, x6 # size_t cm_stride, (x7) - unused # size_t cn_stride, [sp] -> x14 $if DATATYPE == "QC4": # const xnn_f32_qc4w_minmax_params* params) [sp + 8] -> (x8) $else: # const xnn_f32_minmax_params* params) [sp + 8] -> (x8) # d8-d15, x19-x30 need to be preserved if used. x18 is reserved by the OS. # Register usage # A0 x3 v0 # B x5 v20 v24 v21 v25 v22 v26 v23 v27 # C0 x6 v16 v17 # Clamp v4 v5 $if DATATYPE == "QC4": # zerop/mask v6 v7 $if GOI: # B stride x16, x17 BEGIN_FUNCTION xnn_${DATATYPE_SPEC}_gemm${"_goi" if GOI else ""}_minmax_ukernel_1x8__asm_aarch64_neonfma_ld128${"_prfm" if PREFETCH else ""} # Load cn_stride, params pointer LDP x14, x8, [sp] $if DATATYPE == "QC4": # Load min/max/zerop values LD3R {v4.4s, v5.4s, v6.4s}, [x8] NEG v6.4s, v6.4s MOVI v7.8b, 15 $else: # Load min/max values LD2R {v4.4s, v5.4s}, [x8] $if GOI: SUB x16, x2, x2, LSL #3 // B stride on 8th row to top (16 - KC * 7) ADD x16, x16, 16 LSL x17, x2, 3 // B stride for N loop SUB x17, x17, x2 0: $if GOI: # Zero bias for GOI MOVI v16.16b, 0 MOVI v17.16b, 0 $else: # Load initial bias from w into accumulators LDP q16, q17, [x5], 32 # Is there at least 4 floats (16 bytes) SUBS x0, x2, 16 // k = kc - 16 B.LO 3f $if PREFETCH: PRFM PLDL1KEEP, [x5] PRFM PLDL1KEEP, [x5, 64] PRFM PLDL1KEEP, [x5, 128] # Main loop - 4 floats of A (16 bytes) 1: $if GOI: LDR q0, [x3], 16 LD4 {v20.s, v21.s, v22.s, v23.s}[0], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v20.s, v21.s, v22.s, v23.s}[1], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v20.s, v21.s, v22.s, v23.s}[2], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v20.s, v21.s, v22.s, v23.s}[3], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v24.s, v25.s, v26.s, v27.s}[0], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v24.s, v25.s, v26.s, v27.s}[1], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v24.s, v25.s, v26.s, v27.s}[2], [x5], x2 $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] LD4 {v24.s, v25.s, v26.s, v27.s}[3], [x5], x16 $else: $if DATATYPE == "F32": LDR q0, [x3], 16 LDP q20, q24, [x5], 32 LDP q21, q25, [x5], 32 $else: LDR q21, [x5], 16 SXTL v24.8h, v21.8b SXTL2 v25.8h, v21.16b LDR q0, [x3], 16 SXTL v20.4s, v24.4h SXTL v21.4s, v25.4h SXTL2 v24.4s, v24.8h SXTL2 v25.4s, v25.8h SCVTF v20.4s, v20.4s SCVTF v24.4s, v24.4s SCVTF v21.4s, v21.4s SCVTF v25.4s, v25.4s FMLA v16.4s, v20.4s, v0.s[0] FMLA v17.4s, v24.4s, v0.s[0] $if PREFETCH: PRFM PLDL1KEEP, [x5, 128] FMLA v16.4s, v21.4s, v0.s[1] FMLA v17.4s, v25.4s, v0.s[1] $if not GOI: $if DATATYPE == "F32": LDP q22, q26, [x5], 32 LDP q23, q27, [x5], 32 $else: LDR q23, [x5], 16 SXTL v26.8h, v23.8b SXTL2 v27.8h, v23.16b SXTL v22.4s, v26.4h SXTL v23.4s, v27.4h SXTL2 v26.4s, v26.8h SXTL2 v27.4s, v27.8h SCVTF v22.4s, v22.4s SCVTF v26.4s, v26.4s SCVTF v23.4s, v23.4s SCVTF v27.4s, v27.4s SUBS x0, x0, 16 FMLA v16.4s, v22.4s, v0.s[2] FMLA v17.4s, v26.4s, v0.s[2] $if PREFETCH and (DATATYPE == "F32") and not GOI: PRFM PLDL1KEEP, [x5, 128] FMLA v16.4s, v23.4s, v0.s[3] FMLA v17.4s, v27.4s, v0.s[3] B.HS 1b # Is there a remainder?- 2 float of A (8 bytes) TBNZ x0, 3, 4f # Is there a remainder?- 1 float of A (4 bytes) TBNZ x0, 2, 5f 2: $if DATATYPE in ["QC8", "QC4"]: # Scale LDP q22, q26, [x5], 32 FMUL v16.4s, v16.4s, v22.4s FMUL v17.4s, v17.4s, v26.4s SUBS x1, x1, 8 # Clamp FMAX v16.4s, v16.4s, v4.4s FMAX v17.4s, v17.4s, v4.4s FMIN v16.4s, v16.4s, v5.4s FMIN v17.4s, v17.4s, v5.4s # Store full 1 x 8 B.LO 6f $if GOI: ADD x5, x5, x17 // advance to next row of weights STP q16, q17, [x6] ADD x6, x6, x14 SUB x3, x3, x2 // a0 -= kc B.HI 0b RET 3: TBZ x0, 3, 5f # Remainder- 2 float of A (4 bytes) 4: $if GOI: LD2 {v20.s, v21.s}[0], [x5], x2 LD2 {v20.s, v21.s}[1], [x5], x2 LD2 {v20.s, v21.s}[2], [x5], x2 LD2 {v20.s, v21.s}[3], [x5], x2 LD2 {v24.s, v25.s}[0], [x5], x2 LD2 {v24.s, v25.s}[1], [x5], x2 LD2 {v24.s, v25.s}[2], [x5], x2 LD2 {v24.s, v25.s}[3], [x5], x16 SUB x5, x5, 8 $else: # Remainder- 2 floats of A (8 bytes) $if DATATYPE == "F32": LDP q20, q24, [x5], 32 LDP q21, q25, [x5], 32 $else: LDP d24, d25, [x5], 16 SXTL v24.8h, v24.8b SXTL v20.4s, v24.4h SXTL2 v24.4s, v24.8h SCVTF v20.4s, v20.4s SCVTF v24.4s, v24.4s SXTL v25.8h, v25.8b SXTL v21.4s, v25.4h SXTL2 v25.4s, v25.8h SCVTF v21.4s, v21.4s SCVTF v25.4s, v25.4s LDR d0, [x3], 8 FMLA v16.4s, v20.4s, v0.s[0] FMLA v17.4s, v24.4s, v0.s[0] FMLA v16.4s, v21.4s, v0.s[1] FMLA v17.4s, v25.4s, v0.s[1] TBZ x0, 2, 2b # Remainder- 1 float of A (4 bytes) 5: $if GOI: LD1 {v20.s}[0], [x5], x2 LD1 {v20.s}[1], [x5], x2 LD1 {v20.s}[2], [x5], x2 LD1 {v20.s}[3], [x5], x2 LD1 {v24.s}[0], [x5], x2 LD1 {v24.s}[1], [x5], x2 LD1 {v24.s}[2], [x5], x2 LD1 {v24.s}[3], [x5], x16 SUB x5, x5, 12 $else: # Remainder- 2 floats of A (8 bytes) $if DATATYPE == "F32": LDP q20, q24, [x5], 32 $else: LDR d24, [x5], 8 SXTL v24.8h, v24.8b SXTL v20.4s, v24.4h SXTL2 v24.4s, v24.8h SCVTF v20.4s, v20.4s SCVTF v24.4s, v24.4s LDR s0, [x3], 4 FMLA v16.4s, v20.4s, v0.s[0] FMLA v17.4s, v24.4s, v0.s[0] B 2b # Store odd channels 6: TBZ x1, 2, 7f STR q16, [x6], 16 MOV v16.16b, v17.16b 7: TBZ x1, 1, 8f STR d16, [x6], 8 DUP d16, v16.d[1] 8: TBZ x1, 0, 9f STR s16, [x6] 9: RET END_FUNCTION xnn_${DATATYPE_SPEC}_gemm${"_goi" if GOI else ""}_minmax_ukernel_1x8__asm_aarch64_neonfma_ld128${"_prfm" if PREFETCH else ""} #ifdef __ELF__ .section ".note.GNU-stack","",%progbits #endif