OpenCL 一维数组大数乘法

Question

我正在研究将两个整数数组相乘以获得大数。我基本上乘以大约 20 位数字，但对于更多数字，它表现不稳定。我在 CPU 上有一些代码，用于测试的内核和内核在每次工作后给我不同的数字。问题可能是障碍或互斥锁之类的东西，但我对此不是很感兴趣。我该如何解决这种不稳定性？这是代码

kernel void multiply(global int* A,
                 const int M,
                 global int* B,
                 const int N,
                 global int* C){
const int globalRow = get_global_id(0); // Row ID of C (0..M)
int globalCol,i;

// Compute a single element (loop over K)
for (globalCol=0; globalCol<N; globalCol++) {
    int val=A[globalRow]*B[globalCol];

    printf("Row is %d , Col is %d \n",globalRow,globalCol);
    //C[globalCol + globalRow +1]+=val/10;
    C[globalCol + globalRow]+=val%10;

    C[globalCol+1+globalRow]+=val/10;

}

int flag=1;
while (flag) {
    flag=0;
    for (i=M+N-1   ; i>=0 ; i--) {
        if (C[i]>9) {
            C[i+1]+=C[i]/10;
            C[i]=C[i]%10;
            flag=1;
        }
    }
}

}

代码上的标志是 shifting.It 听起来像是访问元素的问题，但我不知道如何解决。

Answer 1

您有竞争条件。相邻线程正在访问相同的内存位置，同时，在对 C 数组执行 += 的两行代码中。由于它们遵循特定模式，您可以在两个 += 操作之间使用 barrier(CLK_GLOBAL_MEM_FENCE) 来消除竞争条件。我认为你在下面的进位代码中有类似的东西。请小心，因为 工作组中的所有线程都必须达到障碍，否则您的内核将挂起或崩溃。

OpenCL 一维数组大数乘法

OpenCL 1D Array Big Number Multiply

c

gpgpu

opencl