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 
25 #if defined(__ARM_FEATURE_SVE)
26 
27 template <>
interleave_block(bfloat16 * & out,const bfloat16 * const * in,size_t width,size_t height,size_t row_offset,bool first)28 void interleave_block<2, 2, VLType::SME, false>(
29   bfloat16 * &out, const bfloat16 * const *in,
30   size_t width, size_t height, size_t row_offset, bool first
31 )
32 {
33   ARM_COMPUTE_UNUSED(first);
34 
35   __asm__ __volatile__(
36       ".inst 0xd503477f  // SMSTART ZA\n"
37       "cnth x20\n"
38       "cntw x16\n"
39       "cntw x15, ALL, MUL #2\n"
40       "cntw x14, ALL, MUL #3\n"
41       "mov x19, %x[width]\n"
42       "inch x19\n"
43       "sub x19, x19, #0x1\n"
44       "udiv x19, x19, x20\n" // n_passes = ceildiv(width, VL<T>)
45       "sub x13, x19, #0x1\n"
46       "lsr x13, x13, #0x1\n" // n_loops = (n_passes - 1) / 2
47       "and x11, x19, #0x1\n" // odd_tail = bool(n_passes & 0x1)
48       "mov x19, %x[width]\n"
49       "sub x10, x20, #0x1\n"
50       "ands x10, x19, x10\n"
51       "csel x10, x10, x20, NE\n"
52       "add x10, x10, #0x1\n"
53       "lsr x10, x10, #0x1\n"
54       "sub x9, x16, #0x2\n"
55       "ptrue p13.s\n"
56       "lsl x20, %x[height], #0x1\n" // height * 2
57       "lsl x19, x16, #0x1\n"
58       "whilelt p12.h, XZR, x20\n"
59       "whilelt p11.h, x19, x20\n"
60       "mov x28, %x[row_offset]\n"
61       "mov x27, %x[out]\n"
62       "mov x26, #0x0\n"
63       "whilelt p10.h, x26, %x[width]\n"
64       "whilelt p9.h, x26, %x[width]\n"
65       "whilelt p8.h, x26, %x[width]\n"
66       "mov x25, %x[in]\n"
67       "add x24, %x[in], x16, LSL #3\n"
68       "ldr x23, [x25, #0x0]\n"
69       "ldr x22, [x24, #0x0]\n"
70       "ldr x21, [x25, #0x8]\n"
71       "ldr x20, [x24, #0x8]\n"
72       "add x25, x25, #0x10\n"
73       "add x24, x24, #0x10\n"
74       "mov x12, #0x0\n"
75       "cbz x9, 2f\n"
76       "1:"  // K loop: Charge: Loop
77       ".inst 0x25286580  // dup p0.h, p9.h/Z, p12.h[w12]\n"
78       ".inst 0xe05c02e0  // ld1h { za0h.h[x12] }, p0/Z, [x23, x28, LSL #1]\n"
79       ".inst 0x25286160  // dup p0.h, p8.h/Z, p11.h[w12]\n"
80       ".inst 0xe05c02c8  // ld1h { za1h.h[x12] }, p0/Z, [x22, x28, LSL #1]\n"
81       ".inst 0x25686580  // dup p0.h, p9.h/Z, p12.h[w12, #2]\n"
82       ".inst 0xe05c02a2  // ld1h { za0h.h[x12, #2] }, p0/Z, [x21, x28, LSL #1]\n"
83       ".inst 0x25686160  // dup p0.h, p8.h/Z, p11.h[w12, #2]\n"
84       ".inst 0xe05c028a  // ld1h { za1h.h[x12, #2] }, p0/Z, [x20, x28, LSL #1]\n"
85       "ldr x23, [x25, #0x0]\n"
86       "ldr x22, [x24, #0x0]\n"
87       "ldr x21, [x25, #0x8]\n"
88       "ldr x20, [x24, #0x8]\n"
89       "add x25, x25, #0x10\n"
90       "add x24, x24, #0x10\n"
91       "add x12, x12, #0x4\n"
92       "cmp x12, x9, LSL #1\n"
93       "blt 1b\n"
94       "2:"  // K loop: Charge: End
95       ".inst 0x25286580  // dup p0.h, p9.h/Z, p12.h[w12]\n"
96       ".inst 0xe05c02e0  // ld1h { za0h.h[x12] }, p0/Z, [x23, x28, LSL #1]\n"
97       ".inst 0x25286160  // dup p0.h, p8.h/Z, p11.h[w12]\n"
98       ".inst 0xe05c02c8  // ld1h { za1h.h[x12] }, p0/Z, [x22, x28, LSL #1]\n"
99       ".inst 0x25686580  // dup p0.h, p9.h/Z, p12.h[w12, #2]\n"
100       ".inst 0xe05c02a2  // ld1h { za0h.h[x12, #2] }, p0/Z, [x21, x28, LSL #1]\n"
101       ".inst 0x25686160  // dup p0.h, p8.h/Z, p11.h[w12, #2]\n"
102       ".inst 0xe05c028a  // ld1h { za1h.h[x12, #2] }, p0/Z, [x20, x28, LSL #1]\n"
103       "mov x25, %x[in]\n"
104       "add x24, %x[in], x16, LSL #3\n"
105       "ldr x23, [x25, #0x0]\n"
106       "ldr x22, [x24, #0x0]\n"
107       "ldr x21, [x25, #0x8]\n"
108       "ldr x20, [x24, #0x8]\n"
109       "add x25, x25, #0x10\n"
110       "add x24, x24, #0x10\n"
111       "inch x28\n"
112       "inch x26\n"
113       "cbz x13, 8f\n"
114       "mov x19, x13\n"
115       "3:"  // K loop: Main loop
116       "whilelt p9.h, x26, %x[width]\n"
117       "whilelt p8.h, x26, %x[width]\n"
118       "mov x13, #0x0\n"
119       "mov x12, #0x0\n"
120       "cbz x9, 5f\n"
121       "4:"  // K loop: Main loop: First: Loop
122       ".inst 0x25396580  // dup p0.h, p9.h/Z, p12.h[w13, #1]\n"
123       ".inst 0xe05c22e1  // ld1h { za0h.h[x13, #1] }, p0/Z, [x23, x28, LSL #1]\n"
124       ".inst 0x25396160  // dup p0.h, p8.h/Z, p11.h[w13, #1]\n"
125       ".inst 0xe05c22c9  // ld1h { za1h.h[x13, #1] }, p0/Z, [x22, x28, LSL #1]\n"
126       ".inst 0x25796580  // dup p0.h, p9.h/Z, p12.h[w13, #3]\n"
127       ".inst 0xe05c22a3  // ld1h { za0h.h[x13, #3] }, p0/Z, [x21, x28, LSL #1]\n"
128       ".inst 0x25796160  // dup p0.h, p8.h/Z, p11.h[w13, #3]\n"
129       ".inst 0xe05c228b  // ld1h { za1h.h[x13, #3] }, p0/Z, [x20, x28, LSL #1]\n"
130       "ldr x23, [x25, #0x0]\n"
131       "ldr x22, [x24, #0x0]\n"
132       "ldr x21, [x25, #0x8]\n"
133       "ldr x20, [x24, #0x8]\n"
134       "add x25, x25, #0x10\n"
135       "add x24, x24, #0x10\n"
136       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
137       ".inst 0xe0bf8360  // st1w { za0v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
138       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
139       ".inst 0x25707541  // dup p1.s, p13.s/Z, p10.s[w12, #1]\n"
140       ".inst 0xe0b08364  // st1w { za1v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
141       ".inst 0x25707540  // dup p0.s, p13.s/Z, p10.s[w12, #1]\n"
142       "add x13, x13, #0x4\n"
143       ".inst 0xe0af8761  // st1w { za0v.s[x12, #1] }, p1/Z, [x27, x15, LSL #2]\n"
144       ".inst 0xe0ae8365  // st1w { za1v.s[x12, #1] }, p0/Z, [x27, x14, LSL #2]\n"
145       "addvl x27, x27, #4\n"
146       "add x12, x12, #0x2\n"
147       "cmp x12, x9\n"
148       "blt 4b\n"
149       "5:"  // K loop: Main loop: First: Tail
150       "mov x25, %x[in]\n"
151       "add x24, %x[in], x16, LSL #3\n"
152       ".inst 0x25396580  // dup p0.h, p9.h/Z, p12.h[w13, #1]\n"
153       ".inst 0xe05c22e1  // ld1h { za0h.h[x13, #1] }, p0/Z, [x23, x28, LSL #1]\n"
154       ".inst 0x25396160  // dup p0.h, p8.h/Z, p11.h[w13, #1]\n"
155       ".inst 0xe05c22c9  // ld1h { za1h.h[x13, #1] }, p0/Z, [x22, x28, LSL #1]\n"
156       ".inst 0x25796580  // dup p0.h, p9.h/Z, p12.h[w13, #3]\n"
157       ".inst 0xe05c22a3  // ld1h { za0h.h[x13, #3] }, p0/Z, [x21, x28, LSL #1]\n"
158       ".inst 0x25796160  // dup p0.h, p8.h/Z, p11.h[w13, #3]\n"
159       ".inst 0xe05c228b  // ld1h { za1h.h[x13, #3] }, p0/Z, [x20, x28, LSL #1]\n"
160       "ldr x23, [x25, #0x0]\n"
161       "ldr x22, [x24, #0x0]\n"
162       "ldr x21, [x25, #0x8]\n"
163       "ldr x20, [x24, #0x8]\n"
164       "add x25, x25, #0x10\n"
165       "add x24, x24, #0x10\n"
166       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
167       ".inst 0xe0bf8360  // st1w { za0v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
168       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
169       ".inst 0x25707541  // dup p1.s, p13.s/Z, p10.s[w12, #1]\n"
170       ".inst 0xe0b08364  // st1w { za1v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
171       ".inst 0x25707540  // dup p0.s, p13.s/Z, p10.s[w12, #1]\n"
172       "whilelt p10.h, x26, %x[width]\n"
173       ".inst 0xe0af8761  // st1w { za0v.s[x12, #1] }, p1/Z, [x27, x15, LSL #2]\n"
174       "inch x26\n"
175       "inch x28\n"
176       ".inst 0xe0ae8365  // st1w { za1v.s[x12, #1] }, p0/Z, [x27, x14, LSL #2]\n"
177       "addvl x27, x27, #4\n"
178       "whilelt p9.h, x26, %x[width]\n"
179       "whilelt p8.h, x26, %x[width]\n"
180       "mov x13, #0x0\n"
181       "mov x12, #0x0\n"
182       "cbz x9, 7f\n"
183       "6:"  // K loop: Main loop: Second: Loop
184       ".inst 0x25296580  // dup p0.h, p9.h/Z, p12.h[w13]\n"
185       ".inst 0xe05c22e0  // ld1h { za0h.h[x13] }, p0/Z, [x23, x28, LSL #1]\n"
186       ".inst 0x25296160  // dup p0.h, p8.h/Z, p11.h[w13]\n"
187       ".inst 0xe05c22c8  // ld1h { za1h.h[x13] }, p0/Z, [x22, x28, LSL #1]\n"
188       ".inst 0x25696580  // dup p0.h, p9.h/Z, p12.h[w13, #2]\n"
189       ".inst 0xe05c22a2  // ld1h { za0h.h[x13, #2] }, p0/Z, [x21, x28, LSL #1]\n"
190       ".inst 0x25696160  // dup p0.h, p8.h/Z, p11.h[w13, #2]\n"
191       ".inst 0xe05c228a  // ld1h { za1h.h[x13, #2] }, p0/Z, [x20, x28, LSL #1]\n"
192       "ldr x23, [x25, #0x0]\n"
193       "ldr x22, [x24, #0x0]\n"
194       "ldr x21, [x25, #0x8]\n"
195       "ldr x20, [x24, #0x8]\n"
196       "add x25, x25, #0x10\n"
197       "add x24, x24, #0x10\n"
198       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
199       ".inst 0xe0bf8368  // st1w { za2v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
200       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
201       ".inst 0x25707541  // dup p1.s, p13.s/Z, p10.s[w12, #1]\n"
202       ".inst 0xe0b0836c  // st1w { za3v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
203       ".inst 0x25707540  // dup p0.s, p13.s/Z, p10.s[w12, #1]\n"
204       "add x13, x13, #0x4\n"
205       ".inst 0xe0af8769  // st1w { za2v.s[x12, #1] }, p1/Z, [x27, x15, LSL #2]\n"
206       ".inst 0xe0ae836d  // st1w { za3v.s[x12, #1] }, p0/Z, [x27, x14, LSL #2]\n"
207       "addvl x27, x27, #4\n"
208       "add x12, x12, #0x2\n"
209       "cmp x12, x9\n"
210       "blt 6b\n"
211       "7:"  // K loop: Main loop: Second: Tail
212       "mov x25, %x[in]\n"
213       "add x24, %x[in], x16, LSL #3\n"
214       ".inst 0x25296580  // dup p0.h, p9.h/Z, p12.h[w13]\n"
215       ".inst 0xe05c22e0  // ld1h { za0h.h[x13] }, p0/Z, [x23, x28, LSL #1]\n"
216       ".inst 0x25296160  // dup p0.h, p8.h/Z, p11.h[w13]\n"
217       ".inst 0xe05c22c8  // ld1h { za1h.h[x13] }, p0/Z, [x22, x28, LSL #1]\n"
218       ".inst 0x25696580  // dup p0.h, p9.h/Z, p12.h[w13, #2]\n"
219       ".inst 0xe05c22a2  // ld1h { za0h.h[x13, #2] }, p0/Z, [x21, x28, LSL #1]\n"
220       ".inst 0x25696160  // dup p0.h, p8.h/Z, p11.h[w13, #2]\n"
221       ".inst 0xe05c228a  // ld1h { za1h.h[x13, #2] }, p0/Z, [x20, x28, LSL #1]\n"
222       "ldr x23, [x25, #0x0]\n"
223       "ldr x22, [x24, #0x0]\n"
224       "ldr x21, [x25, #0x8]\n"
225       "ldr x20, [x24, #0x8]\n"
226       "add x25, x25, #0x10\n"
227       "add x24, x24, #0x10\n"
228       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
229       ".inst 0xe0bf8368  // st1w { za2v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
230       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
231       ".inst 0x25707541  // dup p1.s, p13.s/Z, p10.s[w12, #1]\n"
232       ".inst 0xe0b0836c  // st1w { za3v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
233       ".inst 0x25707540  // dup p0.s, p13.s/Z, p10.s[w12, #1]\n"
234       "whilelt p10.h, x26, %x[width]\n"
235       ".inst 0xe0af8769  // st1w { za2v.s[x12, #1] }, p1/Z, [x27, x15, LSL #2]\n"
236       "inch x26\n"
237       "inch x28\n"
238       ".inst 0xe0ae836d  // st1w { za3v.s[x12, #1] }, p0/Z, [x27, x14, LSL #2]\n"
239       "addvl x27, x27, #4\n"
240       "subs x19, x19, #0x1\n"
241       "bgt 3b\n"
242       "8:"  // K loop: Tails
243       "cbnz x11, 11f\n"
244       "mov x25, %x[in]\n"
245       "whilelt p9.h, x26, %x[width]\n"
246       "whilelt p8.h, x26, %x[width]\n"
247       "mov x13, #0x0\n"
248       "mov x12, #0x0\n"
249       "9:"  // K loop: Tails: Even: First
250       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
251       ".inst 0xe0bf8360  // st1w { za0v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
252       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
253       ".inst 0x25396581  // dup p1.h, p9.h/Z, p12.h[w13, #1]\n"
254       ".inst 0xe0b08364  // st1w { za1v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
255       ".inst 0x25396160  // dup p0.h, p8.h/Z, p11.h[w13, #1]\n"
256       "addvl x27, x27, #2\n"
257       "ldr x23, [x25, #0x0]\n"
258       ".inst 0xe05c26e1  // ld1h { za0h.h[x13, #1] }, p1/Z, [x23, x28, LSL #1]\n"
259       "ldr x22, [x25, x16, LSL #0x3]\n"
260       ".inst 0xe05c22c9  // ld1h { za1h.h[x13, #1] }, p0/Z, [x22, x28, LSL #1]\n"
261       "add x12, x12, #0x1\n"
262       "add x25, x25, #0x8\n"
263       "add x13, x13, #0x2\n"
264       "cmp x12, x16\n"
265       "blt 9b\n"
266       "whilelt p10.h, x26, %x[width]\n"
267       "whilelt p9.h, x26, %x[width]\n"
268       "whilelt p8.h, x26, %x[width]\n"
269       "mov x19, #0x0\n"
270       "mov x12, #0x0\n"
271       "10:"  // K loop: Tails: Even: Second
272       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
273       ".inst 0xe0bf8368  // st1w { za2v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
274       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
275       "add x19, x19, #0x2\n"
276       ".inst 0xe0b0836c  // st1w { za3v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
277       "addvl x27, x27, #2\n"
278       "add x12, x12, #0x1\n"
279       "cmp x12, x10\n"
280       "blt 10b\n"
281       "whilelt p10.h, x26, %x[width]\n"
282       "b 13f\n"
283       "11:"  // K loop: Tails: Odd
284       "mov x12, #0x0\n"
285       "12:"  // K loop: Tails: Odd: Loop
286       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
287       ".inst 0xe0bf8360  // st1w { za0v.s[x12] }, p0/Z, [x27, XZR, LSL #2]\n"
288       ".inst 0x25307540  // dup p0.s, p13.s/Z, p10.s[w12]\n"
289       ".inst 0xe0b08364  // st1w { za1v.s[x12] }, p0/Z, [x27, x16, LSL #2]\n"
290       "addvl x27, x27, #2\n"
291       "add x12, x12, #0x1\n"
292       "cmp x12, x10\n"
293       "blt 12b\n"
294       "13:"  // K loop: End
295       "mov %x[out], x27\n"
296       ".inst 0xd503467f  // SMSTOP\n"
297       : [out] "+&r" (out)
298       : [height] "r" (height), [in] "r" (in), [row_offset] "r" (row_offset), [width] "r" (width)
299       : "cc", "memory", "p0", "p1", "p8", "p9", "p10", "p11", "p12", "p13", "x9", "x10", "x11", "x12", "x13", "x14", "x15", "x16", "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"
300     );
301 }
302 
303 #endif  // defined(__ARM_FEATURE_SVE)
304