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);
屏障不是互斥体。您有一个数据竞争,其中两个工作项试图同时写入同一个变量。
您将不得不使用原子或重新设计您的代码。
我这里有一个测试内核,我正在从 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);
屏障不是互斥体。您有一个数据竞争,其中两个工作项试图同时写入同一个变量。
您将不得不使用原子或重新设计您的代码。