CUDA Convex Hull 程序在大输入时崩溃

Question

我正在尝试在 CUDA 中并行实现 quickHull 算法（用于凸包）。它适用于 input_size

通过注释掉内核的行，我发现当我尝试访问数组元素并且我尝试访问的元素的索引没有超出范围时发生崩溃（双重检查）。以下是它崩溃的内核代码。

for(unsigned int i = old_setIndex; i < old_setIndex + old_setS[tid]; i++) 
{

    int pI = old_set[i];
    if(pI <= -1 || pI > pts.size())
    {               
        printf("Thread %d: i = %d, pI = %d\n", tid, i, pI);
        continue;
    }
    p = pts[pI];

    double d = distance(A,B,p);

    if(d > dist) {
        dist = d;
        furthestPoint = i;
        fpi = pI;
    }
}
//fpi = old_set[furthestPoint]; 
//printf("Thread %d: Furthestpoint = %d\n", tid, furthestPoint);

当我在 for 循环之后取消注释语句（数组访问和 printf）时，我的代码崩溃了。我无法解释该错误，因为 farthestPoint 始终在 old_set 数组大小的范围内。 Old_setS 存储每个线程可以操作的较小数组的大小。即使只是尝试打印 farthestPoint（最后一行）的值而没有上面的数组访问语句，它也会崩溃。

输入大小

请帮我找出崩溃的根源。

Answer 1

您的代码中没有越界内存访问（或者至少没有导致您所看到的症状）。

您的内核正在被显示驱动程序杀死，因为在显示 GPU 上执行需要花费太多时间。所有 CUDA 平台显示驱动程序都包含对 GPU 上的任何操作的时间限制。这样做是为了防止显示器在足够长的时间内冻结，以至于操作系统内核崩溃或用户恐慌并认为机器已经崩溃。在你使用的windows平台上，时间限制是2秒左右。

部分误导您认为问题的根源是数组寻址是代码注释使问题消失。但真正发生的是编译器优化的产物。当您注释掉全局内存写入时，编译器会识别出导致值被存储的计算未被使用，并从它发出的汇编代码中删除所有代码（谷歌“nvcc 死代码删除”以获取更多信息）。这样可以使代码运行得更快，并将其置于显示驱动程序的时间限制之下。

有关解决方法，请参阅this recent stackoverflow question and answer