1#version 450
 2
 3#include "rte.glsl"
 4#include "types.glsl"
 5#include "generic_unary_head.glsl"
 6
 7layout(local_size_x = 512, local_size_y = 1, local_size_z = 1) in;
 8
 9void main() {
10    const uint idx = get_idx();
11
12    if (idx >= p.ne) {
13        return;
14    }
15
16    const uint i13 = fastdiv(idx, p.ne1_012mp, p.ne1_012L);
17    const uint i13_offset = i13 * p.ne12*p.ne11*p.ne10;
18    const uint i12 = fastdiv(idx - i13_offset, p.ne1_01mp, p.ne1_01L);
19    const uint i12_offset = i12*p.ne11*p.ne10;
20    const uint i11 = fastdiv(idx - i13_offset - i12_offset, p.ne1_0mp, p.ne1_0L);
21    const uint i10 = idx - i13_offset - i12_offset - i11*p.ne10;
22
23    if (i10 == i11) {
24        const float val = float(data_a[get_aoffset() + i13*p.nb03 + i12*p.nb02 + 0*p.nb01 + i10*p.nb00]);
25        data_d[get_doffset() + dst_idx(idx)] = D_TYPE(val);
26    } else {
27        data_d[get_doffset() + dst_idx(idx)] = D_TYPE(0);
28    }
29}