Question

我有以下例程，它应该根据索引中找到的索引将元素从src矩阵复制到dst矩阵。索引计算正确但dst未更新。我错过了什么？

__kernel void
src_indexed_copy(__global real *dst, __global const real *src,
   __global const int *index, int src_offset)
{
        int id = get_global_id(ROW_DIM);
        int src_idx = src_offset + index[id];
        dst[id] = src[src_idx];
}

全局工作空间的工作项与索引数组中的索引一样多。

线性代码看起来像这样：

for (k = 0; k < n; k++) {
        dst[k] = src[m * column + index[k]];
}

其中复制了矩阵src中列列的所有索引元素。

这就是我回读缓冲区的方式（在评论中提到）：

rc = clEnqueueReadBuffer(ompctx->clctx.queue, c,
        CL_TRUE, 0, i * sizeof(real), &tmp[0],
        0, NULL, NULL);
if (rc != CL_SUCCESS) {
    log_error("omp", "[%d] readbuf() failed", rc);
    goto err;
}
log_info("omp", "c");
for (k = 0; k < i; k++) {
    log_info("omp", "%6.8f", tmp[k]);
}

Answer 1

主机代码中必定存在错误。请验证：

所有 OpenCL API调用的返回值，
缓冲区创建标志（您需要对目标缓冲区进行WRITE访问），
内核参数对应于内核代码。

两个缓冲区的索引副本

1 个答案: