1#version 450
 2
 3#include "types.glsl"
 4#include "generic_unary_head.glsl"
 5
 6layout(local_size_x = 512, local_size_y = 1, local_size_z = 1) in;
 7
 8void main() {
 9    const uint idx = get_idx();
10
11    if (idx >= p.ne) {
12        return;
13    }
14
15    // Destination multi-index (inlined dst_idx)
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    const uint d_idx = i13*p.nb13 + i12*p.nb12 + i11*p.nb11 + i10*p.nb10;
23
24    // Accumulate from sources
25    A_TYPE acc = A_TYPE(0);
26    for (uint i3 = i13; i3 < p.ne03; i3 += p.ne13) {
27        for (uint i2 = i12; i2 < p.ne02; i2 += p.ne12) {
28            for (uint i1 = i11; i1 < p.ne01; i1 += p.ne11) {
29                for (uint i0 = i10; i0 < p.ne00; i0 += p.ne10) {
30                    acc += data_a[i3*p.nb03 + i2*p.nb02 + i1*p.nb01 + i0*p.nb00];
31                }
32            }
33        }
34    }
35
36    data_d[get_doffset() + d_idx] = D_TYPE(acc);
37}