Function main

code/Chapter9/9-32.c:37–88

        }
}
int main()
{
	isync[0]=1;
	for(i=0;i<N1;i++)
	{
		for(j=0;j<N2;j=j+4)
			a[i][j]=3;
		for(j=1;j<N2;j=j+4)
			a[i][j]=5;
		for(j=2;j<N2;j=j+4)
			a[i][j]=4;
		for(j=3;j<N2;j=j+4)
			a[i][j]=8;
		for(j=4;j<N2;j=j+5)
			a[i][j]=7;
	}
	double start_time,end_time,used_time;
	start_time=omp_get_wtime();
        #pragma omp parallel default(shared) private(i,j,m) shared(a) num_threads(3)
        {
                mthreadnum=1;
                mthreadnum=omp_get_num_threads()+1;
				iam=1;
                iam=omp_get_thread_num()+1;
		    	isync[iam]=0;
                #pragma omp barrier
                int b=9;
                for(int j=1;j<N2-1;j=j+b)
{
                        sync_left();
                        #pragma omp for schedule(static) nowait
                        for(int i=1;i<N1-1;i++)
                        {
                                for(int m=j;m<min(j+b,N2-1);m++)
							{
                                        a[i][m]=0.25*(a[i-1][m]+a[i][m-1]+a[i+1][m]+a[i][m+1]);
						}
                        }
                        sync_right();
                }
        }
	end_time=omp_get_wtime();
    used_time=end_time-start_time;
	
	double sum=0;
	for(i=0;i<N1;i++)
	{
		for(j=0;j<N2;j++)
			sum+=a[i][j];
	}
	printf("sum=%lf,used_time=%lf seconds\n",sum,used_time);
}