1namespace dnnl {
2namespace impl {
3namespace gpu {
4namespace ocl {
5const char *gen9_gemm_nocopy_f16_kernel = R"==(/******************************************************************************* )==""\n"
6R"==(* Copyright 2019-2022 Intel Corporation )==""\n"
7R"==(* )==""\n"
8R"==(* Licensed under the Apache License, Version 2.0 (the "License"); )==""\n"
9R"==(* you may not use this file except in compliance with the License. )==""\n"
10R"==(* You may obtain a copy of the License at )==""\n"
11R"==(* )==""\n"
12R"==(* http: )==""\n"
13R"==(* )==""\n"
14R"==(* Unless required by applicable law or agreed to in writing, software )==""\n"
15R"==(* distributed under the License is distributed on an "AS IS" BASIS, )==""\n"
16R"==(* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. )==""\n"
17R"==(* See the License for the specific language governing permissions and )==""\n"
18R"==(* limitations under the License. )==""\n"
19R"==(*******************************************************************************/ )==""\n"
20R"==(#include "gpu/ocl/gemm/ocl_gemm_attrs.h" )==""\n"
21R"==(#include "gpu/ocl/ocl_post_ops.h" )==""\n"
22R"==(#include "gpu/ocl/ocl_types.h" )==""\n"
23R"==(#pragma OPENCL EXTENSION cl_khr_fp16 : enable )==""\n"
24R"==(#define VLOAD4_ALIGNED(o, p) *((global half4 *)p + o) )==""\n"
25R"==(#if !defined(TRANS_A) )==""\n"
26R"==(#if !defined(TRANS_B) )==""\n"
27R"==(#define NN )==""\n"
28R"==(#define DO_FMA DO_FMA_NN )==""\n"
29R"==(#else )==""\n"
30R"==(#define NT )==""\n"
31R"==(#define DO_FMA DO_FMA_NT )==""\n"
32R"==(#endif )==""\n"
33R"==(#else )==""\n"
34R"==(#if !defined(TRANS_B) )==""\n"
35R"==(#define TN )==""\n"
36R"==(#define DO_FMA DO_FMA_TN )==""\n"
37R"==(#else )==""\n"
38R"==(#define TT )==""\n"
39R"==(#define DO_FMA DO_FMA_TT )==""\n"
40R"==(#endif )==""\n"
41R"==(#endif )==""\n"
42R"==(#define DO_FMA_NN(h, i, i_mod_2, i_div_2, i_mod_16, i_div_16) \ )==""\n"
43R"==(do { \ )==""\n"
44R"==(c[i][0] = mad(sub_group_broadcast(a[h].s##i_mod_2, i_div_2), \ )==""\n"
45R"==(b[0].s##h, c[i][0]); \ )==""\n"
46R"==(c[i][1] = mad(sub_group_broadcast(a[h].s##i_mod_2, i_div_2), \ )==""\n"
47R"==(b[1].s##h, c[i][1]); \ )==""\n"
48R"==(} while (0) )==""\n"
49R"==(#define DO_FMA_NT(h, i, i_mod_2, i_div_2, i_mod_16, i_div_16) \ )==""\n"
50R"==(do { \ )==""\n"
51R"==(c[i][0] = mad(sub_group_broadcast(a[h].s##i_mod_2, i_div_2), b[h].s0, \ )==""\n"
52R"==(c[i][0]); \ )==""\n"
53R"==(c[i][1] = mad(sub_group_broadcast(a[h].s##i_mod_2, i_div_2), b[h].s1, \ )==""\n"
54R"==(c[i][1]); \ )==""\n"
55R"==(} while (0) )==""\n"
56R"==(#define DO_FMA_TN(h, i, i_mod_2, i_div_2, i_mod_16, i_div_16) \ )==""\n"
57R"==(do { \ )==""\n"
58R"==(c[i][0] = mad(sub_group_broadcast(a[i_div_16].s##h, i_mod_16), \ )==""\n"
59R"==(b[0].s##h, c[i][0]); \ )==""\n"
60R"==(c[i][1] = mad(sub_group_broadcast(a[i_div_16].s##h, i_mod_16), \ )==""\n"
61R"==(b[1].s##h, c[i][1]); \ )==""\n"
62R"==(} while (0) )==""\n"
63R"==(#define DO_FMA_TT(h, i, i_mod_2, i_div_2, i_mod_16, i_div_16) \ )==""\n"
64R"==(do { \ )==""\n"
65R"==(c[i][0] = mad(sub_group_broadcast(a[i_div_16].s##h, i_mod_16), \ )==""\n"
66R"==(b[h].s0, c[i][0]); \ )==""\n"
67R"==(c[i][1] = mad(sub_group_broadcast(a[i_div_16].s##h, i_mod_16), \ )==""\n"
68R"==(b[h].s1, c[i][1]); \ )==""\n"
69R"==(} while (0) )==""\n"
70R"==(#if WITH_ELTWISE == 1 )==""\n"
71R"==(#define POST_OP(val) \ )==""\n"
72R"==(do { \ )==""\n"
73R"==(if (last_k_block && last_k_unroll) \ )==""\n"
74R"==(val = fwd_eltwise( \ )==""\n"
75R"==(val, eltwise_alpha, eltwise_beta, eltwise_scale); \ )==""\n"
76R"==(} while (0) )==""\n"
77R"==(#else )==""\n"
78R"==(#define POST_OP(val) )==""\n"
79R"==(#endif )==""\n"
80R"==(#define FMA_I_LOOP(h) \ )==""\n"
81R"==(do { \ )==""\n"
82R"==(DO_FMA(h, 0, 0, 0, 0, 0); \ )==""\n"
83R"==(DO_FMA(h, 1, 1, 0, 1, 0); \ )==""\n"
84R"==(DO_FMA(h, 2, 0, 1, 2, 0); \ )==""\n"
85R"==(DO_FMA(h, 3, 1, 1, 3, 0); \ )==""\n"
86R"==(DO_FMA(h, 4, 0, 2, 4, 0); \ )==""\n"
87R"==(DO_FMA(h, 5, 1, 2, 5, 0); \ )==""\n"
88R"==(DO_FMA(h, 6, 0, 3, 6, 0); \ )==""\n"
89R"==(DO_FMA(h, 7, 1, 3, 7, 0); \ )==""\n"
90R"==(DO_FMA(h, 8, 0, 4, 8, 0); \ )==""\n"
91R"==(DO_FMA(h, 9, 1, 4, 9, 0); \ )==""\n"
92R"==(DO_FMA(h, 10, 0, 5, 10, 0); \ )==""\n"
93R"==(DO_FMA(h, 11, 1, 5, 11, 0); \ )==""\n"
94R"==(DO_FMA(h, 12, 0, 6, 12, 0); \ )==""\n"
95R"==(DO_FMA(h, 13, 1, 6, 13, 0); \ )==""\n"
96R"==(DO_FMA(h, 14, 0, 7, 14, 0); \ )==""\n"
97R"==(DO_FMA(h, 15, 1, 7, 15, 0); \ )==""\n"
98R"==(DO_FMA(h, 16, 0, 8, 0, 1); \ )==""\n"
99R"==(DO_FMA(h, 17, 1, 8, 1, 1); \ )==""\n"
100R"==(DO_FMA(h, 18, 0, 9, 2, 1); \ )==""\n"
101R"==(DO_FMA(h, 19, 1, 9, 3, 1); \ )==""\n"
102R"==(DO_FMA(h, 20, 0, 10, 4, 1); \ )==""\n"
103R"==(DO_FMA(h, 21, 1, 10, 5, 1); \ )==""\n"
104R"==(DO_FMA(h, 22, 0, 11, 6, 1); \ )==""\n"
105R"==(DO_FMA(h, 23, 1, 11, 7, 1); \ )==""\n"
106R"==(DO_FMA(h, 24, 0, 12, 8, 1); \ )==""\n"
107R"==(DO_FMA(h, 25, 1, 12, 9, 1); \ )==""\n"
108R"==(DO_FMA(h, 26, 0, 13, 10, 1); \ )==""\n"
109R"==(DO_FMA(h, 27, 1, 13, 11, 1); \ )==""\n"
110R"==(DO_FMA(h, 28, 0, 14, 12, 1); \ )==""\n"
111R"==(DO_FMA(h, 29, 1, 14, 13, 1); \ )==""\n"
112R"==(DO_FMA(h, 30, 0, 15, 14, 1); \ )==""\n"
113R"==(DO_FMA(h, 31, 1, 15, 15, 1); \ )==""\n"
114R"==(} while (0) )==""\n"
115R"==(#define UPDATE_C_ROW(i, betaZero) \ )==""\n"
116R"==(do { \ )==""\n"
117R"==(if (irem > i) { \ )==""\n"
118R"==(if (jrem > 0) { \ )==""\n"
119R"==(half val \ )==""\n"
120R"==(= ATTR_ALPHA * c[i][0] + ((betaZero) ? 0 : beta * *C); \ )==""\n"
121R"==(POST_OP(val); \ )==""\n"
122R"==(*C = val; \ )==""\n"
123R"==(} \ )==""\n"
124R"==(if (jrem > 16) { \ )==""\n"
125R"==(half val = ATTR_ALPHA * c[i][1] \ )==""\n"
126R"==(+ ((betaZero) ? 0 : beta * *C2); \ )==""\n"
127R"==(POST_OP(val); \ )==""\n"
128R"==(*C2 = val; \ )==""\n"
129R"==(} \ )==""\n"
130R"==(} \ )==""\n"
131R"==(C++; \ )==""\n"
132R"==(C2++; \ )==""\n"
133R"==(} while (0) )==""\n"
134R"==(#define UPDATE_C(betaZero) \ )==""\n"
135R"==(do { \ )==""\n"
136R"==(UPDATE_C_ROW(0, betaZero); \ )==""\n"
137R"==(UPDATE_C_ROW(1, betaZero); \ )==""\n"
138R"==(UPDATE_C_ROW(2, betaZero); \ )==""\n"
139R"==(UPDATE_C_ROW(3, betaZero); \ )==""\n"
140R"==(UPDATE_C_ROW(4, betaZero); \ )==""\n"
141R"==(UPDATE_C_ROW(5, betaZero); \ )==""\n"
142R"==(UPDATE_C_ROW(6, betaZero); \ )==""\n"
143R"==(UPDATE_C_ROW(7, betaZero); \ )==""\n"
144R"==(UPDATE_C_ROW(8, betaZero); \ )==""\n"
145R"==(UPDATE_C_ROW(9, betaZero); \ )==""\n"
146R"==(UPDATE_C_ROW(10, betaZero); \ )==""\n"
147R"==(UPDATE_C_ROW(11, betaZero); \ )==""\n"
148R"==(UPDATE_C_ROW(12, betaZero); \ )==""\n"
149R"==(UPDATE_C_ROW(13, betaZero); \ )==""\n"
150R"==(UPDATE_C_ROW(14, betaZero); \ )==""\n"
151R"==(UPDATE_C_ROW(15, betaZero); \ )==""\n"
152R"==(UPDATE_C_ROW(16, betaZero); \ )==""\n"
153R"==(UPDATE_C_ROW(17, betaZero); \ )==""\n"
154R"==(UPDATE_C_ROW(18, betaZero); \ )==""\n"
155R"==(UPDATE_C_ROW(19, betaZero); \ )==""\n"
156R"==(UPDATE_C_ROW(20, betaZero); \ )==""\n"
157R"==(UPDATE_C_ROW(21, betaZero); \ )==""\n"
158R"==(UPDATE_C_ROW(22, betaZero); \ )==""\n"
159R"==(UPDATE_C_ROW(23, betaZero); \ )==""\n"
160R"==(UPDATE_C_ROW(24, betaZero); \ )==""\n"
161R"==(UPDATE_C_ROW(25, betaZero); \ )==""\n"
162R"==(UPDATE_C_ROW(26, betaZero); \ )==""\n"
163R"==(UPDATE_C_ROW(27, betaZero); \ )==""\n"
164R"==(UPDATE_C_ROW(28, betaZero); \ )==""\n"
165R"==(UPDATE_C_ROW(29, betaZero); \ )==""\n"
166R"==(UPDATE_C_ROW(30, betaZero); \ )==""\n"
167R"==(UPDATE_C_ROW(31, betaZero); \ )==""\n"
168R"==(} while (0) )==""\n"
169R"==(#ifdef NN )==""\n"
170R"==(__attribute__((intel_reqd_sub_group_size(16))) )==""\n"
171R"==(kernel void )==""\n"
172R"==(gen9_gemm_nocopy_f16(global half *A, global half *B, global half *C, )==""\n"
173R"==(long offset_a, long offset_b, long offset_c, int lda, int ldb, int ldc, )==""\n"
174R"==(int m, int n, int k, global float *alpha, float beta, int last_k_block, )==""\n"
175R"==(float eltwise_alpha, float eltwise_beta, float eltwise_scale )==""\n"
176R"==(#ifdef WITH_K_UNROLL )==""\n"
177R"==(, )==""\n"
178R"==(volatile global int *flag, long offset_f) { )==""\n"
179R"==(#else )==""\n"
180R"==() { )==""\n"
181R"==(#endif )==""\n"
182R"==(half2 a[4]; )==""\n"
183R"==(half4 b[2]; )==""\n"
184R"==(half c[32][2]; )==""\n"
185R"==(int idM = get_global_id(1); )==""\n"
186R"==(int idN = get_global_id(0); )==""\n"
187R"==(int lid = get_sub_group_local_id(); )==""\n"
188R"==(int idK = get_global_id(2); )==""\n"
189R"==(int nku = get_global_size(2); )==""\n"
190R"==(int i0 = idM * 32; )==""\n"
191R"==(int j0 = sub_group_broadcast(idN, 0) * 2 + lid; )==""\n"
192R"==(int irem = m - i0; )==""\n"
193R"==(int jrem = n - j0; )==""\n"
194R"==(if (irem < 0) irem = 0; )==""\n"
195R"==(if (jrem < 0) jrem = 0; )==""\n"
196R"==(int irem2 = (irem + 1) >> 1; )==""\n"
197R"==(int last_k_unroll = (idK == nku - 1); )==""\n"
198R"==(#ifdef WITH_K_UNROLL )==""\n"
199R"==(int k0 = idK * UNROLL_K; )==""\n"
200R"==(int kt = k - k0; )==""\n"
201R"==(if (kt < 0) kt = 0; )==""\n"
202R"==(if (kt > UNROLL_K) kt = UNROLL_K; )==""\n"
203R"==(A += offset_a + i0 + k0 * lda; )==""\n"
204R"==(B += offset_b + j0 * ldb + k0; )==""\n"
205R"==(k = kt; )==""\n"
206R"==(#else )==""\n"
207R"==(A += offset_a + i0; )==""\n"
208R"==(B += offset_b + j0 * ldb; )==""\n"
209R"==(#endif )==""\n"
210R"==(C += offset_c + i0 + j0 * ldc; )==""\n"
211R"==(global half *A_ptrs[4] = {A, A + lda, A + 2 * lda, A + 3 * lda}; )==""\n"
212R"==(global half *B_ptrs[2] = {B, B + 16 * ldb}; )==""\n"
213R"==(for (int y = 0; y < 32; y++) )==""\n"
214R"==(for (int z = 0; z < 2; z++) )==""\n"
215R"==(c[y][z] = 0; )==""\n"
216R"==(#ifdef WITH_K_UNROLL )==""\n"
217R"==(flag += offset_f + idM; )==""\n"
218R"==(if (idK == 0 && lid == 0) *flag = 0; )==""\n"
219R"==(#endif )==""\n"
220R"==(int k_align = k & ~3; )==""\n"
221R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
222R"==(if (irem >= 32 && sub_group_broadcast(jrem, 0) >= 32) { )==""\n"
223R"==(#endif )==""\n"
224R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
225R"==(for (int hh = 0; hh < 4; hh++) )==""\n"
226R"==(a[hh] = as_half2(intel_sub_group_block_read( )==""\n"
227R"==((global uint *)(A_ptrs[hh] + h * lda))); )==""\n"
228R"==(for (int j = 0; j < 2; j++) )==""\n"
229R"==(b[j] = VLOAD4_ALIGNED(0, (B_ptrs[j] + h)); )==""\n"
230R"==(FMA_I_LOOP(0); )==""\n"
231R"==(FMA_I_LOOP(1); )==""\n"
232R"==(FMA_I_LOOP(2); )==""\n"
233R"==(FMA_I_LOOP(3); )==""\n"
234R"==(} )==""\n"
235R"==(for (int h = k_align; h < k; h++) { )==""\n"
236R"==(a[0] = as_half2(intel_sub_group_block_read( )==""\n"
237R"==((global uint *)(A_ptrs[0] + h * lda))); )==""\n"
238R"==(for (int j = 0; j < 2; j++) )==""\n"
239R"==(b[j] = B_ptrs[j][h]; )==""\n"
240R"==(FMA_I_LOOP(0); )==""\n"
241R"==(} )==""\n"
242R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
243R"==(} else { )==""\n"
244R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
245R"==(for (int hh = 0; hh < 4; hh++) { )==""\n"
246R"==(if (irem2 > lid) )==""\n"
247R"==(a[hh] = as_half2( )==""\n"
248R"==(*((global uint *)(A_ptrs[hh] + h * lda) + lid)); )==""\n"
249R"==(} )==""\n"
250R"==(for (int j = 0; j < 2; j++) )==""\n"
251R"==(if (jrem > j * 16) b[j] = VLOAD4_ALIGNED(0, (B_ptrs[j] + h)); )==""\n"
252R"==(FMA_I_LOOP(0); )==""\n"
253R"==(FMA_I_LOOP(1); )==""\n"
254R"==(FMA_I_LOOP(2); )==""\n"
255R"==(FMA_I_LOOP(3); )==""\n"
256R"==(} )==""\n"
257R"==(for (int h = k_align; h < k; h++) { )==""\n"
258R"==(if (irem2 > lid) )==""\n"
259R"==(a[0] = as_half2(*((global uint *)(A_ptrs[0] + h * lda) + lid)); )==""\n"
260R"==(for (int j = 0; j < 2; j++) )==""\n"
261R"==(if (jrem > j * 16) b[j] = B_ptrs[j][h]; )==""\n"
262R"==(FMA_I_LOOP(0); )==""\n"
263R"==(} )==""\n"
264R"==(} )==""\n"
265R"==(#endif )==""\n"
266R"==(global half *C2 = C + 16 * ldc; )==""\n"
267R"==(#ifdef WITH_K_UNROLL )==""\n"
268R"==(do { )==""\n"
269R"==(read_mem_fence(CLK_GLOBAL_MEM_FENCE); )==""\n"
270R"==(} while (*flag != idK); )==""\n"
271R"==(if (idK == 0) { )==""\n"
272R"==(if (beta == 0) )==""\n"
273R"==(UPDATE_C(1); )==""\n"
274R"==(else )==""\n"
275R"==(UPDATE_C(0); )==""\n"
276R"==(} else { )==""\n"
277R"==(beta = 1.0; )==""\n"
278R"==(UPDATE_C(0); )==""\n"
279R"==(} )==""\n"
280R"==(if (lid == 0) *flag = idK + 1; )==""\n"
281R"==(#else )==""\n"
282R"==(if (beta == 0) )==""\n"
283R"==(UPDATE_C(1); )==""\n"
284R"==(else )==""\n"
285R"==(UPDATE_C(0); )==""\n"
286R"==(#endif )==""\n"
287R"==(} )==""\n"
288R"==(#endif )==""\n"
289R"==(#ifdef NT )==""\n"
290R"==(__attribute__((intel_reqd_sub_group_size(16))) )==""\n"
291R"==(kernel void )==""\n"
292R"==(gen9_gemm_nocopy_f16(global half *A, global half *B, global half *C, )==""\n"
293R"==(long offset_a, long offset_b, long offset_c, int lda, int ldb, int ldc, )==""\n"
294R"==(int m, int n, int k, global float *alpha, float beta, int last_k_block, )==""\n"
295R"==(float eltwise_alpha, float eltwise_beta, float eltwise_scale )==""\n"
296R"==(#ifdef WITH_K_UNROLL )==""\n"
297R"==(, )==""\n"
298R"==(volatile global int *flag, long offset_f) { )==""\n"
299R"==(#else )==""\n"
300R"==() { )==""\n"
301R"==(#endif )==""\n"
302R"==(half2 a[4]; )==""\n"
303R"==(half2 b[4]; )==""\n"
304R"==(half c[32][2]; )==""\n"
305R"==(int idM = get_global_id(1); )==""\n"
306R"==(int idN = get_global_id(0); )==""\n"
307R"==(int lid = get_sub_group_local_id(); )==""\n"
308R"==(int idK = get_global_id(2); )==""\n"
309R"==(int nku = get_global_size(2); )==""\n"
310R"==(int i0 = idM * 32; )==""\n"
311R"==(int j00 = sub_group_broadcast(idN, 0) * 2; )==""\n"
312R"==(int j0 = j00 + lid; )==""\n"
313R"==(int irem = m - i0; )==""\n"
314R"==(int jrem = n - j0; )==""\n"
315R"==(if (irem < 0) irem = 0; )==""\n"
316R"==(if (jrem < 0) jrem = 0; )==""\n"
317R"==(int irem2 = (irem + 1) >> 1; )==""\n"
318R"==(int last_k_unroll = (idK == nku - 1); )==""\n"
319R"==(#ifdef WITH_K_UNROLL )==""\n"
320R"==(int k0 = idK * UNROLL_K; )==""\n"
321R"==(int kt = k - k0; )==""\n"
322R"==(if (kt < 0) kt = 0; )==""\n"
323R"==(if (kt > UNROLL_K) kt = UNROLL_K; )==""\n"
324R"==(A += offset_a + i0 + k0 * lda; )==""\n"
325R"==(B += offset_b + j00 + k0 * ldb; )==""\n"
326R"==(k = kt; )==""\n"
327R"==(#else )==""\n"
328R"==(A += offset_a + i0; )==""\n"
329R"==(B += offset_b + j00; )==""\n"
330R"==(#endif )==""\n"
331R"==(C += offset_c + i0 + j0 * ldc; )==""\n"
332R"==(global half *A_ptrs[4] = {A, A + lda, A + 2 * lda, A + 3 * lda}; )==""\n"
333R"==(global half *B_ptrs[4] = {B, B + ldb, B + 2 * ldb, B + 3 * ldb}; )==""\n"
334R"==(for (int y = 0; y < 32; y++) )==""\n"
335R"==(for (int z = 0; z < 2; z++) )==""\n"
336R"==(c[y][z] = 0; )==""\n"
337R"==(#ifdef WITH_K_UNROLL )==""\n"
338R"==(flag += offset_f + idM; )==""\n"
339R"==(if (idK == 0 && lid == 0) *flag = 0; )==""\n"
340R"==(#endif )==""\n"
341R"==(int k_align = k & ~3; )==""\n"
342R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
343R"==(if (irem >= 32 && sub_group_broadcast(jrem, 0) >= 32) { )==""\n"
344R"==(#endif )==""\n"
345R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
346R"==(for (int hh = 0; hh < 4; hh++) { )==""\n"
347R"==(a[hh] = as_half2(intel_sub_group_block_read( )==""\n"
348R"==((global uint *)(A_ptrs[hh] + h * lda))); )==""\n"
349R"==(b[hh] = as_half2(intel_sub_group_block_read_us2( )==""\n"
350R"==((global ushort *)(B_ptrs[hh] + h * ldb))); )==""\n"
351R"==(} )==""\n"
352R"==(FMA_I_LOOP(0); )==""\n"
353R"==(FMA_I_LOOP(1); )==""\n"
354R"==(FMA_I_LOOP(2); )==""\n"
355R"==(FMA_I_LOOP(3); )==""\n"
356R"==(} )==""\n"
357R"==(for (int h = k_align; h < k; h++) { )==""\n"
358R"==(a[0] = as_half2(intel_sub_group_block_read( )==""\n"
359R"==((global uint *)(A_ptrs[0] + h * lda))); )==""\n"
360R"==(b[0] = as_half2(intel_sub_group_block_read_us2( )==""\n"
361R"==((global ushort *)(B_ptrs[0] + h * ldb))); )==""\n"
362R"==(FMA_I_LOOP(0); )==""\n"
363R"==(} )==""\n"
364R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
365R"==(} else { )==""\n"
366R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
367R"==(for (int hh = 0; hh < 4; hh++) { )==""\n"
368R"==(if (irem2 > lid) )==""\n"
369R"==(a[hh] = as_half2( )==""\n"
370R"==(*((global uint *)(A_ptrs[hh] + h * lda) + lid)); )==""\n"
371R"==(if (jrem > 0) b[hh].s0 = B_ptrs[hh][h * ldb + lid]; )==""\n"
372R"==(if (jrem > 16) b[hh].s1 = B_ptrs[hh][h * ldb + lid + 16]; )==""\n"
373R"==(} )==""\n"
374R"==(FMA_I_LOOP(0); )==""\n"
375R"==(FMA_I_LOOP(1); )==""\n"
376R"==(FMA_I_LOOP(2); )==""\n"
377R"==(FMA_I_LOOP(3); )==""\n"
378R"==(} )==""\n"
379R"==(for (int h = k_align; h < k; h++) { )==""\n"
380R"==(if (irem2 > lid) )==""\n"
381R"==(a[0] = as_half2(*((global uint *)(A_ptrs[0] + h * lda) + lid)); )==""\n"
382R"==(if (jrem > 0) b[0].s0 = B_ptrs[0][h * ldb + lid]; )==""\n"
383R"==(if (jrem > 16) b[0].s1 = B_ptrs[0][h * ldb + lid + 16]; )==""\n"
384R"==(FMA_I_LOOP(0); )==""\n"
385R"==(} )==""\n"
386R"==(} )==""\n"
387R"==(#endif )==""\n"
388R"==(global half *C2 = C + 16 * ldc; )==""\n"
389R"==(#ifdef WITH_K_UNROLL )==""\n"
390R"==(do { )==""\n"
391R"==(read_mem_fence(CLK_GLOBAL_MEM_FENCE); )==""\n"
392R"==(} while (*flag != idK); )==""\n"
393R"==(if (idK == 0) { )==""\n"
394R"==(if (beta == 0) )==""\n"
395R"==(UPDATE_C(1); )==""\n"
396R"==(else )==""\n"
397R"==(UPDATE_C(0); )==""\n"
398R"==(} else { )==""\n"
399R"==(beta = 1.0; )==""\n"
400R"==(UPDATE_C(0); )==""\n"
401R"==(} )==""\n"
402R"==(if (lid == 0) *flag = idK + 1; )==""\n"
403R"==(#else )==""\n"
404R"==(if (beta == 0) )==""\n"
405R"==(UPDATE_C(1); )==""\n"
406R"==(else )==""\n"
407R"==(UPDATE_C(0); )==""\n"
408R"==(#endif )==""\n"
409R"==(} )==""\n"
410R"==(#endif )==""\n"
411R"==(#ifdef TN )==""\n"
412R"==(__attribute__((intel_reqd_sub_group_size(16))) )==""\n"
413R"==(kernel void )==""\n"
414R"==(gen9_gemm_nocopy_f16(global half *A, global half *B, global half *C, )==""\n"
415R"==(long offset_a, long offset_b, long offset_c, int lda, int ldb, int ldc, )==""\n"
416R"==(int m, int n, int k, global float *alpha, float beta, int last_k_block, )==""\n"
417R"==(float eltwise_alpha, float eltwise_beta, float eltwise_scale )==""\n"
418R"==(#ifdef WITH_K_UNROLL )==""\n"
419R"==(, )==""\n"
420R"==(volatile global int *flag, long offset_f) { )==""\n"
421R"==(#else )==""\n"
422R"==() { )==""\n"
423R"==(#endif )==""\n"
424R"==(half4 a[2]; )==""\n"
425R"==(half4 b[2]; )==""\n"
426R"==(half c[32][2]; )==""\n"
427R"==(int idM = get_global_id(1); )==""\n"
428R"==(int idN = get_global_id(0); )==""\n"
429R"==(int lid = get_sub_group_local_id(); )==""\n"
430R"==(int idK = get_global_id(2); )==""\n"
431R"==(int nku = get_global_size(2); )==""\n"
432R"==(int i0 = idM * 32; )==""\n"
433R"==(int j0 = sub_group_broadcast(idN, 0) * 2 + lid; )==""\n"
434R"==(int irem = m - i0; )==""\n"
435R"==(int jrem = n - j0; )==""\n"
436R"==(if (irem < 0) irem = 0; )==""\n"
437R"==(if (jrem < 0) jrem = 0; )==""\n"
438R"==(int last_k_unroll = (idK == nku - 1); )==""\n"
439R"==(#ifdef WITH_K_UNROLL )==""\n"
440R"==(int k0 = idK * UNROLL_K; )==""\n"
441R"==(int kt = k - k0; )==""\n"
442R"==(if (kt < 0) kt = 0; )==""\n"
443R"==(if (kt > UNROLL_K) kt = UNROLL_K; )==""\n"
444R"==(A += offset_a + (i0 + lid) * lda + k0; )==""\n"
445R"==(B += offset_b + j0 * ldb + k0; )==""\n"
446R"==(k = kt; )==""\n"
447R"==(#else )==""\n"
448R"==(A += offset_a + (i0 + lid) * lda; )==""\n"
449R"==(B += offset_b + j0 * ldb; )==""\n"
450R"==(#endif )==""\n"
451R"==(C += offset_c + i0 + j0 * ldc; )==""\n"
452R"==(global half *A_ptrs[2] = {A, A + 16 * lda}; )==""\n"
453R"==(global half *B_ptrs[2] = {B, B + 16 * ldb}; )==""\n"
454R"==(for (int y = 0; y < 32; y++) )==""\n"
455R"==(for (int z = 0; z < 2; z++) )==""\n"
456R"==(c[y][z] = 0; )==""\n"
457R"==(#ifdef WITH_K_UNROLL )==""\n"
458R"==(flag += offset_f + idM; )==""\n"
459R"==(if (idK == 0 && lid == 0) *flag = 0; )==""\n"
460R"==(#endif )==""\n"
461R"==(int k_align = k & ~3; )==""\n"
462R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
463R"==(if (irem >= 32 && sub_group_broadcast(jrem, 0) >= 32) { )==""\n"
464R"==(#endif )==""\n"
465R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
466R"==(for (int z = 0; z < 2; z++) { )==""\n"
467R"==(a[z] = VLOAD4_ALIGNED(0, (A_ptrs[z] + h)); )==""\n"
468R"==(b[z] = VLOAD4_ALIGNED(0, (B_ptrs[z] + h)); )==""\n"
469R"==(} )==""\n"
470R"==(FMA_I_LOOP(0); )==""\n"
471R"==(FMA_I_LOOP(1); )==""\n"
472R"==(FMA_I_LOOP(2); )==""\n"
473R"==(FMA_I_LOOP(3); )==""\n"
474R"==(} )==""\n"
475R"==(for (int h = k_align; h < k; h++) { )==""\n"
476R"==(for (int z = 0; z < 2; z++) { )==""\n"
477R"==(a[z] = A_ptrs[z][h]; )==""\n"
478R"==(b[z] = B_ptrs[z][h]; )==""\n"
479R"==(} )==""\n"
480R"==(FMA_I_LOOP(0); )==""\n"
481R"==(} )==""\n"
482R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
483R"==(} else { )==""\n"
484R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
485R"==(for (int z = 0; z < 2; z++) { )==""\n"
486R"==(if (irem > (lid + z * 16)) )==""\n"
487R"==(a[z] = VLOAD4_ALIGNED(0, (A_ptrs[z] + h)); )==""\n"
488R"==(if (jrem > z * 16) b[z] = VLOAD4_ALIGNED(0, (B_ptrs[z] + h)); )==""\n"
489R"==(} )==""\n"
490R"==(FMA_I_LOOP(0); )==""\n"
491R"==(FMA_I_LOOP(1); )==""\n"
492R"==(FMA_I_LOOP(2); )==""\n"
493R"==(FMA_I_LOOP(3); )==""\n"
494R"==(} )==""\n"
495R"==(for (int h = k_align; h < k; h++) { )==""\n"
496R"==(for (int z = 0; z < 2; z++) { )==""\n"
497R"==(if (irem > (lid + z * 16)) a[z] = A_ptrs[z][h]; )==""\n"
498R"==(if (jrem > z * 16) b[z] = B_ptrs[z][h]; )==""\n"
499R"==(} )==""\n"
500R"==(FMA_I_LOOP(0); )==""\n"
501R"==(} )==""\n"
502R"==(} )==""\n"
503R"==(#endif )==""\n"
504R"==(global half *C2 = C + 16 * ldc; )==""\n"
505R"==(#ifdef WITH_K_UNROLL )==""\n"
506R"==(do { )==""\n"
507R"==(read_mem_fence(CLK_GLOBAL_MEM_FENCE); )==""\n"
508R"==(} while (*flag != idK); )==""\n"
509R"==(if (idK == 0) { )==""\n"
510R"==(if (beta == 0) )==""\n"
511R"==(UPDATE_C(1); )==""\n"
512R"==(else )==""\n"
513R"==(UPDATE_C(0); )==""\n"
514R"==(} else { )==""\n"
515R"==(beta = 1.0; )==""\n"
516R"==(UPDATE_C(0); )==""\n"
517R"==(} )==""\n"
518R"==(if (lid == 0) *flag = idK + 1; )==""\n"
519R"==(#else )==""\n"
520R"==(if (beta == 0) )==""\n"
521R"==(UPDATE_C(1); )==""\n"
522R"==(else )==""\n"
523R"==(UPDATE_C(0); )==""\n"
524R"==(#endif )==""\n"
525R"==(} )==""\n"
526R"==(#endif )==""\n"
527R"==(#ifdef TT )==""\n"
528R"==(__attribute__((intel_reqd_sub_group_size(16))) )==""\n"
529R"==(kernel void )==""\n"
530R"==(gen9_gemm_nocopy_f16(global half *A, global half *B, global half *C, )==""\n"
531R"==(long offset_a, long offset_b, long offset_c, int lda, int ldb, int ldc, )==""\n"
532R"==(int m, int n, int k, global float *alpha, float beta, int last_k_block, )==""\n"
533R"==(float eltwise_alpha, float eltwise_beta, float eltwise_scale )==""\n"
534R"==(#ifdef WITH_K_UNROLL )==""\n"
535R"==(, )==""\n"
536R"==(volatile global int *flag, long offset_f) { )==""\n"
537R"==(#else )==""\n"
538R"==() { )==""\n"
539R"==(#endif )==""\n"
540R"==(half4 a[2]; )==""\n"
541R"==(half2 b[4]; )==""\n"
542R"==(half c[32][2]; )==""\n"
543R"==(int idM = get_global_id(1); )==""\n"
544R"==(int idN = get_global_id(0); )==""\n"
545R"==(int lid = get_sub_group_local_id(); )==""\n"
546R"==(int idK = get_global_id(2); )==""\n"
547R"==(int nku = get_global_size(2); )==""\n"
548R"==(int i0 = idM * 32; )==""\n"
549R"==(int j00 = sub_group_broadcast(idN, 0) * 2; )==""\n"
550R"==(int j0 = j00 + lid; )==""\n"
551R"==(int irem = m - i0; )==""\n"
552R"==(int jrem = n - j0; )==""\n"
553R"==(if (irem < 0) irem = 0; )==""\n"
554R"==(if (jrem < 0) jrem = 0; )==""\n"
555R"==(int last_k_unroll = (idK == nku - 1); )==""\n"
556R"==(#ifdef WITH_K_UNROLL )==""\n"
557R"==(int k0 = idK * UNROLL_K; )==""\n"
558R"==(int kt = k - k0; )==""\n"
559R"==(if (kt < 0) kt = 0; )==""\n"
560R"==(if (kt > UNROLL_K) kt = UNROLL_K; )==""\n"
561R"==(A += offset_a + (i0 + lid) * lda + k0; )==""\n"
562R"==(B += offset_b + j00 + k0 * ldb; )==""\n"
563R"==(k = kt; )==""\n"
564R"==(#else )==""\n"
565R"==(A += offset_a + (i0 + lid) * lda; )==""\n"
566R"==(B += offset_b + j00; )==""\n"
567R"==(#endif )==""\n"
568R"==(C += offset_c + i0 + j0 * ldc; )==""\n"
569R"==(global half *A_ptrs[2] = {A, A + 16 * lda}; )==""\n"
570R"==(global half *B_ptrs[4] = {B, B + ldb, B + 2 * ldb, B + 3 * ldb}; )==""\n"
571R"==(for (int y = 0; y < 32; y++) )==""\n"
572R"==(for (int z = 0; z < 2; z++) )==""\n"
573R"==(c[y][z] = 0; )==""\n"
574R"==(#ifdef WITH_K_UNROLL )==""\n"
575R"==(flag += offset_f + idM; )==""\n"
576R"==(if (idK == 0 && lid == 0) *flag = 0; )==""\n"
577R"==(#endif )==""\n"
578R"==(int k_align = k & ~3; )==""\n"
579R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
580R"==(if (irem >= 32 && sub_group_broadcast(jrem, 0) >= 32) { )==""\n"
581R"==(#endif )==""\n"
582R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
583R"==(for (int z = 0; z < 2; z++) )==""\n"
584R"==(a[z] = VLOAD4_ALIGNED(0, (A_ptrs[z] + h)); )==""\n"
585R"==(for (int hh = 0; hh < 4; hh++) )==""\n"
586R"==(b[hh] = as_half2(intel_sub_group_block_read_us2( )==""\n"
587R"==((global ushort *)(B_ptrs[hh] + h * ldb))); )==""\n"
588R"==(FMA_I_LOOP(0); )==""\n"
589R"==(FMA_I_LOOP(1); )==""\n"
590R"==(FMA_I_LOOP(2); )==""\n"
591R"==(FMA_I_LOOP(3); )==""\n"
592R"==(} )==""\n"
593R"==(for (int h = k_align; h < k; h++) { )==""\n"
594R"==(for (int z = 0; z < 2; z++) )==""\n"
595R"==(a[z] = A_ptrs[z][h]; )==""\n"
596R"==(b[0] = as_half2(intel_sub_group_block_read_us2( )==""\n"
597R"==((global ushort *)(B_ptrs[0] + h * ldb))); )==""\n"
598R"==(FMA_I_LOOP(0); )==""\n"
599R"==(} )==""\n"
600R"==(#ifndef ALLOW_READ_OVERRUNS )==""\n"
601R"==(} else { )==""\n"
602R"==(for (int h = 0; h < k_align; h += 4) { )==""\n"
603R"==(for (int z = 0; z < 2; z++) )==""\n"
604R"==(if (irem > (lid + z * 16)) )==""\n"
605R"==(a[z] = VLOAD4_ALIGNED(0, (A_ptrs[z] + h)); )==""\n"
606R"==(for (int hh = 0; hh < 4; hh++) { )==""\n"
607R"==(if (jrem > 0) b[hh].s0 = B_ptrs[hh][h * ldb + lid]; )==""\n"
608R"==(if (jrem > 16) b[hh].s1 = B_ptrs[hh][h * ldb + lid + 16]; )==""\n"
609R"==(} )==""\n"
610R"==(FMA_I_LOOP(0); )==""\n"
611R"==(FMA_I_LOOP(1); )==""\n"
612R"==(FMA_I_LOOP(2); )==""\n"
613R"==(FMA_I_LOOP(3); )==""\n"
614R"==(} )==""\n"
615R"==(for (int h = k_align; h < k; h++) { )==""\n"
616R"==(for (int z = 0; z < 2; z++) )==""\n"
617R"==(a[z] = A_ptrs[z][h]; )==""\n"
618R"==(if (jrem > 0) b[0].s0 = B_ptrs[0][h * ldb + lid]; )==""\n"
619R"==(if (jrem > 16) b[0].s1 = B_ptrs[0][h * ldb + lid + 16]; )==""\n"
620R"==(FMA_I_LOOP(0); )==""\n"
621R"==(} )==""\n"
622R"==(} )==""\n"
623R"==(#endif )==""\n"
624R"==(global half *C2 = C + 16 * ldc; )==""\n"
625R"==(#ifdef WITH_K_UNROLL )==""\n"
626R"==(do { )==""\n"
627R"==(read_mem_fence(CLK_GLOBAL_MEM_FENCE); )==""\n"
628R"==(} while (*flag != idK); )==""\n"
629R"==(if (idK == 0) { )==""\n"
630R"==(if (beta == 0) )==""\n"
631R"==(UPDATE_C(1); )==""\n"
632R"==(else )==""\n"
633R"==(UPDATE_C(0); )==""\n"
634R"==(} else { )==""\n"
635R"==(beta = 1.0; )==""\n"
636R"==(UPDATE_C(0); )==""\n"
637R"==(} )==""\n"
638R"==(if (lid == 0) *flag = idK + 1; )==""\n"
639R"==(#else )==""\n"
640R"==(if (beta == 0) )==""\n"
641R"==(UPDATE_C(1); )==""\n"
642R"==(else )==""\n"
643R"==(UPDATE_C(0); )==""\n"
644R"==(#endif )==""\n"
645R"==(} )==""\n"
646R"==(#endif )==""\n"
647R"==()==";
648}
649}
650}
651}