关于双核并行后的错误

[memoryjs]

我用双核并行马太棋盘，用intel编译器后，运行速度是快了，但是出来的结果是错误的，如果不用并行，那么时间会长，但是答案是对的 我用的是 #pragma omp task #pragma omp parallel #pragma omp single 这三个命令并行的 是不是我代码用错 还是什么别的？而且我并行别的程序也是用这三个指令都有错误 ...

[denghui0815]

估计你有需要同步的数据访问没有做同步

[fanbin23]

Thread Checker
Thread Inspector
Sun Studio Thread Analyzer
贴代码
....

[无天]

代码?

[memoryjs]

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>
void read_file(int *m,int *n,int *p);
int write_file(int *m,int *n);
void draw_chessboard(int *m,int *n);//根据m，n画棋盘
void find_location(int *m,int *n,int *p,int *i,int*j);//根据数值，找出坐标
void move_chess(int i,int j,int k,int m,int n);
int c[32][32],a[32][32];//c是画的棋盘，a是用来记录每一步
int ok=0;//0为不成功，1为成功
main()
{
int m,n,p,i,j;
double d;
time_t begin,end;
a[0][0]=0;
read_file(&m,&n,&p);
printf("%d %d\n",m,n);
printf("%d\n",p);
draw_chessboard(&m,&n);
find_location(&m,&n,&p,&i,&j);
a[i][j]=1;//初始化，使得得到的坐标为第一步
begin=clock();
move_chess(i,j,2,m,n);

for(i=0;i<m;i++)
for(j=0;j<n;j++)
{ printf(" %d ",a[i][j]);
if(j==n-1)
printf("\n");
}
end=clock();
d=(double)(end-begin)/CLOCKS_PER_SEC;
printf("your run time is:%f\n",d);
}
void move_chess(int i,int j,int k,int m,int n)
{
int x,g,h,p=-1;//g记录i的变化，h记录j的变化，p每次运行都要加1，所以定义为负数，用来访问数组
int direction[8][2]={{2,1},{1,2},{-1,2},{-2,1},{-2,-1},{-1,-2},{1,-2},{2,-1}};//定义8个方向
do{
p++;
ok=0;
g=i+direction[p][0];
h=j+direction[p][1];
#pragma omp parallel
#pragma omp single
if(g>=0&&g<m&&h>=0&&h<n&&a[g][h]==0){
a[g][h]=k;
if(k<m*n){
#pragma omp task
move_chess(g,h,k+1,m,n);
if(ok==0)
a[g][h]=0;//如果八个方向都不行，从函数退出，并且使原来尝试的一步清零

}
else
for(x=0;x<8;x++)

if((g+direction[x][0])<m&&(g+direction[x][0])>=0&&(h+direction[x][1])<n&&(h+direction[x][1])>=0&&a[g+direction[x][0]][h+direction[x][1]]==1) {
ok=1;break;
}
if(ok==0) a[g][h]=0;

}


}while(ok==0&&p!=7);

}
void find_location(int *m,int *n,int *p,int *i,int *j)
{
int k=0;
for((*i)=0;(*i)<*m;(*i)++)
{
for((*j)=0;(*j)<*n;(*j)++)
if(c[*i][*j]==*p){
k=1;
break;
}
if(k==1)
break;
}
// printf("%d %d\n",*i,*j);

}
void read_file(int *m,int *n,int *p)
{
FILE *fp;
if((fp=fopen("hamilton.in","r"))==NULL){
printf("Opent the hamilton.in file fail!\n");
exit(0);
}
fscanf(fp,"%d",m);
fscanf(fp,"%d",n);
fscanf(fp,"%d",p);
}
void draw_chessboard(int *m,int *n)
{
int i,j,s,count=0;
s=(*m)*(*n);
for(i=0;i<*m;i++)
for(j=*n-1;j>=0;j--)

{
c[i][j]=s-1;
s--;
}
/* for(i=0;i<*m;i++)
for(j=0;j<*n;j++)
{
printf(" %d ",c[i][j]);
count++;
if(count==8){
printf("\n");
count=0;
}
}*/
}
帮忙看看 哪里错了

[memoryjs]

我是用文件读如的 文件名为hamilton.in 没有tet后缀名 里面的内容是 6 6 1 即是一个6*6上的一个棋盘上 在1这个位置为起点 数字大了 运行时间太长 我用单核这个时间为5秒左右 用双核并行就错了 你可以试试

[memoryjs]

我贴了代码了 你帮忙看看

[fanbin23]

没仔细看，但有这么几个建议：
1. 按照你对程序的理解把parallel region中所有用到的变量显示的指明他们的data sharing clause，例如g,h,k,m,n和数组a。如果你用linux，或许你可以试一下Sun Studio OpenMP compiler的autoscoping工具
2. move_chess似乎是个递归函数，或许你应该考虑把parallel single结构拿到递归外面去，只在递归函数里面使用task

[memoryjs]

虽然不是很懂 但是我还是试试看 那么并行程序就是用openmp的命令？还有TBB？

[loadbalancing]

结果是错误的，则并行没有意义了。首先得保证结果是对的。

[memoryjs]

这个道理我也知道啊 但是我不并行是对的 并行就错了 我都不知道是哪里错了 要不你也帮忙看看

[fanbin23]

看起来你还没理解什么是并行和OpenMP……找一些基本的材料看看吧，具体是什么材料……我现在也想不起来该推荐什么
OpenMP和TBB都是编写并行程序的方式而已

[memoryjs]

的确不太理解 我找了一些资料 感觉帮助不大 用了之后速度没有提高多少 要不你推荐一下？

[memoryjs]

就是弄不懂openMp 和TBB怎么用

[fanbin23]

看材料其实没啥用，你可以就从这个程序入手，想清楚为什么要这么写，为什么这么用那几个OpenMP directives。弄明白了基本上OpenMP你就懂了一大半了

[fanbin23]

OpenMP官网上有个tutorial链接，另外随便搜搜能找到一些slides，稍微有点概念就行了，看多了也没啥用
https://computing.llnl.gov/tutorials/openMP/

TBB就不知道了