技术沙龙-OpenMP并行技术-China-VO.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 实例1：蒙特卡罗算法 利用蒙特卡罗算法计算半径为1单元的 球体体积：(下面为相应的串行代码) #includestdio.h #includestdlib.h #includetime.h int main() { long long max long long i,count=0; double x,y,z,bulk,start_time,end_time; time_t t; start_time=clock(); // 产生以当前时间开始的随机种子 srand((unsigned) time(t)); for(i=0;imax;i++) { x=rand(); x=x/32767; y=rand(); y=y/32767; z=rand(); z=z/32767; if((x*x+y*y+z*z)=1) count++; } bulk=8*(double(count)/max); end_time= clock(); printf(”Sphere bulk is %f \n“, bulk); printf(“Running time is %f \n“, end_time-start_time); return 0; } * * OpenMP Parallel for with a reduction #includestdio.h #includestdlib.h #includetime.h #include omp.h int main() { long long max long long i,count=0; double x,y,z,bulk,start_time,end_time; time_t t; start_time=omp_get_wtime(); //产生以当前时间开始的随机种子 srand((unsigned) time(t)); #pragma omp parallel for private(x,y,z) reduction(+:count) for(i=0;imax;i++) {x=rand(); x=x/32767; y=rand(); y=y/32767; z=rand(); z=z/32767; if((x*x+y*y+z*z)=1) count++; } bulk=8*(double(count)/max); end_time=omp_get_wtime(); printf(”Sphere bulk is %f \n“, bulk); printf(“Running time is %f \n“, end_time-start_time); return 0; } * * 作业 1. 循环是否可以直接并行 do i=2,n a(i)=2*a(i-1) end do ix = base do i=1,n a(ix) = a(ix)*b(i) ix = ix + stride end do do i=1,n b(i)= (a(i)-a(i-1))*0.5 end do * * 2.上机调试培训内容中提供的求解pi的各OpenMP并行程序。 3.编写一个矩阵-向量相乘的 OpenMP并行程序。 4.分析critical、atomic、lock、flush的用法 谢谢大家！ * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 调度子句SCHEDULE例图 * * 数据竞争问题 下面的循环无法正确执行： #pragma omp parallel for for(k=0;k100;k++) { x=