多维嵌套的OpenMP循环

Question

在OpenMP中并行化多维尴尬平行循环的正确方法是什么？尺寸的数量在编译时是已知的，但哪些尺寸很大并不是。他们中的任何一个可能是一个，两个或一百万。当然我不想ňomp parallel的为N维环...

思考：

• 问题的概念很简单。只有最外面的'大'循环需要并行化，但循环维度在编译时是未知的，并且可能会改变。

• 将动态设置omp_set_num_threads(1)和#pragma omp for schedule(static, huge_number)使某些循环并行化为无操作？这会产生不希望的副作用/开销吗？感觉像一个kludge。

• OpenMP Specification（2.10，A.38，A.39）讲述了符合和不符合嵌套并行性之间的区别，但没有提出解决这个问题的最佳方法。

• 重新排序循环是可能的，但可能会导致大量缓存未命中。展开是可能的，但不是微不足道的。有另一种方法吗？

这里是想什么我并行：

for(i0=0; i0<n[0]; i0++) { 
    for(i1=0; i1<n[1]; i1++) { 
    ... 
     for(iN=0; iN<n[N]; iN++) { 
     <embarrasingly parallel operations> 
     } 
    ... 
    } 
} 

谢谢！

Answer 1

collapse指令可能是你要找的，如描述here。这基本上会形成一个单一的循环，然后这个循环被并行化，并被设计用于这些种类的情况。所以你会这样做：

#pragma omp parallel for collapse(N) 
for(int i0=0; i0<n[0]; i0++) { 
    for(int i1=0; i1<n[1]; i1++) { 
    ... 
     for(int iN=0; iN<n[N]; iN++) { 
     <embarrasingly parallel operations> 
     } 
    ... 
    } 
} 

并且都设置好了。