简单的MPI计划

时间:2011-05-10 04:54:05

标签: c mpi

我正在尝试做一些MPI,这是一个使用MPI_Send,MPI_Recv(是的,阻止)的简单程序。等级0向所有其他进程发送消息,其他进程接收消息。但是,我的MPI_Send永远不会返回。我在这里错过了什么吗?谢谢!

#include<stdio.h>
#include<string.h>
#include<unistd.h>
#include"mpi.h"

int main(int argc, char **argv)
{

int rank, size, i;
MPI_Status status;
char message[20];
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &size);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
strcpy(message, "Hello, world");
if (rank==0){
   for(i=1;i<size;i++){

         printf("I am sending %s\n", message);
         MPI_Send(message, 23, MPI_BYTE, 0, 7, MPI_COMM_WORLD);
         printf("Sending node=%d, message=%s\n", rank, message);
 }
}
else{
                     MPI_Recv(message, 23, MPI_BYTE, MPI_ANY_SOURCE, 7, MPI_COMM_WORLD,&status);
                     printf("Receiving node=%d, message=%s\n", rank, message);
 }
MPI_Finalize();
return 0;
}

2 个答案:

答案 0 :(得分:11)

您正在向自己发送消息:

MPI_Send(message, 23, MPI_BYTE, 0, 7, MPI_COMM_WORLD);

我相信这应该是

MPI_Send(message, 23, MPI_BYTE, i, 7, MPI_COMM_WORLD);

答案 1 :(得分:4)

正如@cnicutar指出的那样,进程0正在向自己发送消息。然而,在我看来,MPI_Bcast将是一个更好的选择。它肯定会简单得多,因为它既是“发送”又是“接收”方法。