Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 1 | // Copyright 2019 Google LLC |
| 2 | // |
| 3 | // This source code is licensed under the BSD-style license found in the |
| 4 | // LICENSE file in the root directory of this source tree. |
| 5 | |
| 6 | #include <xnnpack/assembly.h> |
| 7 | |
| 8 | .syntax unified |
| 9 | |
| 10 | // void xnn_f32_gemm_ukernel_4x8__aarch32_neon_ld64( |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 11 | // size_t mr, r0 |
| 12 | // size_t nc, r1 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 13 | // size_t kc, r2 -> r5 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 14 | // const uint8_t*restrict a, r3 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 15 | // size_t a_stride, sp + 96 -> (r7) |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 16 | // const void*restrict w, sp + 100 -> r9 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 17 | // uint8_t*restrict c, sp + 104 -> r11 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 18 | // size_t cm_stride, sp + 108 -> (r6) |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 19 | // size_t cn_stride, sp + 112 -> r7 |
Marat Dukhan | eb09a6b | 2020-04-08 17:34:32 -0700 | [diff] [blame] | 20 | // const union xnn_f32_minmax_params params[restrict static 1]) sp + 116 -> (r7) |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 21 | |
| 22 | |
| 23 | // inner loop registers |
| 24 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 25 | // A0 r3 d0 |
| 26 | // A1 r12 d1 |
| 27 | // A2 r10 d2 |
| 28 | // A3 r0 d3 |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 29 | |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 30 | // B r9 d8, d9, d10, d11 |
| 31 | // B d12, d13, d14, d15 |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 32 | |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 33 | // C0 r11 d16-d17 q8 d18-d19 q9 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 34 | // C1 r4 d20-d21 q10 d22-d23 q11 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 35 | // C2 r8 d24-d25 q12 d26-d27 q13 |
| 36 | // C3 r6 d28-d29 q14 d30-d31 q15 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 37 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 38 | // Clamp (r5) d4 d5 d6 d7 |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 39 | |
| 40 | BEGIN_FUNCTION xnn_f32_gemm_ukernel_4x8__aarch32_neon_ld64 |
| 41 | .arm |
| 42 | #ifndef __APPLE__ |
| 43 | .arch armv7-a |
| 44 | .fpu neon |
| 45 | #endif |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 46 | // Push 96 bytes |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 47 | PUSH {r4, r5, r6, r7, r8, r9, r10, r11} // 32 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 48 | VPUSH {d8-d15} // +64 = 96 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 49 | |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 50 | LDR r7, [sp, 96] // a_stride |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 51 | LDR r11, [sp, 104] // c |
| 52 | LDR r6, [sp, 108] // cm_stride |
Frank Barchard | 9b499d6 | 2020-02-12 11:05:00 -0800 | [diff] [blame] | 53 | LDR r9, [sp, 100] // w |
Marat Dukhan | eb09a6b | 2020-04-08 17:34:32 -0700 | [diff] [blame] | 54 | LDR r5, [sp, 116] // params |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 55 | |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 56 | // Clamp A and C pointers |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 57 | CMP r0, 2 // if mr >= 2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 58 | ADD r12, r3, r7 // a1 = a0 + a_stride |
| 59 | ADD r4, r11, r6 // c1 = c0 + cm_stride |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 60 | MOVLO r12, r3 // a1 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 61 | MOVLO r4, r11 // c1 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 62 | // if mr > 2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 63 | ADD r10, r12, r7 // a2 = a1 + a_stride |
| 64 | ADD r8, r4, r6 // c2 = c1 + cm_stride |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 65 | MOVLS r10, r12 // a2 |
| 66 | MOVLS r8, r4 // c2 |
| 67 | |
| 68 | CMP r0, 4 // if mr >=4 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 69 | ADD r0, r10, r7 // a3 = a2 + a_stride |
| 70 | ADD r6, r8, r6 // c3 = c2 + cm_stride |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 71 | MOVLO r0, r10 // a3 |
| 72 | MOVLO r6, r8 // c3 |
| 73 | |
Marat Dukhan | eb09a6b | 2020-04-08 17:34:32 -0700 | [diff] [blame] | 74 | // Load min/max values |
Frank Barchard | b738ad2 | 2019-12-10 00:02:50 -0800 | [diff] [blame] | 75 | VLD1.32 {d4[], d5[]}, [r5]! |
| 76 | LDR r7, [sp, 112] // cn_stride |
| 77 | VLD1.32 {d6[], d7[]}, [r5] |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 78 | |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 79 | 1: |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 80 | # Load initial bias from w into accumulators |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 81 | VLDM r9!, {d16-d19} // Bias |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 82 | SUBS r5, r2, 8 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 83 | VMOV q10, q8 |
| 84 | VMOV q11, q9 |
| 85 | VMOV q12, q8 |
| 86 | VMOV q13, q9 |
| 87 | VMOV q14, q8 |
| 88 | VMOV q15, q9 |
Frank Barchard | 9b499d6 | 2020-02-12 11:05:00 -0800 | [diff] [blame] | 89 | BLO 5f // less than 2 channels? |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 90 | |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 91 | // Main loop - 2 floats of A (8 bytes) |
| 92 | 2: |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 93 | VLD1.32 {d0}, [r3]! // A0 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 94 | VLDM r9!, {d8-d11} // B0 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 95 | VLD1.32 {d1}, [r12]! // A1 |
| 96 | VLD1.32 {d2}, [r10]! // A2 |
| 97 | VLD1.32 {d3}, [ r0]! // A3 |
| 98 | |
| 99 | VMLA.F32 q8, q4, d0[0] |
| 100 | VMLA.F32 q9, q5, d0[0] |
| 101 | VMLA.F32 q10, q4, d1[0] |
| 102 | VMLA.F32 q11, q5, d1[0] |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 103 | VLDM r9!, {d12-d15} // B1 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 104 | VMLA.F32 q12, q4, d2[0] |
| 105 | VMLA.F32 q13, q5, d2[0] |
| 106 | VMLA.F32 q14, q4, d3[0] |
| 107 | VMLA.F32 q15, q5, d3[0] |
| 108 | VMLA.F32 q8, q6, d0[1] |
| 109 | VMLA.F32 q9, q7, d0[1] |
| 110 | VMLA.F32 q10, q6, d1[1] |
| 111 | VMLA.F32 q11, q7, d1[1] |
| 112 | VMLA.F32 q12, q6, d2[1] |
| 113 | VMLA.F32 q13, q7, d2[1] |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 114 | SUBS r5, r5, 8 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 115 | VMLA.F32 q14, q6, d3[1] |
| 116 | VMLA.F32 q15, q7, d3[1] |
| 117 | BHS 2b |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 118 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 119 | // Is there a remainder?- 1 floats of A (4 bytes) |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 120 | TST r5, 4 |
Frank Barchard | 9b499d6 | 2020-02-12 11:05:00 -0800 | [diff] [blame] | 121 | BNE 5f |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 122 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 123 | 4: |
| 124 | // Clamp |
Marat Dukhan | a51cf48 | 2020-04-08 16:16:19 -0700 | [diff] [blame] | 125 | VMAX.F32 q8, q8, q2 |
Frank Barchard | 0090f5b | 2019-12-16 17:02:57 -0800 | [diff] [blame] | 126 | SUBS r1, r1, 8 |
Marat Dukhan | a51cf48 | 2020-04-08 16:16:19 -0700 | [diff] [blame] | 127 | VMAX.F32 q9, q9, q2 |
| 128 | VMAX.F32 q10, q10, q2 |
| 129 | VMAX.F32 q11, q11, q2 |
| 130 | VMAX.F32 q12, q12, q2 |
| 131 | VMAX.F32 q13, q13, q2 |
| 132 | VMAX.F32 q14, q14, q2 |
| 133 | VMAX.F32 q15, q15, q2 |
| 134 | VMIN.F32 q8, q8, q3 |
| 135 | VMIN.F32 q9, q9, q3 |
| 136 | VMIN.F32 q10, q10, q3 |
| 137 | VMIN.F32 q11, q11, q3 |
| 138 | VMIN.F32 q12, q12, q3 |
| 139 | VMIN.F32 q13, q13, q3 |
| 140 | VMIN.F32 q14, q14, q3 |
| 141 | VMIN.F32 q15, q15, q3 |
Frank Barchard | 72d6afb | 2019-11-30 19:21:26 -0800 | [diff] [blame] | 142 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 143 | // Store full 4 x 8 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 144 | BLO 10f |
| 145 | VST1.32 {d16-d19}, [r11], r7 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 146 | SUB r0, r0, r2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 147 | VST1.32 {d20-d23}, [r4], r7 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 148 | SUB r10, r10, r2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 149 | VST1.32 {d24-d27}, [r8], r7 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 150 | SUB r12, r12, r2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 151 | VST1.32 {d28-d31}, [r6], r7 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 152 | SUB r3, r3, r2 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 153 | BHI 1b |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 154 | |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 155 | VPOP {d8-d15} |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 156 | POP {r4, r5, r6, r7, r8, r9, r10, r11} |
| 157 | BX lr |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 158 | |
Frank Barchard | 9b499d6 | 2020-02-12 11:05:00 -0800 | [diff] [blame] | 159 | 5: |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 160 | // Remainder- 1 floats of A (4 bytes) |
| 161 | VLDM r3!, {s0} // A0 |
| 162 | VLDM r9!, {d8-d11} // B0 |
| 163 | VLDM r12!, {s2} // A1 |
| 164 | VLDM r10!, {s4} // A2 |
| 165 | VLDM r0!, {s6} // A3 |
| 166 | VMLA.F32 q8, q4, d0[0] |
| 167 | VMLA.F32 q9, q5, d0[0] |
| 168 | VMLA.F32 q10, q4, d1[0] |
| 169 | VMLA.F32 q11, q5, d1[0] |
| 170 | VMLA.F32 q12, q4, d2[0] |
| 171 | VMLA.F32 q13, q5, d2[0] |
| 172 | VMLA.F32 q14, q4, d3[0] |
| 173 | VMLA.F32 q15, q5, d3[0] |
| 174 | B 4b |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 175 | |
Frank Barchard | 9b499d6 | 2020-02-12 11:05:00 -0800 | [diff] [blame] | 176 | // Store odd width |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 177 | 10: |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 178 | TST r1, 4 |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 179 | BEQ 11f |
| 180 | VST1.32 {d16-d17}, [r11]! |
| 181 | VMOV q8, q9 |
| 182 | VST1.32 {d20-d21}, [r4]! |
| 183 | VMOV q10, q11 |
| 184 | VST1.32 {d24-d25}, [r8]! |
| 185 | VMOV q12, q13 |
| 186 | VST1.32 {d28-d29}, [r6]! |
| 187 | VMOV q14, q15 |
Frank Barchard | 61cad89 | 2019-12-02 13:06:15 -0800 | [diff] [blame] | 188 | |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 189 | 11: |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 190 | TST r1, 2 |
| 191 | BEQ 12f |
| 192 | VST1.32 {d16}, [r11]! |
| 193 | VMOV d16, d17 |
| 194 | VST1.32 {d20}, [r4]! |
| 195 | VMOV d20, d21 |
| 196 | VST1.32 {d24}, [r8]! |
| 197 | VMOV d24, d25 |
| 198 | VST1.32 {d28}, [r6]! |
| 199 | VMOV d28, d29 |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 200 | |
| 201 | 12: |
Frank Barchard | f917cbd | 2019-12-04 23:02:08 -0800 | [diff] [blame] | 202 | TST r1, 1 |
| 203 | BEQ 13f |
| 204 | VST1.32 {d16[0]}, [r11] |
| 205 | VST1.32 {d20[0]}, [r4] |
| 206 | VST1.32 {d24[0]}, [r8] |
| 207 | VST1.32 {d28[0]}, [r6] |
| 208 | |
| 209 | 13: |
| 210 | VPOP {d8-d15} |
| 211 | POP {r4, r5, r6, r7, r8, r9, r10, r11} |
| 212 | BX lr |
Frank Barchard | 8b0f026 | 2019-11-27 23:18:40 -0800 | [diff] [blame] | 213 | |
| 214 | END_FUNCTION xnn_f32_gemm_ukernel_4x8__aarch32_neon_ld64 |
| 215 | |
| 216 | #ifdef __ELF__ |
| 217 | .section ".note.GNU-stack","",%progbits |
| 218 | #endif |
| 219 | |
| 220 | |
| 221 | |