1 | namespace dnnl { |
2 | namespace impl { |
3 | namespace gpu { |
4 | namespace ocl { |
5 | const char *ref_softmax_kernel = R"==(/******************************************************************************* )==" "\n" |
6 | R"==(* Copyright 2019-2022 Intel Corporation )==" "\n" |
7 | R"==(* )==" "\n" |
8 | R"==(* Licensed under the Apache License, Version 2.0 (the "License"); )==" "\n" |
9 | R"==(* you may not use this file except in compliance with the License. )==" "\n" |
10 | R"==(* You may obtain a copy of the License at )==" "\n" |
11 | R"==(* )==" "\n" |
12 | R"==(* http: )==" "\n" |
13 | R"==(* )==" "\n" |
14 | R"==(* Unless required by applicable law or agreed to in writing, software )==" "\n" |
15 | R"==(* distributed under the License is distributed on an "AS IS" BASIS, )==" "\n" |
16 | R"==(* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. )==" "\n" |
17 | R"==(* See the License for the specific language governing permissions and )==" "\n" |
18 | R"==(* limitations under the License. )==" "\n" |
19 | R"==(*******************************************************************************/ )==" "\n" |
20 | R"==(#include "gpu/ocl/ocl_types.h" )==" "\n" |
21 | R"==(#define CONCAt2(a, b) a##b )==" "\n" |
22 | R"==(#define CONCAT2(a, b) CONCAt2(a, b) )==" "\n" |
23 | R"==(#define IS_HALF_half 1 )==" "\n" |
24 | R"==(#define IS_HALF(dt) CONCAT2(IS_HALF_, dt) )==" "\n" |
25 | R"==(#if IS_FWD )==" "\n" |
26 | R"==(#define DD(i) CONCAt2(DST_D, i) )==" "\n" |
27 | R"==(#elif IS_BWD )==" "\n" |
28 | R"==(#define DD(i) CONCAt2(SRC_D, i) )==" "\n" |
29 | R"==(#else )==" "\n" |
30 | R"==(#error unsupported data parameter )==" "\n" |
31 | R"==(#endif )==" "\n" |
32 | R"==(#define OFF(dim, idx) \ )==" "\n" |
33 | R"==((dim % CONCAT2(DATA_B, idx)) * CONCAT2(DATA_SB, idx) \ )==" "\n" |
34 | R"==(+ (dim / CONCAT2(DATA_B, idx)) * CONCAT2(DATA_S, idx) )==" "\n" |
35 | R"==(#if SOFTMAX_AXIS_IDX == 0 )==" "\n" |
36 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
37 | R"==(OFF(softmax_dim, 0) + OFF(dim0, 1) + OFF(dim1, 2) + OFF(dim2, 3) \ )==" "\n" |
38 | R"==(+ OFF(dim3, 4) + OFF(dim4, 5) )==" "\n" |
39 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
40 | R"==(softmax_dim >= DD(0) || dim0 >= DD(1) || dim1 >= DD(2) || dim2 >= DD(3) \ )==" "\n" |
41 | R"==(|| dim3 >= DD(4) || dim4 >= DD(5) )==" "\n" |
42 | R"==(#elif SOFTMAX_AXIS_IDX == 1 )==" "\n" |
43 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
44 | R"==(OFF(dim0, 0) + OFF(softmax_dim, 1) + OFF(dim1, 2) + OFF(dim2, 3) \ )==" "\n" |
45 | R"==(+ OFF(dim3, 4) + OFF(dim4, 5) )==" "\n" |
46 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
47 | R"==(dim0 >= DD(0) || softmax_dim >= DD(1) || dim1 >= DD(2) || dim2 >= DD(3) \ )==" "\n" |
48 | R"==(|| dim3 >= DD(4) || dim4 >= DD(5) )==" "\n" |
49 | R"==(#elif SOFTMAX_AXIS_IDX == 2 )==" "\n" |
50 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
51 | R"==(OFF(dim0, 0) + OFF(dim1, 1) + OFF(softmax_dim, 2) + OFF(dim2, 3) \ )==" "\n" |
52 | R"==(+ OFF(dim3, 4) + OFF(dim4, 5) )==" "\n" |
53 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
54 | R"==(dim0 >= DD(0) || dim1 >= DD(1) || softmax_dim >= DD(2) || dim2 >= DD(3) \ )==" "\n" |
55 | R"==(|| dim3 >= DD(4) || dim4 >= DD(5) )==" "\n" |
56 | R"==(#elif SOFTMAX_AXIS_IDX == 3 )==" "\n" |
57 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
58 | R"==(OFF(dim0, 0) + OFF(dim1, 1) + OFF(dim2, 2) + OFF(softmax_dim, 3) \ )==" "\n" |
59 | R"==(+ OFF(dim3, 4) + OFF(dim4, 5) )==" "\n" |
60 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
61 | R"==(dim0 >= DD(0) || dim1 >= DD(1) || dim2 >= DD(2) || softmax_dim >= DD(3) \ )==" "\n" |
62 | R"==(|| dim3 >= DD(4) || dim4 >= DD(5) )==" "\n" |
63 | R"==(#elif SOFTMAX_AXIS_IDX == 4 )==" "\n" |
64 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
65 | R"==(OFF(dim0, 0) + OFF(dim1, 1) + OFF(dim2, 2) + OFF(dim3, 3) \ )==" "\n" |
66 | R"==(+ OFF(softmax_dim, 4) + OFF(dim4, 5) )==" "\n" |
67 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
68 | R"==(dim0 >= DD(0) || dim1 >= DD(1) || dim2 >= DD(2) || dim3 >= DD(3) \ )==" "\n" |
69 | R"==(|| softmax_dim >= DD(4) || dim4 >= DD(5) )==" "\n" |
70 | R"==(#elif SOFTMAX_AXIS_IDX == 5 )==" "\n" |
71 | R"==(#define DATA_OFF(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
72 | R"==(OFF(dim0, 0) + OFF(dim1, 1) + OFF(dim2, 2) + OFF(dim3, 3) + OFF(dim4, 4) \ )==" "\n" |
73 | R"==(+ OFF(softmax_dim, 5) )==" "\n" |
74 | R"==(#define NEEDS_PADDING(dim0, dim1, dim2, dim3, dim4, softmax_dim) \ )==" "\n" |
75 | R"==(dim0 >= DD(0) || dim1 >= DD(1) || dim2 >= DD(2) || dim3 >= DD(3) \ )==" "\n" |
76 | R"==(|| dim4 >= DD(4) || softmax_dim >= DD(5) )==" "\n" |
77 | R"==(#else )==" "\n" |
78 | R"==(#error unsupported softmax dimension )==" "\n" |
79 | R"==(#endif )==" "\n" |
80 | R"==(#if IS_FWD )==" "\n" |
81 | R"==(#if SUB_GROUP_SIZE == 16 )==" "\n" |
82 | R"==(__attribute__((reqd_work_group_size(GROUP_SIZE, 1, 1))) )==" "\n" |
83 | R"==(__attribute__((intel_reqd_sub_group_size(SUB_GROUP_SIZE))) )==" "\n" |
84 | R"==(#endif )==" "\n" |
85 | R"==(__kernel void )==" "\n" |
86 | R"==(ref_softmax_fwd_generic(__global SRC_DATA_T *src, __global DATA_T *dst, )==" "\n" |
87 | R"==(__global float *src_scale, __global float *dst_scale) { )==" "\n" |
88 | R"==(const int dim[] = { )==" "\n" |
89 | R"==((get_global_id(0) / GROUP_SIZE) % BLOCK_0, )==" "\n" |
90 | R"==(get_global_id(1) % BLOCK_1, )==" "\n" |
91 | R"==(get_global_id(2) % BLOCK_2, )==" "\n" |
92 | R"==((get_global_id(0) / GROUP_SIZE) / BLOCK_0, )==" "\n" |
93 | R"==(get_global_id(1) / BLOCK_1, )==" "\n" |
94 | R"==(get_global_id(2) / BLOCK_2, )==" "\n" |
95 | R"==(}; )==" "\n" |
96 | R"==(float scale = 1.0f; )==" "\n" |
97 | R"==(#if WITH_SRC_SCALES )==" "\n" |
98 | R"==(scale *= src_scale[0]; )==" "\n" |
99 | R"==(#endif )==" "\n" |
100 | R"==(#if WITH_DST_SCALES )==" "\n" |
101 | R"==(scale /= dst_scale[0]; )==" "\n" |
102 | R"==(#endif )==" "\n" |
103 | R"==(#if SUB_GROUP_SIZE == 16 )==" "\n" |
104 | R"==(const int local_id = get_local_id(0); )==" "\n" |
105 | R"==(const int begin = local_id * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
106 | R"==(const int end = (local_id == GROUP_SIZE - 1) )==" "\n" |
107 | R"==(? SOFTMAX_AXIS )==" "\n" |
108 | R"==(: (local_id + 1) * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
109 | R"==(#if SOFTMAX_AXIS - (GROUP_SIZE - 1) * (SOFTMAX_AXIS / GROUP_SIZE) \ )==" "\n" |
110 | R"==(> SOFTMAX_AXIS / GROUP_SIZE )==" "\n" |
111 | R"==(const int buf_size )==" "\n" |
112 | R"==(= SOFTMAX_AXIS - (GROUP_SIZE - 1) * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
113 | R"==(#else )==" "\n" |
114 | R"==(const int buf_size = SOFTMAX_AXIS / GROUP_SIZE; )==" "\n" |
115 | R"==(#endif )==" "\n" |
116 | R"==(#else )==" "\n" |
117 | R"==(const int begin = 0; )==" "\n" |
118 | R"==(const int end = SOFTMAX_AXIS; )==" "\n" |
119 | R"==(const int buf_size = SOFTMAX_AXIS; )==" "\n" |
120 | R"==(#endif )==" "\n" |
121 | R"==(#if IS_HALF(SRC_DATA_T) == 1 && IS_HALF(DST_DATA_TYPE) == 1 )==" "\n" |
122 | R"==(typedef half acc_t; )==" "\n" |
123 | R"==(const acc_t acc_max = HALF_MAX; )==" "\n" |
124 | R"==(const acc_t acc_zero = 0.h; )==" "\n" |
125 | R"==(#else )==" "\n" |
126 | R"==(typedef float acc_t; )==" "\n" |
127 | R"==(const acc_t acc_max = FLT_MAX; )==" "\n" |
128 | R"==(const acc_t acc_zero = 0.f; )==" "\n" |
129 | R"==(#endif )==" "\n" |
130 | R"==(acc_t d[buf_size]; )==" "\n" |
131 | R"==(acc_t max_ = -acc_max; )==" "\n" |
132 | R"==(acc_t denom_ = acc_zero; )==" "\n" |
133 | R"==(if (!(NEEDS_PADDING(dim[0], dim[1], dim[2], dim[3], dim[4], begin))) { )==" "\n" |
134 | R"==(for (int i = begin; i < end && i < DD(SOFTMAX_AXIS_IDX); ++i) { )==" "\n" |
135 | R"==(size_t data_off )==" "\n" |
136 | R"==(= DATA_OFF(dim[0], dim[1], dim[2], dim[3], dim[4], i); )==" "\n" |
137 | R"==(d[i - begin] = SRC_TO_REF(src[data_off]); )==" "\n" |
138 | R"==(max_ = max(max_, d[i - begin]); )==" "\n" |
139 | R"==(} )==" "\n" |
140 | R"==(} )==" "\n" |
141 | R"==(#if SUB_GROUP_SIZE == 16 )==" "\n" |
142 | R"==(#if GROUP_SIZE == SUB_GROUP_SIZE )==" "\n" |
143 | R"==(max_ = sub_group_reduce_max(max_); )==" "\n" |
144 | R"==(#else )==" "\n" |
145 | R"==(max_ = work_group_reduce_max(max_); )==" "\n" |
146 | R"==(#endif )==" "\n" |
147 | R"==(#endif )==" "\n" |
148 | R"==(if (!(NEEDS_PADDING(dim[0], dim[1], dim[2], dim[3], dim[4], begin))) { )==" "\n" |
149 | R"==(for (int i = begin; i < end && i < DD(SOFTMAX_AXIS_IDX); ++i) { )==" "\n" |
150 | R"==(#if LOGSOFTMAX )==" "\n" |
151 | R"==(denom_ += exp(d[i - begin] - max_); )==" "\n" |
152 | R"==(#else )==" "\n" |
153 | R"==(d[i - begin] = exp(d[i - begin] - max_); )==" "\n" |
154 | R"==(denom_ += d[i - begin]; )==" "\n" |
155 | R"==(#endif )==" "\n" |
156 | R"==(} )==" "\n" |
157 | R"==(} )==" "\n" |
158 | R"==(#if SUB_GROUP_SIZE == 16 )==" "\n" |
159 | R"==(#if GROUP_SIZE == SUB_GROUP_SIZE )==" "\n" |
160 | R"==(denom_ = sub_group_reduce_add(denom_); )==" "\n" |
161 | R"==(#else )==" "\n" |
162 | R"==(denom_ = work_group_reduce_add(denom_); )==" "\n" |
163 | R"==(#endif )==" "\n" |
164 | R"==(#endif )==" "\n" |
165 | R"==(#if LOGSOFTMAX )==" "\n" |
166 | R"==(denom_ = log(denom_); )==" "\n" |
167 | R"==(#else )==" "\n" |
168 | R"==(denom_ = 1.0 / denom_; )==" "\n" |
169 | R"==(#endif )==" "\n" |
170 | R"==(for (int i = begin; i < end; ++i) { )==" "\n" |
171 | R"==(size_t data_off = DATA_OFF(dim[0], dim[1], dim[2], dim[3], dim[4], i); )==" "\n" |
172 | R"==(if (NEEDS_PADDING(dim[0], dim[1], dim[2], dim[3], dim[4], i)) { )==" "\n" |
173 | R"==(dst[data_off] = REF_TO_DST(acc_zero); )==" "\n" |
174 | R"==(} else { )==" "\n" |
175 | R"==(float unscaled; )==" "\n" |
176 | R"==(#if LOGSOFTMAX )==" "\n" |
177 | R"==(unscaled = d[i - begin] - max_ - denom_; )==" "\n" |
178 | R"==(#else )==" "\n" |
179 | R"==(unscaled = d[i - begin] * denom_; )==" "\n" |
180 | R"==(#endif )==" "\n" |
181 | R"==(#if DT_S8 == 1 || DT_U8 == 1 )==" "\n" |
182 | R"==(dst[data_off] = REF_TO_DST(round(scale * unscaled)); )==" "\n" |
183 | R"==(#else )==" "\n" |
184 | R"==(dst[data_off] = REF_TO_DST(scale * unscaled); )==" "\n" |
185 | R"==(#endif )==" "\n" |
186 | R"==(} )==" "\n" |
187 | R"==(} )==" "\n" |
188 | R"==(} )==" "\n" |
189 | R"==(#endif )==" "\n" |
190 | R"==(#if IS_BWD )==" "\n" |
191 | R"==(__attribute__((reqd_work_group_size(GROUP_SIZE, 1, 1))) )==" "\n" |
192 | R"==(__attribute__((intel_reqd_sub_group_size(SUB_GROUP_SIZE))) )==" "\n" |
193 | R"==(__kernel void )==" "\n" |
194 | R"==(ref_softmax_bwd_generic(__global DST_DATA_T *dst, __global SRC_DATA_T *diff_src, )==" "\n" |
195 | R"==(__global DST_DATA_T *diff_dst) { )==" "\n" |
196 | R"==(const int dim[] = { )==" "\n" |
197 | R"==((get_global_id(0) / GROUP_SIZE) % BLOCK_0, )==" "\n" |
198 | R"==(get_global_id(1) % BLOCK_1, )==" "\n" |
199 | R"==(get_global_id(2) % BLOCK_2, )==" "\n" |
200 | R"==((get_global_id(0) / GROUP_SIZE) / BLOCK_0, )==" "\n" |
201 | R"==(get_global_id(1) / BLOCK_1, )==" "\n" |
202 | R"==(get_global_id(2) / BLOCK_2, )==" "\n" |
203 | R"==(}; )==" "\n" |
204 | R"==(const int local_id = get_local_id(0); )==" "\n" |
205 | R"==(const int begin = local_id * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
206 | R"==(const int end = (local_id == GROUP_SIZE - 1) )==" "\n" |
207 | R"==(? SOFTMAX_AXIS )==" "\n" |
208 | R"==(: (local_id + 1) * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
209 | R"==(#if SOFTMAX_AXIS - (GROUP_SIZE - 1) * (SOFTMAX_AXIS / GROUP_SIZE) \ )==" "\n" |
210 | R"==(> SOFTMAX_AXIS / GROUP_SIZE )==" "\n" |
211 | R"==(const int buf_size )==" "\n" |
212 | R"==(= SOFTMAX_AXIS - (GROUP_SIZE - 1) * (SOFTMAX_AXIS / GROUP_SIZE); )==" "\n" |
213 | R"==(#else )==" "\n" |
214 | R"==(const int buf_size = SOFTMAX_AXIS / GROUP_SIZE; )==" "\n" |
215 | R"==(#endif )==" "\n" |
216 | R"==(#if IS_HALF(SRC_DATA_T) == 1 && IS_HALF(DST_DATA_TYPE) == 1 )==" "\n" |
217 | R"==(typedef half acc_t; )==" "\n" |
218 | R"==(const acc_t acc_zero = 0.h; )==" "\n" |
219 | R"==(#else )==" "\n" |
220 | R"==(typedef float acc_t; )==" "\n" |
221 | R"==(const acc_t acc_zero = 0.f; )==" "\n" |
222 | R"==(#endif )==" "\n" |
223 | R"==(acc_t diff_d[buf_size]; )==" "\n" |
224 | R"==(acc_t d[buf_size]; )==" "\n" |
225 | R"==(acc_t sbr = acc_zero; )==" "\n" |
226 | R"==(if (!(NEEDS_PADDING(dim[0], dim[1], dim[2], dim[3], dim[4], begin))) { )==" "\n" |
227 | R"==(for (int i = begin; i < end && i < DD(SOFTMAX_AXIS_IDX); ++i) { )==" "\n" |
228 | R"==(size_t idx = DATA_OFF(dim[0], dim[1], dim[2], dim[3], dim[4], i); )==" "\n" |
229 | R"==(diff_d[i - begin] = DST_TO_REF(diff_dst[idx]); )==" "\n" |
230 | R"==(d[i - begin] = DST_TO_REF(dst[idx]); )==" "\n" |
231 | R"==(#if LOGSOFTMAX )==" "\n" |
232 | R"==(sbr += diff_d[i - begin]; )==" "\n" |
233 | R"==(#else )==" "\n" |
234 | R"==(sbr += diff_d[i - begin] * d[i - begin]; )==" "\n" |
235 | R"==(#endif )==" "\n" |
236 | R"==(} )==" "\n" |
237 | R"==(} )==" "\n" |
238 | R"==(#if GROUP_SIZE == SUB_GROUP_SIZE )==" "\n" |
239 | R"==(sbr = sub_group_reduce_add(sbr); )==" "\n" |
240 | R"==(#else )==" "\n" |
241 | R"==(sbr = work_group_reduce_add(sbr); )==" "\n" |
242 | R"==(#endif )==" "\n" |
243 | R"==(for (int i = begin; i < end; ++i) { )==" "\n" |
244 | R"==(size_t idx = DATA_OFF(dim[0], dim[1], dim[2], dim[3], dim[4], i); )==" "\n" |
245 | R"==(if (NEEDS_PADDING(dim[0], dim[1], dim[2], dim[3], dim[4], i)) { )==" "\n" |
246 | R"==(diff_src[idx] = REF_TO_SRC(acc_zero); )==" "\n" |
247 | R"==(} else { )==" "\n" |
248 | R"==(#if LOGSOFTMAX )==" "\n" |
249 | R"==(diff_src[idx] )==" "\n" |
250 | R"==(= REF_TO_SRC(diff_d[i - begin] - exp(d[i - begin]) * sbr); )==" "\n" |
251 | R"==(#else )==" "\n" |
252 | R"==(acc_t inner_data = diff_d[i - begin] - sbr; )==" "\n" |
253 | R"==(diff_src[idx] = REF_TO_SRC(d[i - begin] * inner_data); )==" "\n" |
254 | R"==(#endif )==" "\n" |
255 | R"==(} )==" "\n" |
256 | R"==(} )==" "\n" |
257 | R"==(} )==" "\n" |
258 | R"==(#endif )==" "\n" |
259 | R"==()==" ; |
260 | } |
261 | } |
262 | } |
263 | } |