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

Question

我正在尝试在 CUDA 中并行实现 quickHull 算法(用于凸包)。它适用于 input_size <= 100 万。当我尝试 1000 万点时，程序崩溃了。我的显卡大小是 1982 MB，算法中的所有数据结构总共需要不超过 600 MB 的输入大小，这不到可用空间的 50%。

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

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)时，我的代码崩溃了。我无法解释这个错误，因为 furthestPoint 总是在 old_set 数组大小的范围内。 Old_setS 存储每个线程可以操作的较小数组的大小。即使只是尝试打印 furthestPoint(最后一行)的值而没有上面的数组访问语句，它也会崩溃。

对于输入大小 <= 100 万，上面的代码没有问题。如果是 1000 万，我是否会溢出设备中的某些缓冲区？

请帮助我找到崩溃的根源。

Answer 1

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

发生的事情是您的内核被显示驱动程序杀死，因为它在您的显示 GPU 上执行花费了太多时间。所有 CUDA 平台显示驱动程序都包含对 GPU 上的任何操作的时间限制。这样做是为了防止显示卡住足够长的时间，以免操作系统内核崩溃或用户崩溃并认为机器崩溃了。在你使用的windows平台上，时间限制在2秒左右。

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

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