1 /*
2 * Copyright (c) 2022 Arm Limited.
3 *
4 * SPDX-License-Identifier: MIT
5 *
6 * Permission is hereby granted, free of charge, to any person obtaining a copy
7 * of this software and associated documentation files (the "Software"), to
8 * deal in the Software without restriction, including without limitation the
9 * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
10 * sell copies of the Software, and to permit persons to whom the Software is
11 * furnished to do so, subject to the following conditions:
12 *
13 * The above copyright notice and this permission notice shall be included in
14 * all copies or substantial portions of the Software.
15 *
16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
17 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
18 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
19 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
20 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
21 * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS
22 * IN THE SOFTWARE.
23 */
24 #if defined(__aarch64__) && (defined(FP16_KERNELS) || defined(__ARM_FEATURE_FP16_VECTOR_ARITHMETIC))
25
26 #include <cstddef>
27
28 namespace arm_gemm {
29
a64_ffinterleaved_fp16_mla_8x24(const __fp16 * Apanel,const __fp16 * Bpanel,size_t B_stride,__fp16 * Cpanel,int ablocks,size_t N,int K)30 void a64_ffinterleaved_fp16_mla_8x24(
31 const __fp16 *Apanel,
32 const __fp16 *Bpanel,
33 size_t B_stride,
34 __fp16 *Cpanel,
35 int ablocks,
36 size_t N,
37 int K) {
38
39 struct KernelArgs {
40 size_t K = {};
41 const __fp16 *Bpanel = {};
42 size_t N = {};
43 size_t B_stride = {};
44 const __fp16 *cur_B_ptr = {};
45 } ka;
46
47 ka.K = (K/1) - 1;
48 ka.Bpanel = Bpanel;
49 ka.N = N;
50 ka.B_stride = B_stride;
51
52 __asm__ __volatile__(
53 "1:" // Height loop
54 "ldr x24, [%x[args_ptr], %[offsetof_Bpanel]]\n"
55 "ldr x23, [%x[args_ptr], %[offsetof_N]]\n"
56 "str x24, [%x[args_ptr], %[offsetof_cur_B_ptr]]\n"
57 "mov x22, %x[Apanel]\n"
58 "2:" // Width loop
59 "ldr x24, [%x[args_ptr], %[offsetof_cur_B_ptr]]\n"
60 "ldr x19, [%x[args_ptr], %[offsetof_B_stride]]\n"
61 "add x21, x24, x19, LSL #1\n"
62 "add x20, x21, x19, LSL #1\n"
63 "add x19, x20, x19, LSL #1\n"
64 "str x19, [%x[args_ptr], %[offsetof_cur_B_ptr]]\n"
65 "cmp x23, #0x10\n"
66 "mov %x[Apanel], x22\n"
67 "bgt 3f\n"
68 "cmp x23, #0x8\n"
69 "mov x20, x24\n"
70 "bgt 3f\n"
71 "mov x21, x24\n"
72 "3:" // B setup done
73 "ldr q0, [%x[Apanel], #0x0]\n"
74 "ldr q2, [x24, #0x0]\n"
75 "movi v8.16b, #0x0\n"
76 "ldr q3, [x21, #0x0]\n"
77 "ldr q4, [x20, #0x0]\n"
78 "movi v9.16b, #0x0\n"
79 "ldr x19, [%x[args_ptr], %[offsetof_K]]\n"
80 "cmp x19, #0x2\n"
81 "movi v10.16b, #0x0\n"
82 "movi v11.16b, #0x0\n"
83 "movi v12.16b, #0x0\n"
84 "movi v13.16b, #0x0\n"
85 "movi v14.16b, #0x0\n"
86 "movi v15.16b, #0x0\n"
87 "movi v16.16b, #0x0\n"
88 "movi v17.16b, #0x0\n"
89 "movi v18.16b, #0x0\n"
90 "movi v19.16b, #0x0\n"
91 "movi v20.16b, #0x0\n"
92 "movi v21.16b, #0x0\n"
93 "movi v22.16b, #0x0\n"
94 "movi v23.16b, #0x0\n"
95 "movi v24.16b, #0x0\n"
96 "movi v25.16b, #0x0\n"
97 "movi v26.16b, #0x0\n"
98 "movi v27.16b, #0x0\n"
99 "movi v28.16b, #0x0\n"
100 "movi v29.16b, #0x0\n"
101 "movi v30.16b, #0x0\n"
102 "movi v31.16b, #0x0\n"
103 "blt 5f\n"
104 "4:" // main loop head
105 "ldr q1, [%x[Apanel], #0x10]\n"
106 "ldr q5, [x24, #0x10]\n"
107 "fmla v8.8h, v2.8h, v0.h[0]\n"
108 "ldr q6, [x21, #0x10]\n"
109 "ldr q7, [x20, #0x10]\n"
110 "fmla v11.8h, v2.8h, v0.h[1]\n"
111 "fmla v14.8h, v2.8h, v0.h[2]\n"
112 "fmla v17.8h, v2.8h, v0.h[3]\n"
113 "sub x19, x19, #0x2\n"
114 "fmla v20.8h, v2.8h, v0.h[4]\n"
115 "fmla v23.8h, v2.8h, v0.h[5]\n"
116 "cmp x19, #0x2\n"
117 "fmla v26.8h, v2.8h, v0.h[6]\n"
118 "fmla v29.8h, v2.8h, v0.h[7]\n"
119 "add %x[Apanel], %x[Apanel], #0x20\n"
120 "fmla v9.8h, v3.8h, v0.h[0]\n"
121 "fmla v12.8h, v3.8h, v0.h[1]\n"
122 "add x24, x24, #0x20\n"
123 "ldr q2, [x24, #0x0]\n"
124 "fmla v15.8h, v3.8h, v0.h[2]\n"
125 "fmla v18.8h, v3.8h, v0.h[3]\n"
126 "fmla v21.8h, v3.8h, v0.h[4]\n"
127 "fmla v24.8h, v3.8h, v0.h[5]\n"
128 "add x21, x21, #0x20\n"
129 "fmla v27.8h, v3.8h, v0.h[6]\n"
130 "fmla v30.8h, v3.8h, v0.h[7]\n"
131 "ldr q3, [x21, #0x0]\n"
132 "fmla v10.8h, v4.8h, v0.h[0]\n"
133 "fmla v13.8h, v4.8h, v0.h[1]\n"
134 "add x20, x20, #0x20\n"
135 "fmla v16.8h, v4.8h, v0.h[2]\n"
136 "fmla v19.8h, v4.8h, v0.h[3]\n"
137 "fmla v22.8h, v4.8h, v0.h[4]\n"
138 "fmla v25.8h, v4.8h, v0.h[5]\n"
139 "fmla v28.8h, v4.8h, v0.h[6]\n"
140 "fmla v31.8h, v4.8h, v0.h[7]\n"
141 "ldr q0, [%x[Apanel], #0x0]\n"
142 "ldr q4, [x20, #0x0]\n"
143 "fmla v8.8h, v5.8h, v1.h[0]\n"
144 "fmla v11.8h, v5.8h, v1.h[1]\n"
145 "fmla v14.8h, v5.8h, v1.h[2]\n"
146 "fmla v17.8h, v5.8h, v1.h[3]\n"
147 "fmla v20.8h, v5.8h, v1.h[4]\n"
148 "fmla v23.8h, v5.8h, v1.h[5]\n"
149 "fmla v26.8h, v5.8h, v1.h[6]\n"
150 "fmla v29.8h, v5.8h, v1.h[7]\n"
151 "fmla v9.8h, v6.8h, v1.h[0]\n"
152 "fmla v12.8h, v6.8h, v1.h[1]\n"
153 "fmla v15.8h, v6.8h, v1.h[2]\n"
154 "fmla v18.8h, v6.8h, v1.h[3]\n"
155 "fmla v21.8h, v6.8h, v1.h[4]\n"
156 "fmla v24.8h, v6.8h, v1.h[5]\n"
157 "fmla v27.8h, v6.8h, v1.h[6]\n"
158 "fmla v30.8h, v6.8h, v1.h[7]\n"
159 "fmla v10.8h, v7.8h, v1.h[0]\n"
160 "fmla v13.8h, v7.8h, v1.h[1]\n"
161 "fmla v16.8h, v7.8h, v1.h[2]\n"
162 "fmla v19.8h, v7.8h, v1.h[3]\n"
163 "fmla v22.8h, v7.8h, v1.h[4]\n"
164 "fmla v25.8h, v7.8h, v1.h[5]\n"
165 "fmla v28.8h, v7.8h, v1.h[6]\n"
166 "fmla v31.8h, v7.8h, v1.h[7]\n"
167 "bge 4b\n"
168 "5:" // main loop skip
169 "fmla v8.8h, v2.8h, v0.h[0]\n"
170 "fmla v11.8h, v2.8h, v0.h[1]\n"
171 "add %x[Apanel], %x[Apanel], #0x10\n"
172 "fmla v14.8h, v2.8h, v0.h[2]\n"
173 "fmla v17.8h, v2.8h, v0.h[3]\n"
174 "add x24, x24, #0x10\n"
175 "fmla v20.8h, v2.8h, v0.h[4]\n"
176 "fmla v23.8h, v2.8h, v0.h[5]\n"
177 "add x21, x21, #0x10\n"
178 "fmla v26.8h, v2.8h, v0.h[6]\n"
179 "fmla v29.8h, v2.8h, v0.h[7]\n"
180 "add x20, x20, #0x10\n"
181 "fmla v9.8h, v3.8h, v0.h[0]\n"
182 "fmla v12.8h, v3.8h, v0.h[1]\n"
183 "fmla v15.8h, v3.8h, v0.h[2]\n"
184 "fmla v18.8h, v3.8h, v0.h[3]\n"
185 "fmla v21.8h, v3.8h, v0.h[4]\n"
186 "fmla v24.8h, v3.8h, v0.h[5]\n"
187 "fmla v27.8h, v3.8h, v0.h[6]\n"
188 "fmla v30.8h, v3.8h, v0.h[7]\n"
189 "fmla v10.8h, v4.8h, v0.h[0]\n"
190 "fmla v13.8h, v4.8h, v0.h[1]\n"
191 "fmla v16.8h, v4.8h, v0.h[2]\n"
192 "fmla v19.8h, v4.8h, v0.h[3]\n"
193 "fmla v22.8h, v4.8h, v0.h[4]\n"
194 "fmla v25.8h, v4.8h, v0.h[5]\n"
195 "fmla v28.8h, v4.8h, v0.h[6]\n"
196 "fmla v31.8h, v4.8h, v0.h[7]\n"
197 "cbz x19, 6f\n"
198 "ldr q0, [%x[Apanel], #0x0]\n"
199 "ldr q5, [x24, #0x0]\n"
200 "fmla v8.8h, v5.8h, v0.h[0]\n"
201 "ldr q6, [x21, #0x0]\n"
202 "ldr q7, [x20, #0x0]\n"
203 "fmla v11.8h, v5.8h, v0.h[1]\n"
204 "fmla v14.8h, v5.8h, v0.h[2]\n"
205 "fmla v17.8h, v5.8h, v0.h[3]\n"
206 "add %x[Apanel], %x[Apanel], #0x10\n"
207 "fmla v20.8h, v5.8h, v0.h[4]\n"
208 "fmla v23.8h, v5.8h, v0.h[5]\n"
209 "fmla v26.8h, v5.8h, v0.h[6]\n"
210 "fmla v29.8h, v5.8h, v0.h[7]\n"
211 "fmla v9.8h, v6.8h, v0.h[0]\n"
212 "fmla v12.8h, v6.8h, v0.h[1]\n"
213 "fmla v15.8h, v6.8h, v0.h[2]\n"
214 "fmla v18.8h, v6.8h, v0.h[3]\n"
215 "fmla v21.8h, v6.8h, v0.h[4]\n"
216 "fmla v24.8h, v6.8h, v0.h[5]\n"
217 "fmla v27.8h, v6.8h, v0.h[6]\n"
218 "fmla v30.8h, v6.8h, v0.h[7]\n"
219 "fmla v10.8h, v7.8h, v0.h[0]\n"
220 "fmla v13.8h, v7.8h, v0.h[1]\n"
221 "fmla v16.8h, v7.8h, v0.h[2]\n"
222 "fmla v19.8h, v7.8h, v0.h[3]\n"
223 "fmla v22.8h, v7.8h, v0.h[4]\n"
224 "fmla v25.8h, v7.8h, v0.h[5]\n"
225 "fmla v28.8h, v7.8h, v0.h[6]\n"
226 "fmla v31.8h, v7.8h, v0.h[7]\n"
227 "6:" // multiply loop done
228 "subs x23, x23, #0x18\n"
229 "str q8, [%x[Cpanel], #0x0]\n"
230 "str q9, [%x[Cpanel], #0x10]\n"
231 "str q10, [%x[Cpanel], #0x20]\n"
232 "str q11, [%x[Cpanel], #0x30]\n"
233 "str q12, [%x[Cpanel], #0x40]\n"
234 "str q13, [%x[Cpanel], #0x50]\n"
235 "str q14, [%x[Cpanel], #0x60]\n"
236 "str q15, [%x[Cpanel], #0x70]\n"
237 "str q16, [%x[Cpanel], #0x80]\n"
238 "str q17, [%x[Cpanel], #0x90]\n"
239 "str q18, [%x[Cpanel], #0xa0]\n"
240 "str q19, [%x[Cpanel], #0xb0]\n"
241 "str q20, [%x[Cpanel], #0xc0]\n"
242 "str q21, [%x[Cpanel], #0xd0]\n"
243 "str q22, [%x[Cpanel], #0xe0]\n"
244 "str q23, [%x[Cpanel], #0xf0]\n"
245 "str q24, [%x[Cpanel], #0x100]\n"
246 "str q25, [%x[Cpanel], #0x110]\n"
247 "str q26, [%x[Cpanel], #0x120]\n"
248 "str q27, [%x[Cpanel], #0x130]\n"
249 "str q28, [%x[Cpanel], #0x140]\n"
250 "str q29, [%x[Cpanel], #0x150]\n"
251 "str q30, [%x[Cpanel], #0x160]\n"
252 "str q31, [%x[Cpanel], #0x170]\n"
253 "add %x[Cpanel], %x[Cpanel], #0x180\n"
254 "bgt 2b\n"
255 "subs %x[ablocks], %x[ablocks], #0x1\n"
256 "bne 1b\n"
257 : [Apanel] "+&r" (Apanel), [Cpanel] "+&r" (Cpanel), [ablocks] "+&r" (ablocks)
258 : [args_ptr] "r" (&ka), [offsetof_B_stride] "I" (offsetof(KernelArgs, B_stride)), [offsetof_Bpanel] "I" (offsetof(KernelArgs, Bpanel)), [offsetof_K] "I" (offsetof(KernelArgs, K)), [offsetof_N] "I" (offsetof(KernelArgs, N)), [offsetof_cur_B_ptr] "I" (offsetof(KernelArgs, cur_B_ptr))
259 : "cc", "memory", "v0", "v1", "v2", "v3", "v4", "v5", "v6", "v7", "v8", "v9", "v10", "v11", "v12", "v13", "v14", "v15", "v16", "v17", "v18", "v19", "v20", "v21", "v22", "v23", "v24", "v25", "v26", "v27", "v28", "v29", "v30", "v31", "x19", "x20", "x21", "x22", "x23", "x24"
260 );
261 }
262
263 } // namespace arm_gemm
264 #endif // defined(__aarch64__) && (defined(FP16_KERNELS) || defined(__ARM_FEATURE_FP16_VECTOR_ARITHMETIC))
265