不知道全局大小的 OpenCL 1D 范围循环

Question

我想知道如何迭代具有任意数量工作项的循环（每个组无关紧要）

我有 3 个 arrays，其中一个是二维的（矩阵）。 第一个数组包含一组整数。 矩阵填充了另一组（重复和随机）整数。 第三个只是存储结果。 我需要从矩阵中的第一个数组中搜索最远的一对数字出现次数。 总结一下：

• A：带有随机数的矩阵
• num：包含要在A中搜索的数字的数组
• d：每个数字与num之间距离最大的数组

该算法很简单（因为我不需要优化它），我只比较计算出的曼哈顿距离并保留最大值。

为了简单起见，它执行以下操作（类似 C 的伪代码）：

for(number in num){
  maxDistance = 0
  for(row in A){
    for(column in A){
      //calculateDistance is a function to another nested loop like this
      //it returns the max found distance if it is, and 0 otherwise
      currentDistance = calculateDistance(row, column, max)
      if(currentDistance > maxDistance){
        maxDistance = currentDistance
      }
    }
  }
}

如您所见，迭代之间没有相关数据。 我试图为每个工作项分配矩阵 A 的一部分，但仍然不能说服我。

重要提示：kernel 必须仅使用一个维度来解决问题。

有任何想法吗？ 如何使用全局 ID 一次进行多个搜索？

编辑：

我添加了代码以消除任何疑问。 这是 kernel：

__kernel void maxDistances(int N, __constant int *A, int n, __constant int *numbers, __global int *distances)
{
    
    //N is matrix row and col size
    //A the matrix
    //n the total count of numbers to be searched
    //numbers is the array containing the numbers
    //distances is the array containing the computed distances
    

    size_t id = get_global_id(0);   
    

    int slice = (N*N)/get_global_size(0);

    for(int idx_num = 0; idx_num < n; idx_num++)
    {
        int number = numbers[idx_num];
        int currentDistance = 0;
        int maxDistance = 0;


        for(int c = id*slice; c < (id+1)*slice; c++)
        {
            int i = c/N;
            int j = c%N;
            if(*CELL(A,N,i,j) == number){
                coord_t coords;
                coords.i = i;
                coords.j = j;

                //bestDistance is a function with 2 nested loop iterating over
                //rows and column to retrieve the farest pair of the number
                currentDistance = bestDistance(N,A,coords,number, maxDistance);
                if(currentDistance > maxDistance)
                {
                    maxDistance = currentDistance;
                }
            }
        }

        distances[idx_num] = maxDistance;
    }
}

Answer 1

这个答案可能被认为是不完整的，不过，我将发布它以结束这个问题。

我的问题不是代码，kernel（或那个算法），而是机器。 上面的代码是正确的并且可以完美运行。 在我在另一台机器上尝试我的程序后，它执行并计算了解决方案，完全没有问题。

因此，简而言之，问题出在 OpenCL 设备或很可能是主机库。