MPI归并

/* Reduces values on all processes to a single value 将所有进程的值归并为一个值 */ #include "stdio.h" #include "mpi.h" #include "stdlib.h" int main(int argc, char * argv[]) { int myid, nprocs; int i; int send[3]; int recv[3]; MPI_Init(&argc, &argv); /*init MPI*/ MPI_Comm_rank(MPI_COMM_WORLD, &myid); /* get current process id */ MPI_Comm_size(MPI_COMM_WORLD, &nprocs); /* get number of processes */ printf("--------------------------------------\n"); printf("myrank=%d\n", myid); /*给每个进程的send缓存区填充数据*/ for (i = 0; i < 3; i++) { send[i] = myid + i; printf("send[i]=%d ", send[i]); } printf("\n"); /*------------------------------MPI_Reduce详解-----------------------------*/ //功能： // 将每个进程的send缓存区数据进行MPI_Op[运算]归并，结果储存在根进程的recv数 // 据缓存中 // 示例：如send缓存为int send[N],recv缓存为int recv[N],对send缓存 // 数据进行求和归并。 // 则根进程recv[i]=proc0(send[i])+proc1(send[i])+...+procN(send[i]) // N位该通信域进程的总数,proc(0-N)为进程0到进程N //函数参数： // MPI_Reduce(发送数据缓存地址，接受数据缓存地址，数据个数，数据类型，数据操作， // 根进程，通信域) MPI_Reduce(send, recv, 3, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD); /*结果输出*/ if (myid == 0) { for(i = 0;i<3;i++) printf("recv buf %d\n", recv[i]); } MPI_Finalize(); return 0; }