核心在嵌套循环中是否在OpenMP中分区？

Question

让我们假设我们有这样的代码：

for(size_t i=0; i<5; i++){ 
    //do something expensive here 
    for(size_t i=0; i<5; i++){ 
    //do somehting expensive here too 
    } 
} 

让我们假设我们有至少25个物理内核。我们如何利用所有这些？ A collapse将是不可能的，因为第一个昂贵的计算。在这种情况下，最有效的解决方案是什么？

Answer 1

只解耦环路。如果您当前的代码是

#pragma omp parallel for 
for(int i = 0; i < 5; ++i) 
{ 
    expensive_1(i); 
    for(int j = 0; j < 5; ++j) 
    { 
    expensive_2(j); 
    } 
} 

然后使它成这样，可能节省的expensive_1的结果，如果需要，他们在expensive_2。

#pragma omp parallel for 
for(int i = 0; i < 5; ++i) 
{ 
    expensive_1(i); 
} 

#pragma omp parallel for collapse(2) 
for(int i = 0; i < 5; ++i) 
{ 
    for(int j = 0; j < 5; ++j) 
    { 
    expensive_2(j); 
    } 
}