使用 OpenMP 并行化递归的基本情况计算

Question

我正在尝试学习 OpenMP 的概念，但偶然发现了一个案例，我很难掌握如何使用这个库来解决这个问题。

假设我们有以下递归函数

// ...
void recurse(int tmp[], int p, const int size)
{
   if (p == size)
   {
      // Computationally heavy, should be executed in its own "thread"
      performTask(tmp); // Note: Only requires read access
   }
   else
   {
      for(int i = 0; i < size; i++)
      {
         // Alter tmp and continue recursion
         tmp[p] = i;
         recurse(tmp, p+1, size);
      }
   }
}
// ...
int main(int argc, char * argv[])
{
    int tmp[10];
    recurse(tmp, 0, 10);
    return 0;
}

如何在使用 OpenMP 在主线程中生成新结构的同时并行执行 performTask？

我知道有一种叫做“任务”的东西，我认为这就是我应该在这里使用的东西，但是我想出的所有东西根本没有任何性能提升。请指出正确的方向。

编辑：为了更好地解释，我将示例程序更具体。

Answer 1

下面的代码不能按原样工作，但希望它能为您指明正确的方向：

// ...
void recurse(int tmp[], int p, const int size)
{
   if (p == size)
   {
      // Computationally heavy, should be executed in its own "thread"
      // perform task using the thread pool
#pragma omp task     
      performTask(tmp); // Note: Only requires read access
   }
   else
   {
      for(int i = 0; i < size; i++)
      {
         // Alter tmp and continue recursion
         tmp[p] = i;
         recurse(tmp, p+1, size);
      }
   }
}
// ...
int main(int argc, char * argv[])
{    
    int tmp[10];
    // start threads
#pragma omp parallel
{
    // use single thread to construct `tmp` values
#pragma omp single nowait
    recurse(tmp, 0, 10);
}
    return 0;
}

代码基于Comparing Nested Parallel Regions and Tasking in OpenMP 3.0。