第03章 消息传递编程.pptVIP

二、ＭＰＩ MPI集合通信 广播 MPI_Bcast(Address,Count,Datatype,Root,Comm) 汇集 MPI_Gather(SendAddr,SendCount,SendDatatype,RecvAddr, RecvCount,RecvDatatype,Root,Comm) 散射 MPI_Scatter(SendAddr,SendCount,SendDatatype,RecvAddr, RecvCount,RecvDatatype,Root,Comm) 全交换 MPI_Alltoall(SendAddr,SendCount,SendDatatype,RecvAddr, RecvCount,RecvDatatype,Comm) 聚集 MPI_Reduce(SendAddr,RecvAddr,Count,Datatype,Op,Root,Comm) MPI_Scan(SendAddr,RecvAddr,Count,Datatype,Op,Comm) 路障 MPI_Barrier(Comm) 设有三个进程，执行如下代码： int i, j, my_rank, group_size, A[3], B[3], tag=1, root=0; MPI_Comm comm; MPI_Init( argc, argv ); comm = MPI_COMM_WORLD; MPI_Comm_rank( comm, my_rank ); MPI_Comm_size( comm, group_size ); for (i=0; i3; i++) A[i] = B[i] = my_rank*group_size+i; MPI_Bcast( A, 3, MPI_INT, root, comm ); MPI_Scatter( A, 1, MPI_INT, B, 1, MPI_INT, root, comm ); MPI_Scan( A, B, 1, MPI_INT, MPI_SUM, comm ); B A 1 1 2 0 B 2 0 A P1 4 4 5 3 5 3 P2 7 7 8 6 B 8 6 A P3 B A 1 1 2 0 B 2 0 A P1 4 1 5 3 2 0 P2 7 1 8 6 B 2 0 A P3 B A 1 1 2 0 B 2 0 A P1 4 4 5 3 5 3 P2 7 7 8 9 B 8 6 A P3 B A 1 1 2 0 B 2 0 A P1 4 4 5 1 5 3 P2 7 7 8 2 B 8 6 A P3 二、ＭＰＩ MPI集合通信 通信子中所有进程必须调用集合通信例程。若某些成员调用集合例程而其余的没有调用，该代码是错误的，其行为是未定义的。 集合通信采用标准锁定通信方式，即进程完成集合操作中的参与动作后便立刻从集合例程中返回，此时发送缓冲区为安全可重用的，而其他进程可能还未开始相应的集合通信。 集合通信的同步性依赖于具体的实现，因此要求用户负责代码的正确性。 在集合通信中的Count和Datatype必须是兼容的。 集合例程没有消息标志参数，消息信封由通信子和源/目的进程说明构成。 #include “mpi.h” int foo(int i){ ... } main( int argc, char *argv[] ) { int i, tmp, sum=0, group_size, my_rank, N; MPI_Status status; MPI_Init( argc, argv ); MPI_Comm_size( MPI_COMM_WORLD, group_size ); MPI_Comm_rank( MPI_COMM_WORLD, my_rank ); if (my_rank==0) { printf( “N=“ ); scanf( “%d”, N ); for (i=1; igroup_size; i++) S1: MPI_Send( N, 1, MPI_INT, i, i, MPI_COMM_WORLD; ); for (i=my_rank; iN; i+=group_size) sum += foo(i); for (i=1; igroup_size; i++) { S2: MPI_Recv( tmp, 1, MPI_INT, i, i, MPI_COMM_WORLD, status ); sum += tmp; } printf (“\nThe result=%d\n”, sum ); } else { S3: MPI_Recv( N, 1, MPI_INT, 0, my_r