14
14
#extension GL_KHR_shader_subgroup_arithmetic : require
15
15
#extension GL_EXT_debug_printf : enable
16
16
17
- layout(local_size_x = 32) in;
17
+ layout (local_size_x_id = 0) in;
18
+ layout (local_size_y_id = 1) in;
19
+ layout (constant_id = 1) const uint nsg = 2;
18
20
19
21
layout(binding = 0) readonly buffer tensorInA { uint8_t inA[]; };
20
22
layout(binding = 1) readonly buffer tensorInB { float inB[]; };
@@ -38,34 +40,38 @@ layout(push_constant) uniform parameter {
38
40
}
39
41
pcs;
40
42
41
- #define ELS_PER_BLOCK 32
42
- #define QS_OFFSET 2
43
- #define BLOCK_SIZE ((ELS_PER_BLOCK / 2) + QS_OFFSET)
43
+ const uint els_per_block = 32;
44
+ const uint qs_offset = 2;
45
+ const uint block_size = (els_per_block / 2) + qs_offset;
46
+
44
47
45
48
void main() {
46
- uvec3 gid = gl_GlobalInvocationID;
49
+ uvec3 gid = gl_WorkGroupID;
50
+ uvec3 lid = gl_LocalInvocationID;
51
+ gid.y = gid.y * nsg + lid.y;
47
52
48
53
uint bc_ab = pcs.ne12 > pcs.ne02 ? gid.z / (pcs.ne12 / pcs.ne02) : gid.z;
49
54
uint bc_ba = pcs.ne02 > pcs.ne12 ? gid.z / (pcs.ne02 / pcs.ne12) : gid.z;
50
55
51
-
52
56
const uint x = (gid.x*pcs.nb01 + bc_ab*pcs.nb02) + pcs.inAOff; // Based from inA
53
57
const uint y = (gid.y*pcs.nb11 + bc_ba*pcs.nb12) / 4 + pcs.inBOff; // based from inB
54
58
float sum = 0.0f;
55
- for (uint i = 0; i < pcs.ne00; i+=ELS_PER_BLOCK) {
56
- for (uint j = 0; j < ELS_PER_BLOCK / 2; j++) {
57
- const uint block_number = i / ELS_PER_BLOCK;
58
- const uint block_offset = block_number * BLOCK_SIZE;
59
+ for (uint i = gl_SubgroupInvocationID * 2; i < pcs.ne00; i+=gl_SubgroupSize * 2) {
60
+ const uint block_number = i / els_per_block;
61
+ const uint block_offset = block_number * block_size;
59
62
const float d = u8BufToFloat16(inA, x + block_offset);
63
+ const uint j = (i % els_per_block) / 2;
60
64
const uint byte_position_in_block = j;
61
- const int q0 = (inA[x+block_offset+QS_OFFSET +byte_position_in_block] & 0x0F) - 8;
62
- const int q1 = (inA[x+block_offset+QS_OFFSET +byte_position_in_block] >> 4) - 8;
65
+ const int q0 = (inA[x+block_offset+qs_offset +byte_position_in_block] & 0x0F) - 8;
66
+ const int q1 = (inA[x+block_offset+qs_offset +byte_position_in_block] >> 4) - 8;
63
67
const float dq0 = d * q0;
64
68
const float dq1 = d * q1;
65
- sum += (dq0 * float(inB[y+i+j])) + \
66
- (dq1 * float(inB[y+i+j+(ELS_PER_BLOCK/2) ]));
67
- }
69
+ const uint block_base = block_number * els_per_block;
70
+ sum += (dq0 * float(inB[y+block_base+j ])) + \
71
+ (dq1 * float(inB[y+block_base+j+(els_per_block/2)]));
68
72
}
69
73
70
- out_[gid.z*(pcs.nb2/4) + gid.y*(pcs.nb1/4) + gid.x + pcs.outOff] = sum;
74
+ const float all_sum = subgroupAdd(sum);
75
+ if (subgroupElect())
76
+ out_[gid.z*(pcs.nb2/4) + gid.y*(pcs.nb1/4) + gid.x + pcs.outOff] = all_sum;
71
77
}
0 commit comments