程序设计与技术实验报告.doc

程序设计与技术实验报告 一、实验题目 用MPI技术计算π的值。 二、实验目的 熟悉mpi编程，加深对mpi编程的理解。 三、实验内容 搭建mpi的编程环境，并计算出π的值。 四、源代码 #include mpi.h #include stdio.h #include math.h double f(a) double a; { return (4.0/(1.0+a*a)); } int main(argc,argv) int argc; char *argv[]; { int done=0,n=100,myid,numprocs,i; double PI25DT=3.141592653589793238462643; double mypi,pi,h,sum,x,a,startwtime,endwtime; MPI_Init(argc,argv); MPI_Comm_size(MPI_COMM_WORLD,numprocs); MPI_Comm_rank(MPI_COMM_WORLD,myid); if (myid == 0) startwtime = MPI_Wtime(); h=1.0/(double) n; sum=0.0; for (i=myid+1;i=n;i+=numprocs) { x=h*((double)i-0.5); sum+=f(x); } mypi=h*sum; MPI_Reduce(mypi,pi,1,MPI_DOUBLE,MPI_SUM,0,MPI_COMM_WORLD); if (myid==0) { printf(pi is %.16f, Error is %.16f\n, pi, fabs(pi - PI25DT)); endwtime=MPI_Wtime(); printf(wall clock time = %f\n, endwtime-startwtime); } MPI_Finalize(); } 五、实验结果截图 六．实验心得 这个实验对于我这种正在学C的同学太难了，以前从未使用过这个软件，这是第一次使用，问了很多人才写成这样。编写一段程序不是一朝一夕就能完成的，它是靠平时一点一滴的积累，靠自己的耐心，永不放弃的纠正错误，不能嫌它过于繁琐。感谢赵老师的耐心指点，我发现需要我学习和改进的地方还很多，我会更加努力！