OpenCL 内核 _local 内存行为不正确

OpenCL Kernel _local memory not behaving correctly

我这里有一个测试内核,我正在从 Julia OpenCL API 调用它。我从 Julia 调用它这一事实并不重要,这只是我在 运行 OpenCL 中使用的,这里是代码:

using OpenCL
const cl = OpenCL

device, ctx, queue = cl.create_compute_context()

C_buff = cl.Buffer(Float32, ctx, :w, 2)

const testkernel = """
kernel void test(global float *C)
{
    int gid = get_global_id(0);
    int lid = get_local_id(0);

    local float x;

    if (lid == 0)
    {
        x = 0.0f;
    }
    barrier(CLK_LOCAL_MEM_FENCE);

    x += 1.0f;

    barrier(CLK_LOCAL_MEM_FENCE);

    if (lid == 0)
    {
        C[gid / 2] = x;
    }
}
"""

program = cl.Program(ctx, source=testkernel) |> cl.build!
kernel = cl.Kernel(program, "test")
cl.call(queue, kernel, 4, 2, C_buff)
cl.read(queue, C_buff)

我想不通的是这个 return 是一个向量 [1.0,1.0],而它似乎应该 return 向量 [2.0,2.0]。因为基本上我有 4 个工作项分成两个工作组(每个工作组包含 2 个工作项)。

为每个工作组实例化一个本地浮点数x,每个工作组中的第一个工作项将其设置为0。然后工作组中的每个工作项将其加1 ,并且由于每个工作组中有两个工作项,因此它应该是 2,但是当我 return C 时,我得到的是一个向量。

barrier(CLK_LOCAL_MEM_FENCE);

x += 1.0f;

barrier(CLK_LOCAL_MEM_FENCE);

屏障不是互斥体。您有一个数据竞争,其中两个工作项试图同时写入同一个变量。

您将不得不使用原子或重新设计您的代码。