1#version 450
 2
 3#extension GL_EXT_control_flow_attributes : enable
 4
 5#include "types.glsl"
 6#include "generic_binary_head.glsl"
 7#include "dequant_funcs.glsl"
 8
 9layout(local_size_x = 512, local_size_y = 1, local_size_z = 1) in;
10
11void main() {
12    const uint i00 = (gl_GlobalInvocationID.x)*2;
13
14#ifdef NEEDS_INIT_IQ_SHMEM
15    init_iq_shmem(gl_WorkGroupSize);
16#endif
17
18    if (i00 >= p.ne00) {
19        return;
20    }
21
22    uint gid_z = gl_GlobalInvocationID.z;
23    while (gid_z < p.ne11 * p.ne12) {
24        uint gid_y = gl_GlobalInvocationID.y;
25        while (gid_y < p.ne10) {
26            const uint i10 = gid_y;
27            const uint i11 = gid_z / p.ne12;
28            const uint i12 = gid_z % p.ne12;
29
30            const uint i01 = data_b[i10*p.nb10 + i11*p.nb11 + i12*p.nb12];
31
32            const uint a_offset = i01*p.nb01 + i11*p.nb02 + i12*p.nb03;
33            const uint d_offset = i10*p.nb21 + i11*p.nb22 + i12*p.nb23;
34
35            const uint ib = a_offset + i00/QUANT_K; // block index
36            const uint iqs = (i00%QUANT_K)/QUANT_R; // quant index
37            const uint iybs = i00 - i00%QUANT_K; // dst block start index
38            const uint y_offset = QUANT_R == 1 ? 1 : QUANT_K/2;
39
40            vec2 v = dequantize(ib, iqs, 0);
41            const vec2 dm = get_dm(ib, 0);
42            v = v * dm.x + dm.y;
43
44            data_d[d_offset + iybs + iqs           ] = D_TYPE(v.x);
45            data_d[d_offset + iybs + iqs + y_offset] = D_TYPE(v.y);
46
47            gid_y += gl_WorkGroupSize.y * gl_NumWorkGroups.y;
48        }
49        gid_z += gl_WorkGroupSize.z * gl_NumWorkGroups.z;
50    }
51}