pragma omp for inside pragma omp master or single

Question

我在这里坐着一些尝试使孤儿工作的东西，并通过减少#pragma omp parallel的调用来减少开销。 我想要是一样的东西：

#pragma omp parallel default(none) shared(mat,mat2,f,max_iter,tol,N,conv) private(diff,k) 
{ 
#pragma omp master // I'm not against using #pragma omp single or whatever will work 
{ 
while(diff>tol) { 
    do_work(mat,mat2,f,N); 
    swap(mat,mat2); 
    if(!(k%100)) // Only test stop criteria every 100 iteration 
     diff = conv[k] = do_more_work(mat,mat2); 
    k++; 
} // end while 
} // end master 
} // end parallel 

的do_work取决于先前迭代这样的while循环是有运行顺序。 但是我想能够运行'do_work'平行的，所以它看起来是这样的：

void do_work(double *mat, double *mat2, double *f, int N) 
{ 
int i,j; 
double scale = 1/4.0; 
#pragma omp for schedule(runtime) // Just so I can test different settings without having to recompile 
for(i=0;i<N;i++) 
    for(j=0;j<N;j++) 
     mat[i*N+j] = scale*(mat2[(i+1)*N+j]+mat2[(i-1)*N+j] + ... + f[i*N+j]); 
} 

我希望这可以实现一些方法，我只是不知道如何。所以我能得到的任何帮助都非常感谢（如果你告诉我这是不可能的）。顺便说一句我正在打开mp 3.0，gcc编译器和sun studio编译器。

Answer 1

原始代码中的外部平行区域只包含一个连续片段（#pragma omp master），这在没有意义的情况下会导致纯串行执行（无并行）。由于do_work()取决于以前的迭代，但您想要并行运行它，您必须使用同步。 openmp工具是一个（显式或隐式）同步屏障。

例如（代码类似于你）：

#pragma omp parallel 
for(int j=0; diff>tol; ++j) // must be the same condition for each thread! 
#pragma omp for    // note: implicit synchronisation after for loop 
    for(int i=0; i<N; ++i) 
    work(j,i); 

注意，隐式同步确保没有线程进入下一j如果任何线程仍然在当前j工作。

替代

for(int j=0; diff>tol; ++j) 
#pragma omp parallel for 
    for(int i=0; i<N; ++i) 
    work(j,i); 

应不太有效的，因为它在每次迭代创建线程的一个新的团队，而不是仅仅同步。