mpi分布內(nèi)存并行程序開發(fā)-資料下載頁
2025-07-17 15:46本頁面
【正文】 me,amp。namelen)。 fprint(stdout,”Process %d of %d on % s\n”,myid,numprocs, processor_name)。 n=0。 if (myid==0) { printf(“Please give N=”)。 scanf(amp。n)。 startwtime=MPI_Wtime()。 for (j=1。jnumprocs。j++) { MPI_Send(amp。n,1,MPI_INT,j,99,MPI_COMM_WORLD)。 } } 41 else { MPI_Recv(amp。n,1,MPI_INT,MPI_ANY_SOURCE,99,MPI_COMM_WORLD,amp。status)。 } h=(double) n。 sum=。 for(i=myid+1。i=n。i+=numprocs) /* 每一個進(jìn)程計(jì)算一部分矩形的面積，若進(jìn)程總數(shù) numprocs為 4， 將 01區(qū)間劃分為 100個矩形，則各個進(jìn)程分別計(jì)算矩形塊 0進(jìn)程 1， 5， 9， 13， …… ， 97 1進(jìn)程 2， 6， 10， 14， …… ， 98 2進(jìn)程 3， 7， 11， 15， …… ， 99 3進(jìn)程 4， 8， 12， 16， …… ， 100 */ 42 { x=h*((double))。 sum+=f(x)。 } mypi=h*sum。 /*各進(jìn)程并行計(jì)算得到的部分和 */ /*將部分和累加得到所有矩形的面積，該面積和即為近似 PI值 */ if (myid != 0) MPI_Send(amp。mypi,1,MPI_DOUBLE,0,myid,MPI_COMM_WORLD) ； else { pi=。 pi=pi+mypi。 43 for (j=1。jnumprocs。j++) { MPI_Recv(amp。mypi,1,MPI_DOUBLE,MPI_ANY_SOURCE, MPI_ANY_TAG,MPI_COMM_WORLD,amp。status)。 pi=pi+mypi。 } printf(“pi is approximately %.16f,Error is %.16f\n”, pi,fabs(piPI25DT))。 endwtime=MPI_Wtime()。 printf(“wall clock time=% f\n”,endwtimestartwtime)。 fflush(stdout)。 } MPI_Finalize()。 } 44 謝謝！
高考資料相關(guān)推薦
鄂ICP備17016276號-1