c++ - 为什么在使用嵌套的 OpenMP pragma 时 c++11 线程变得不可连接？

Question

以下代码应该非常简单，但在尝试使用嵌套的 OpenMP 代码在线程上执行 .join() 时，似乎最终陷入了挂起状态。使用 GCC 编译器 4.7.2 x64 和来自 http://sourceforge.net/projects/mingwbuilds 的 pthreads使用 g++ threadexample.cpp -Wall -std=c++11 -fopenmp -o threads

// threadexample.cpp
#include <iostream>
#include <thread>
#include <omp.h>

using namespace std;

void hello(int a) {

    #pragma omp parallel for
        for (int i=0;i<5;++i) {
            #pragma omp critical
            cout << "Hello from " << a << "! " << "OMP thread iter " << i << endl;
        }

    cout << "About to return from hello function" << endl;
}

int main (int argc, char ** argv) {

    thread t1(hello, 1); //fork
    cout << "t1 away!" << endl;
    thread t2(hello, 2);
    cout << "t2 away!" << endl;

    t1.join(); //join
    cout << "thread 1 joined" << endl;
    t2.join();
    cout << "thread 2 joined" << endl;

    return 0;
}

Answer 1

混合使用 OpenMP 和任何其他线程库(pthreads、Win32 线程等)可能不是一个好主意。 OpenMP 运行时的编写可能假设它完全控制线程并且可能不支持同时运行的parallel 区域(例如，它可能使用像信号量这样的全局变量来控制线程池) .

实现此目的的更好的纯 OpenMP 方法是:

#include <iostream>
#include <omp.h>

using namespace std;

void hello(int a) {

    #pragma omp parallel for
    for (int i=0;i<5;++i) {
        #pragma omp critical
        cout << "Hello from " << a << "! " << "OMP thread iter " << i << endl;
    }

    cout << "About to return from hello function" << endl;
}

int main (int argc, char ** argv) {

    omp_set_nested(1);

    #pragma omp parallel sections num_threads(2)
    {
       #pragma omp section
       {
           hello(1);
       }
       #pragma omp section
       {
           hello(2);
       }
    }

    return 0;
}

需要调用 omp_set_nested() 以启用默认情况下禁用的嵌套并行性。