1#version 450
 2
 3#include "types.glsl"
 4
 5layout (push_constant) uniform parameter
 6{
 7    uint ne;
 8    uint ne00; uint ne01; uint ne02; uint ne03; uint nb00; uint nb01; uint nb02; uint nb03;
 9    uint ne10; uint ne11; uint ne12; uint ne13; uint nb10; uint nb11; uint nb12; uint nb13;
10    uint misalign_offsets;
11    uint circular;
12
13    uint lp0; uint rp0;
14    uint lp1; uint rp1;
15    uint lp2; uint rp2;
16    uint lp3; uint rp3;
17} p;
18
19uint get_aoffset() { return p.misalign_offsets >> 16; }
20uint get_doffset() { return p.misalign_offsets & 0xFFFF; }
21
22uint wrap_around(int coord, uint size) {
23    return (uint(coord + int(size))) % size; // add size to avoid issues with negative
24}
25
26layout (binding = 0) readonly buffer A {A_TYPE data_a[];};
27layout (binding = 1) writeonly buffer D {D_TYPE data_d[];};
28
29layout(local_size_x = 512, local_size_y = 1, local_size_z = 1) in;
30
31void main() {
32    const uint idx = gl_GlobalInvocationID.z * 262144 + gl_GlobalInvocationID.y * 512 + gl_GlobalInvocationID.x;
33
34    if (idx >= p.ne) {
35        return;
36    }
37
38    const uint i3 = idx / (p.ne12*p.ne11*p.ne10);
39    const uint i3_offset = i3 * p.ne12*p.ne11*p.ne10;
40    const uint i2 = (idx - i3_offset) / (p.ne11*p.ne10);
41    const uint i2_offset = i2*p.ne11*p.ne10;
42    const uint i1 = (idx - i3_offset - i2_offset) / p.ne10;
43    const uint i0 = idx - i3_offset - i2_offset - i1*p.ne10;
44
45    const uint src0_idx = (i3 - p.lp3)*p.nb03 + (i2 - p.lp2)*p.nb02 + (i1 - p.lp1)*p.nb01 + (i0 - p.lp0)*p.nb00;
46    const uint dst_idx = i3*p.nb13 + i2*p.nb12 + i1*p.nb11 + i0*p.nb10;
47
48    if (p.circular != 0u) {
49        const uint ci0 = wrap_around(int(i0) - int(p.lp0), p.ne00);
50        const uint ci1 = wrap_around(int(i1) - int(p.lp1), p.ne01);
51        const uint ci2 = wrap_around(int(i2) - int(p.lp2), p.ne02);
52        const uint ci3 = wrap_around(int(i3) - int(p.lp3), p.ne03);
53        const uint circular_src_idx = ci3*p.nb03 + ci2*p.nb02 + ci1*p.nb01 + ci0*p.nb00;
54        data_d[get_doffset() + dst_idx] = D_TYPE(data_a[get_aoffset() + circular_src_idx]);
55    } else {
56        const bool is_src0 = i0 >= p.lp0 && i0 < p.ne10 - p.rp0 &&
57                             i1 >= p.lp1 && i1 < p.ne11 - p.rp1 &&
58                             i2 >= p.lp2 && i2 < p.ne12 - p.rp2 &&
59                             i3 >= p.lp3 && i3 < p.ne13 - p.rp3;
60        data_d[get_doffset() + dst_idx] = D_TYPE(is_src0 ? data_a[get_aoffset() + src0_idx] : 0.0f);
61    }
62
63
64}