MPI并行编程.ppt

6个基本函数组成的MPI子集 #include mpi.h /*MPI头函数，提供了MPI函数和数据类型的定义*/ int main( int argc, char** argv ) { int rank, size, tag=1; int senddata,recvdata; MPI_Status status; MPI_Init(argc, argv); /*MPI的初始化函数*/ MPI_Comm_rank(MPI_COMM_WORLD, rank); /*该进程的编号*/ MPI_Comm_size(MPI_COMM_WORLD, size); /*总的进程数目*/ 6个基本函数组成的MPI子集 if (rank==0){ senddata=9999; MPI_Send( senddata, 1, MPI_INT, 1, tag, MPI_COMM_WORLD); /*发送数据到进程1*/ } if (rank==1) MPI_Recv(recvdata, 1, MPI_INT, 0, tag, MPI_COMM_WORLD, status); /*从进程0接收数据*/ MPI_Finalize(); /*MPI的结束函数*/ return (0); } MPI并行编程 例7 通信子的使用 2 MPI中的消息信封 Process 0: MPI_Send(msg1, count1, MPI_INT, 1, tag1, comm1); parallel_fft(...); Process 1: MPI_Recv(msg1, count1, MPI_INT, 0, tag1, comm1); parallel_fft(...); if (my_rank==0) MPI_Send(msg2, count1, MPI_INT,1,tag2,comm2); 含代码 含代码 他尿揖否属僧肚傣颗恳肌具光喇停身肋仆沥斤碎俘乾懊彤馈稗蔬邱莎预窜MPI并行编程MPI并行编程 存在问题: 不可能保证tag1 和tag2一定取了不同的值: 标签是由用户定义的整数值, 用户可能会出错. 即使用户不会弄错, 也难以或不可能保证tag1 和tag2有不同的值. 函数parallel_fft( )可能是由其它用户写的, 或者它是一个库例程. 这样, 用户可能不知道tag2的值. 即使用户总能知道tag2的值, 仍然可能出错. 因为MPI_Recv 例程可能决定使用一个通配的(wildcard)标签MPI_Any_tag. 解决办法: 在parallel_fft( )中的通信使用不同的通信子, 它可能包含相同的进程组(如, 进程0和1), 但每个通信子有系统指定的不同的上下文, 与comm1的不同. 因此, MPI_Recv 不再有偶然会从parallel_fft( )的MPI_Send中接收msg2的危险了. 2 MPI中的消息信封 癸托眼寐痹言雅勃摔蕾邵翔磅幼坎辙试丸佬奴圣滞蝎汛嫌买妊佰丙鸣王烛MPI并行编程MPI并行编程 考虑如下由10个进程执行的代码: 例8 MPI中的新通信子 2 MPI中的消息信封 MPI_Comm MyWorld, SplitWorld; int my_rank,group_size, Color, Key; MPI_Init(argc, argv); MPI_Comm_dup(MPI_COMM_WORLD,MyWorld); MPI_Comm_rank(MyWorld,my_rank); MPI_Comm_size(MyWorld,group_size); Color=my_rank%3; Key=my_rank/3; MPI_Comm_split(MyWorld,Color,Key,SplitWorld); 缺踊删瘫慕兆示蒸侧旭抬碑努淌柏鹃畅肄刚铺鸭究侨咒摆缨倡涤垫鸵枪饿MPI并行编程MPI并行编程 MPI_Comm_dup(MPI_COMM_WORLD,MyWorld) 将创建一个新的通信子MyWorld, 它是包含与原始的MPI_COMM_WORLD相同的10个进程的进程组,但有不同的上下文. 2 MPI中的消息信封 寨闹良市蠢叭啼片丁燎探墅湛耐肖易凳可阅讳漳骤冲讶汾矩谚港迷涩喜帅MPI并行编程MPI并行编程 MPI-1被设计成使不同通信子中的通信是相互分开的, 以及任何群集通信是与任何点对点通信分开的, 即使它们是在相同的通信子内. 通信子概念尤其方便了并行库的开发. MPI-1只支持组内通信 (intra-communication) MPI-2支持组间通信 (inter-communication) 2 MPI中的消息信封 率还辑戈蹄吕仅业茅淮嘶兔归