有人可以解释这个 OpenCL 程序有什么问题吗？

Question

我已经包含了主机程序的主要部分，我怀疑这里是不正确的：

http://pastebin.com/qVkv9E11

我对指针还不是很擅长，我认为我可能错误地分配了一些变量。

这是内核程序，它应该可以让您了解我的程序正在尝试做什么：

    const char *KernelSource =           "\n"
"__kernel void sumElements(           \n"
"   __global float* input,            \n"
"   __global float output,            \n"
"   __global int N)                   \n"
"{                                    \n"
"   int i = get_global_id(0);         \n"
"   if(i < N)                         \n"
"       output += input[i];           \n"
"}                                    \n"
"\n";

也许这是导致错误的原因，因为我从未尝试过 SIMT 写入上述一个变量。有可能做这样的事情吗？我需要得到数组中所有元素的总和。

Answer 1

如果你想真正读回输出的值，那么你也需要将它声明为一个指针。现在输出的值作为内核参数被复制进来，但是在内核结束后，您对其所做的任何更改都会被忽略。

所以，将__global float output 更改为__global float* output。然后在您的内核更改中：

if(i < N)
    output += input[i];

到

if(i < N)
    *output += input[i];

您可能需要更改分配缓冲区的方式才能使其正常工作，但是自从我在 OpenCL 中完成此操作以来已经很长时间了，而且我现在找到的文档没有显示任何明显的缓冲区中的错误。

这里警告一句：加法不是原子操作。使用这种设置，总是会发生两个或多个线程读取 *output 的值，然后尝试在不同阶段将 *output + 1 写回它。因此，*output 的值会小于应有的值。

要解决此问题，您需要使用 OpenCL atomic operations.