mpi并行程序設(shè)計(jì)自學(xué)教程-資料下載頁(yè)
2025-01-12 20:03本頁(yè)面
【正文】 (rank=2) . . Recv()。 . . Recv()。 . . . . . Send()。 . . . . . . Send() . . . 問(wèn)題 :進(jìn)程 1和 2誰(shuí)先開(kāi)始發(fā)送消息 ?誰(shuí)先完成發(fā)送？ ? % 2022年 2月 42/55 運(yùn)行 greetings [dair@node01 ~]$ mpicc –o greeting [dair@node01 ~]$ mpirun np 4 greeting Greetings from process 1! Greetings from process 2! Greetings from process 3! [dair@node01 ~]$ 計(jì)算機(jī)打印字符 我們輸入的命令 2022年 2月 43/55 最基本的 MPI MPI調(diào)用借口的總數(shù)雖然龐大，但根據(jù)實(shí)際編寫(xiě) MPI的經(jīng)驗(yàn)，常用的 MPI調(diào)用的個(gè)數(shù)非常有限。上面介紹的是 6個(gè)最基本的 MPI函數(shù)。 1. MPI_Init(…)。 2. MPI_Comm_size(…)。 3. MPI_Comm_rank(…)。 4. MPI_Send(…)。 5. MPI_Recv(…)。 6. MPI_Finalize()。 MPI_Init(…)。 … 并行代碼 。 … MPI_Fainalize()。 只能有串行代碼 。 2022年 2月 44/55 現(xiàn)在您已經(jīng)能夠用 MPI進(jìn)行并行編程了 ! 2022年 2月 45/55 實(shí)例分析 :求 PI 2022年 2月 46/55 串行代碼 h=(double)n。 sum=。 for (i=1。 i=n。 i++) { x=h*((double)i – )。 sum += f(x)。 } pi=h*sum。 double f(double a) { return ((+a*a))。 } 2022年 2月 47/55 并行代碼 h=(double)n。 sum=。 for (i=myid+1。 i=n。 i+=numprocs) { x=h*((double)i – )。 sum += f(x)。 } mypi=h*sum。 MPI_Reduce(amp。mypi, amp。pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD)。 double f(double a) { return ((+a*a))。 } 2022年 2月 48/55 include include include double f( double )。 double f( double a )。 { return ( / ( + a*a))。 } 2022年 2月 49/55 int main( int argc, char *argv[]) { int done = 0, n, myid, numprocs, i。 double PI25DT = 。 double mypi, pi, h, sum, x。 double startwtime = , endwtime。 int namelen。 char processor_name[MPI_MAX_PROCESSOR_NAME]。 MPI_Init(amp。argc,amp。argv)。 MPI_Comm_size(MPI_COMM_WORLD,amp。numprocs)。 MPI_Comm_rank(MPI_COMM_WORLD,amp。myid)。 MPI_Get_processor_name(processor_name,amp。namelen)。 fprintf(stderr,Process %d on %s\n, myid, processor_name)。 2022年 2月 50/55 n = 0。 while (!done) { if (myid == 0) { if (n==0) n=100。 else n=0。 startwtime = MPI_Wtime()。 } MPI_Bcast(amp。n, 1, MPI_INT, 0, MPI_COMM_WORLD)。 2022年 2月 51/55 if (n == 0) done = 1。 else { h = / (double) n。 sum = 。 for (i = myid + 1。 i = n。 i += numprocs) { x = h * ((double)i )。 sum += f(x)。 } mypi = h * sum。 2022年 2月 52/55 MPI_Reduce(amp。mypi, amp。pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD)。 if (myid == 0) { printf(pi is approximately %.16f, Error is %.16f\n, pi, fabs(pi PI25DT))。 endwtime = MPI_Wtime()。 printf(wall clock time = %f\n, endwtimestartwtime)。 } } } MPI_Finalize()。 return 0。 }
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1