1 /*
2  * Copyright (c) 2021-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 all
14  * 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 FROM,
21  * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
22  * SOFTWARE.
23  */
24 
25 #include <cstdint>
26 #include <cstddef>
27 
28 #if defined(ARM_COMPUTE_ENABLE_SVE) && defined(__ARM_FP16_ARGS)
29 
30 namespace arm_conv {
31 namespace pooling {
32 
33 
sve_fp16_nhwc_max_generic_depthfirst_impl(const uint64_t,const uint64_t n_valid_cells,uint64_t n_channels,const __fp16 * const * const inptrs,__fp16 * outptr)34 void sve_fp16_nhwc_max_generic_depthfirst_impl(
35   const uint64_t,
36   const uint64_t n_valid_cells,
37   uint64_t n_channels,
38   const __fp16 *const *const inptrs,
39   __fp16 *outptr
40 )
41 {
42   __asm__ __volatile__(
43     "mov x28, #0x0\n"
44     "cnth x27\n"
45     "cnth x26, ALL, MUL #2\n"
46     "cnth x25, ALL, MUL #3\n"
47     "whilelt p4.h, x28, %x[n_channels]\n"
48     "whilelt p3.h, x27, %x[n_channels]\n"
49     "whilelt p2.h, x26, %x[n_channels]\n"
50     "whilelt p1.h, x25, %x[n_channels]\n"
51     "ptrue p0.b\n"
52     "b.none 7f\n"
53     "1:"  // 4-vectors of channels
54     "lsr x24, %x[n_valid_cells], #0x2\n"
55     "mov z8.h, #0xfc00\n"
56     "mov z7.h, #0xfc00\n"
57     "mov x19, %x[inptrs]\n"
58     "mov z6.h, #0xfc00\n"
59     "mov z5.h, #0xfc00\n"
60     "cbz x24, 4f\n"
61     "ldp x23, x22, [x19, #0x0]\n"
62     "ldp x21, x20, [x19, #0x10]\n"
63     "subs x24, x24, #0x1\n"
64     "add x19, x19, #0x20\n"
65     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
66     "ld1h { z3.h }, p4/Z, [x22, x28, LSL #1]\n"
67     "ld1h { z2.h }, p4/Z, [x21, x28, LSL #1]\n"
68     "ld1h { z1.h }, p4/Z, [x20, x28, LSL #1]\n"
69     "ld1h { z0.h }, p3/Z, [x23, x27, LSL #1]\n"
70     "ld1h { z31.h }, p3/Z, [x22, x27, LSL #1]\n"
71     "ld1h { z22.h }, p3/Z, [x21, x27, LSL #1]\n"
72     "ld1h { z30.h }, p3/Z, [x20, x27, LSL #1]\n"
73     "ld1h { z29.h }, p2/Z, [x23, x26, LSL #1]\n"
74     "ld1h { z28.h }, p2/Z, [x22, x26, LSL #1]\n"
75     "ld1h { z21.h }, p2/Z, [x21, x26, LSL #1]\n"
76     "ld1h { z27.h }, p2/Z, [x20, x26, LSL #1]\n"
77     "ld1h { z26.h }, p1/Z, [x23, x25, LSL #1]\n"
78     "ld1h { z25.h }, p1/Z, [x22, x25, LSL #1]\n"
79     "ld1h { z20.h }, p1/Z, [x21, x25, LSL #1]\n"
80     "ld1h { z24.h }, p1/Z, [x20, x25, LSL #1]\n"
81     "beq 3f\n"
82     "2:"  // 4-vectors of channels: 4 inputs loop
83     "movprfx z19, z4\n fmax z19.h, p0/M, z19.h, z3.h\n"
84     "movprfx z23, z2\n fmax z23.h, p0/M, z23.h, z1.h\n"
85     "ldp x23, x22, [x19, #0x0]\n"
86     "ldp x21, x20, [x19, #0x10]\n"
87     "movprfx z18, z0\n fmax z18.h, p0/M, z18.h, z31.h\n"
88     "fmax z22.h, p0/M, z22.h, z30.h\n"
89     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
90     "ld1h { z3.h }, p4/Z, [x22, x28, LSL #1]\n"
91     "movprfx z17, z29\n fmax z17.h, p0/M, z17.h, z28.h\n"
92     "fmax z21.h, p0/M, z21.h, z27.h\n"
93     "ld1h { z2.h }, p4/Z, [x21, x28, LSL #1]\n"
94     "ld1h { z1.h }, p4/Z, [x20, x28, LSL #1]\n"
95     "movprfx z16, z26\n fmax z16.h, p0/M, z16.h, z25.h\n"
96     "fmax z20.h, p0/M, z20.h, z24.h\n"
97     "ld1h { z0.h }, p3/Z, [x23, x27, LSL #1]\n"
98     "ld1h { z31.h }, p3/Z, [x22, x27, LSL #1]\n"
99     "fmax z19.h, p0/M, z19.h, z23.h\n"
100     "fmax z18.h, p0/M, z18.h, z22.h\n"
101     "ld1h { z22.h }, p3/Z, [x21, x27, LSL #1]\n"
102     "ld1h { z30.h }, p3/Z, [x20, x27, LSL #1]\n"
103     "fmax z17.h, p0/M, z17.h, z21.h\n"
104     "fmax z16.h, p0/M, z16.h, z20.h\n"
105     "ld1h { z29.h }, p2/Z, [x23, x26, LSL #1]\n"
106     "ld1h { z28.h }, p2/Z, [x22, x26, LSL #1]\n"
107     "subs x24, x24, #0x1\n"
108     "fmax z8.h, p0/M, z8.h, z19.h\n"
109     "ld1h { z21.h }, p2/Z, [x21, x26, LSL #1]\n"
110     "ld1h { z27.h }, p2/Z, [x20, x26, LSL #1]\n"
111     "fmax z7.h, p0/M, z7.h, z18.h\n"
112     "fmax z6.h, p0/M, z6.h, z17.h\n"
113     "ld1h { z26.h }, p1/Z, [x23, x25, LSL #1]\n"
114     "ld1h { z25.h }, p1/Z, [x22, x25, LSL #1]\n"
115     "fmax z5.h, p0/M, z5.h, z16.h\n"
116     "add x19, x19, #0x20\n"
117     "ld1h { z20.h }, p1/Z, [x21, x25, LSL #1]\n"
118     "ld1h { z24.h }, p1/Z, [x20, x25, LSL #1]\n"
119     "bgt 2b\n"
120     "3:"  // 4-vectors of channels: 4 inputs tail
121     "movprfx z19, z4\n fmax z19.h, p0/M, z19.h, z3.h\n"
122     "movprfx z23, z2\n fmax z23.h, p0/M, z23.h, z1.h\n"
123     "movprfx z18, z0\n fmax z18.h, p0/M, z18.h, z31.h\n"
124     "fmax z22.h, p0/M, z22.h, z30.h\n"
125     "movprfx z17, z29\n fmax z17.h, p0/M, z17.h, z28.h\n"
126     "fmax z21.h, p0/M, z21.h, z27.h\n"
127     "movprfx z16, z26\n fmax z16.h, p0/M, z16.h, z25.h\n"
128     "fmax z20.h, p0/M, z20.h, z24.h\n"
129     "fmax z19.h, p0/M, z19.h, z23.h\n"
130     "fmax z18.h, p0/M, z18.h, z22.h\n"
131     "fmax z17.h, p0/M, z17.h, z21.h\n"
132     "fmax z16.h, p0/M, z16.h, z20.h\n"
133     "fmax z8.h, p0/M, z8.h, z19.h\n"
134     "fmax z7.h, p0/M, z7.h, z18.h\n"
135     "fmax z6.h, p0/M, z6.h, z17.h\n"
136     "fmax z5.h, p0/M, z5.h, z16.h\n"
137     "4:"  // 4-vectors of channels: After loop
138     "ands x20, %x[n_valid_cells], #0x3\n"
139     "beq 6f\n"
140     "5:"  // 4-vectors of channels: Single input loop
141     "ldr x23, [x19], #0x8\n"
142     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
143     "subs x20, x20, #0x1\n"
144     "fmax z8.h, p0/M, z8.h, z4.h\n"
145     "ld1h { z0.h }, p3/Z, [x23, x27, LSL #1]\n"
146     "ld1h { z29.h }, p2/Z, [x23, x26, LSL #1]\n"
147     "fmax z7.h, p0/M, z7.h, z0.h\n"
148     "fmax z6.h, p0/M, z6.h, z29.h\n"
149     "ld1h { z26.h }, p1/Z, [x23, x25, LSL #1]\n"
150     "fmax z5.h, p0/M, z5.h, z26.h\n"
151     "bgt 5b\n"
152     "6:"  // 4-vectors of channels: Single input loop: End
153     "st1h { z8.h }, p4, [%x[outptr], x28, LSL #1]\n"
154     "inch x28, ALL, MUL #4\n"
155     "st1h { z7.h }, p3, [%x[outptr], x27, LSL #1]\n"
156     "inch x27, ALL, MUL #4\n"
157     "st1h { z6.h }, p2, [%x[outptr], x26, LSL #1]\n"
158     "inch x26, ALL, MUL #4\n"
159     "st1h { z5.h }, p1, [%x[outptr], x25, LSL #1]\n"
160     "inch x25, ALL, MUL #4\n"
161     "whilelt p1.h, x25, %x[n_channels]\n"
162     "b.any 1b\n"
163     "7:"  // Single vector of channels
164     "whilelt p4.h, x28, %x[n_channels]\n"
165     "b.none 14f\n"
166     "8:"  // Single vector of channels: Loop
167     "lsr x24, %x[n_valid_cells], #0x2\n"
168     "mov z8.h, #0xfc00\n"
169     "mov x19, %x[inptrs]\n"
170     "cbz x24, 11f\n"
171     "ldp x23, x22, [x19, #0x0]\n"
172     "ldp x21, x20, [x19, #0x10]\n"
173     "subs x24, x24, #0x1\n"
174     "add x19, x19, #0x20\n"
175     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
176     "ld1h { z3.h }, p4/Z, [x22, x28, LSL #1]\n"
177     "ld1h { z2.h }, p4/Z, [x21, x28, LSL #1]\n"
178     "ld1h { z1.h }, p4/Z, [x20, x28, LSL #1]\n"
179     "beq 10f\n"
180     "9:"  // Single vector of channels: Loop: 4 inputs loop
181     "movprfx z19, z4\n fmax z19.h, p0/M, z19.h, z3.h\n"
182     "movprfx z23, z2\n fmax z23.h, p0/M, z23.h, z1.h\n"
183     "ldp x23, x22, [x19, #0x0]\n"
184     "ldp x21, x20, [x19, #0x10]\n"
185     "fmax z19.h, p0/M, z19.h, z23.h\n"
186     "subs x24, x24, #0x1\n"
187     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
188     "ld1h { z3.h }, p4/Z, [x22, x28, LSL #1]\n"
189     "fmax z8.h, p0/M, z8.h, z19.h\n"
190     "add x19, x19, #0x20\n"
191     "ld1h { z2.h }, p4/Z, [x21, x28, LSL #1]\n"
192     "ld1h { z1.h }, p4/Z, [x20, x28, LSL #1]\n"
193     "bgt 9b\n"
194     "10:"  // Single vector of channels: Loop: 4 inputs tail
195     "movprfx z19, z4\n fmax z19.h, p0/M, z19.h, z3.h\n"
196     "movprfx z23, z2\n fmax z23.h, p0/M, z23.h, z1.h\n"
197     "fmax z19.h, p0/M, z19.h, z23.h\n"
198     "fmax z8.h, p0/M, z8.h, z19.h\n"
199     "11:"  // Single vector of channels: Loop: After loop
200     "ands x20, %x[n_valid_cells], #0x3\n"
201     "beq 13f\n"
202     "12:"  // Single vector of channels: Loop: Single input loop
203     "ldr x23, [x19], #0x8\n"
204     "ld1h { z4.h }, p4/Z, [x23, x28, LSL #1]\n"
205     "subs x20, x20, #0x1\n"
206     "fmax z8.h, p0/M, z8.h, z4.h\n"
207     "bgt 12b\n"
208     "13:"  // Single vector of channels: Loop: Single input loop: End
209     "st1h { z8.h }, p4, [%x[outptr], x28, LSL #1]\n"
210     "inch x28\n"
211     "whilelt p4.h, x28, %x[n_channels]\n"
212     "b.any 8b\n"
213     "14:"  // End
214     :
215     : [inptrs] "r" (inptrs), [n_channels] "r" (n_channels), [n_valid_cells] "r" (n_valid_cells), [outptr] "r" (outptr)
216     : "cc", "memory", "p0", "p1", "p2", "p3", "p4", "x19", "x20", "x21", "x22", "x23", "x24", "x25", "x26", "x27", "x28", "z0", "z1", "z2", "z3", "z4", "z5", "z6", "z7", "z8", "z16", "z17", "z18", "z19", "z20", "z21", "z22", "z23", "z24", "z25", "z26", "z27", "z28", "z29", "z30", "z31"
217   );
218 }
219 
220 }  // namespace pooling
221 }  // namespace arm_conv
222 
223 #endif  // defined(ARM_COMPUTE_ENABLE_SVE) && defined(__ARM_FP16_ARGS)
224