c++ - glCallList 上的 OpenMP 段错误的 GLFW3(以及 glfw 中的其他地方)

Question

我遇到了一些应该同时使用 glfw3 和 OpenMP 的代码的问题(顺便说一句，在 OSX 上使用 g++ 5.3.0)。这是我一直致力于编写的一些数学软件的简化版，到目前为止似乎进展顺利。代码本身庞大，甚至我的示例也非常庞大(我劫持了我的一个已经让 glfw 工作的旧教程项目)，所以我将附上代码所在的 git 存储库用于复制。我研究了多线程引起的其他 glfw3 问题，并提出了 this和 this .

glfwAppTutorial git repo!

齿轮应用程序中存在问题。当我在没有 OpenMP 的情况下编译它时，它可以工作。当我用 OpenMP 编译它并且只使用 1 个线程时，它可以工作。当我转向使用 2 个线程时，它会在 glfwPlatformPollEvents 调用上出现段错误。我在创建上下文时通过 std::this_thread::get_id() 调用检查了线程 ID，然后在我要更新它时检查了线程 ID，但是当我移动到多个 omp 线程时它似乎仍然出现段错误，但仅在我添加一组特定 命令之后。它比那更奇怪。这段开创性的代码是我开始做一些计算量大的事情的地方，我当然想并行化这些事情。代码如下所示(主要版本为 glfwAppTutorial/gears/src/gears.cpp::glloop，为清楚起见删除了无关的内容):

// gears.cpp excerpt
#pragma omp parallel
{
    int tid;
    double *fx, *fy, *fz;
    tid = omp_get_thread_num();

    fx = frc_ + (3*tid)*10;
    for (int i = 0; i < 3*10; ++i) {
        fx[i] = 0.0; // BOOM!!!!!!!!!, including this segfaults, but not here...
    }
    fy = frc_ + (3*tid+1)*10;
    fz = frc_ + (3*tid+2)*10;

    std::cout << "Running on thread: " << tid << std::endl;

    #pragma omp for schedule(runtime) nowait
    for (int i = 0; i < 10; ++i) {
        // simulate work
        std::mt19937_64 eng{std::random_device{}()};  // or seed however you want
        std::uniform_int_distribution<> dist{10, 100};
        std::this_thread::sleep_for(std::chrono::milliseconds{dist(eng)});
    }
    // synchronize if we need to
    #pragma omp barrier

} // pragma omp parallel

前面定义 frc_ 的地方(并归零为所有 0.0 double ):

// gears.hpp excerpt
nthreads = omp_get_num_threads();
frc_ = new double[3*10*nthreads];
for (int i = 0; i < 3*10*nthreads; ++i) {
    frc_[i] = 0.0;
}

当然，在这个例子中，一切都在(glfwApp.cpp)下运行:

//glfwApp.cpp excerpt
void glfwApp::glfwloop() {
    while(!glfwWindowShouldClose(_window)) {
        this->glloop();

        glfwSwapBuffers(_window);
        glfwPollEvents();
    }
}

导致段错误的行是 fx[i] = 0.0，但该行不是实际的段错误源。 gleCallList 或 glfwPlatformPoll 中的程序段错误。那么有人知道会发生什么吗？不使用 glfw3 的纯计算版本很好，我已经使用 valgrind(或在 OSX 上，Instruments)来检查可能存在的问题。我所做的只是一些应该相当简单的指针算法，但由于某些原因不是。我将尝试创建一个更更简单的示例，因为我意识到这可能无法遵循，但几天来它一直让我抓狂。

这是来自 OSX 的堆栈跟踪，因为 Xcode 和 gdb 现在很愚蠢

Exception Type:        EXC_BAD_ACCESS (SIGSEGV)
Exception Codes:       KERN_INVALID_ADDRESS at 0x0000000000000110

VM Regions Near 0x110:
--> 
    __TEXT                 0000000106019000-0000000106020000 [   28K] r-x/rwx SM=COW  /Users/USER/*

Thread 0 Crashed:: Dispatch queue: com.apple.main-thread
0   libGPUSupportMercury.dylib      0x00007fff8d03f57c gpusLoadCurrentVertexArray + 132
1   com.apple.AMDRadeonX4000GLDriver    0x000000010b66c695 gldUpdateDispatch + 2819
2   GLEngine                        0x00007fff8b69bad4 gleDoDrawDispatchCore + 533
3   GLEngine                        0x00007fff8b63e636 gleSetupAndDrawArraysOrElementsOutOfLine_ListExec + 886
4   GLEngine                        0x00007fff8b5c8a6a gleCallList + 170
5   gears                           0x000000010601bc25 gears::draw() + 199
6   gears                           0x000000010601c285 gears::glloop() + 33
7   gears                           0x000000010601dd2a glfwApp::glfwloop() + 60
8   gears                           0x000000010601dcea glfwApp::start() + 24
9   gears                           0x000000010601d04a main + 171
10  libdyld.dylib                   0x00007fff967475ad start + 1

Answer 1

除非从并行区域的动态范围内调用 gears 的构造函数，否则 omp_get_num_threads() 返回 1，因此 frc_ 是空间不够大，无法容纳多个线程的数据，并且堆会被第一段中显示的代码损坏。尝试替换 omp_get_num_threads() here使用 omp_get_max_threads() 看看是否有帮助。此外，请确保 omp_set_dynamic(0); 存在于您的代码中的某处，并且在其间没有调用 omp_set_num_threads()，否则团队规模可能与并行区域到另一个区域，并使用 omp_get_num_threads() 的值。