【发布时间】：2018-03-23 07:49:48
【问题描述】：
我正在尝试使用 openMP 来加快列表排名的并行版本。我的实现如下：
int ListRankingParallel(int *R1,int *S, int N)
{
int i;
int *Q = (int*)malloc(N * sizeof(int));
#pragma omp parallel for private(i)
for (i=0; i<N; i++){
if( S[i] != -1)R1[i] = 1;
else R1[i] = 0;
Q[i] = S[i];
}
#pragma omp parallel for private(i)
for(i=0; i<N; i++)
while (Q[i] != -1 & Q[Q[i]] != -1) {
R1[i] = R1[i] + R1[Q[i]];
Q[i] = Q[Q[i]];
}
free(Q);
return *R1;
}
我的榜单排名连载版是
int ListRankingSerial(int *R2,int *S, int N)
{
int temp;
int j,i;
for( i=0; i<N; i++){
j = 0;
temp = S[i];
while(S[i]!=-1)
{
j++;
S[i] = S[S[i]];
}
R2[i] = j;
S[i] = temp;
}
return *R2;
}
当我分别运行它们时，使用
get_walltime(&S1);
ListRankingParallel(R1,S,N);
get_walltime(&E1);
get_walltime(&S3);
ListRankingSerial(R3,S,N);
get_walltime(&E3);
如果我在我的 Mac 上运行我的代码，并行版本的运行速度明显快于串行版本。但是，如果我在另一个 linux 集群上运行它，并行版本比串行版本慢一倍。
在我的 Mac 上，我使用
编译我的代码gcc-7 -fopenmp <file name>.c
在集群上，使用
gcc -fopenmp <file name>.c
如果你想测试我的代码，请使用：
int main(){
int N = 1e+5;
int *S = (int*)malloc(N * sizeof(int));
int *R1 = (int*)malloc(N * sizeof(int));
int *R3 = (int*)malloc(N * sizeof(int));
double S1,S2,S3,E1,E2,E3;
int i;
for( i = 0; i < N; i++)
S[i] = i+1;
S[N-1] = -1;
get_walltime(&S1);
ListRankingParallel(R1,S,N);
get_walltime(&E1);
printf("%f\n",E1-S1);
get_walltime(&S3);
ListRankingSerial(R3,S,N);
get_walltime(&E3);
printf("%f\n",E3-S3);
}
谁能给我一些建议？谢谢！
【问题讨论】：
-
您通过以非同步方式访问同一变量来创建竞争条件。尝试在更新数组之前插入
#pragma omp atomic行。 -
@AlexQuilliam 我还是不太明白如何处理数组竞争条件。您能否更具体地说明我的代码的修改？谢谢！
-
列表排名是一个众所周知的难以并行化的问题。不要指望简单的 OpenMP 代码可以做到这一点。您确实需要一种避免并发写入同一数组元素的策略。使用 atomic 执行此操作是一种方法，但代价高昂，因此您不会看到任何性能提升。有很多关于 LR 的文献，请阅读。此外，现在您应该在
for循环中声明i而不是private(i)。 -
顺便说一句，在 C 语言中，
#include <stdlib.h>比 to cast the return value ofmalloc()and family 更好。 -
除了你的并行代码的竞争条件之外，永远不要对没有优化开关编译的代码进行性能测量。至少尝试使用
-O3和更好的-O3 -march=native -mtune=native