1#version 450
 2
 3#include "types.glsl"
 4#include "generic_binary_head.glsl"
 5
 6layout(local_size_x = 512, local_size_y = 1, local_size_z = 1) in;
 7
 8void main() {
 9    const uint idx = gl_GlobalInvocationID.z * 262144 + gl_GlobalInvocationID.y * 512 + gl_GlobalInvocationID.x;
10    const int dim = p.param3;
11
12    if (idx >= p.ne) {
13        return;
14    }
15
16    const uint i3 = idx / (p.ne22*p.ne21*p.ne20);
17    const uint i3_offset = i3 * p.ne22*p.ne21*p.ne20;
18    const uint i2 = (idx - i3_offset) / (p.ne21*p.ne20);
19    const uint i2_offset = i2*p.ne21*p.ne20;
20    const uint i1 = (idx - i3_offset - i2_offset) / p.ne20;
21    const uint i0 = idx - i3_offset - i2_offset - i1*p.ne20;
22
23    uint o[4] = {0, 0, 0, 0};
24    o[dim] = dim == 0 ? p.ne00 : (dim == 1 ? p.ne01 : (dim == 2 ? p.ne02 : p.ne03));
25
26    const uint src0_idx = i3*p.nb03 + i2*p.nb02 + i1*p.nb01 + i0*p.nb00;
27    const uint src1_idx = (i3 - o[3])*p.nb13 + (i2 - o[2])*p.nb12 + (i1 - o[1])*p.nb11 + (i0 - o[0])*p.nb10;
28    const uint dst_idx = i3*p.nb23 + i2*p.nb22 + i1*p.nb21 + i0*p.nb20;
29
30    const bool is_src0 = i0 < p.ne00 && i1 < p.ne01 && i2 < p.ne02 && i3 < p.ne03;
31
32#ifndef OPTIMIZATION_ERROR_WORKAROUND
33    data_d[get_doffset() + dst_idx] = D_TYPE(is_src0 ? data_a[get_aoffset() + src0_idx] : data_b[get_boffset() + src1_idx]);
34#else
35    if (is_src0) {
36        data_d[get_doffset() + dst_idx] = data_a[get_aoffset() + src0_idx];
37    } else {
38        data_d[get_doffset() + dst_idx] = data_b[get_boffset() + src1_idx];
39    }
40#endif
41}