我是MPI的新手,我只是在编写一个基本的发送和接收模块,在这个模块中我将发送12个月到n个处理器并接收每个月并打印其值。所以我能够正确地发送值并且能够接收所有这些值但是我的程序被卡住了,即最后没有打印“程序完成后”。你能帮忙吗?
#include <stdio.h>
#include <string.h>
#include "mpi.h"
#include<math.h>
int main(int argc, char* argv[]){
int my_rank; /* rank of process */
int p; /* number of processes */
int tag=0; /* tag for messages */
MPI_Status status ; /* return status for receive */
int i;
int pro;
/* start up MPI */
MPI_Init(&argc, &argv);
// find out process rank
MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
//find out number of processes
MPI_Comm_size(MPI_COMM_WORLD, &p);
if (my_rank==0)
{
for(i=1;i<=12;i++)
{
pro = (i-1)%p;
MPI_Send(&i, 1, MPI_INT,pro, tag, MPI_COMM_WORLD);
printf("Value of Processor is %d Month %d\n",pro,i);
}
}
//else{
for(int n=0;n<=p;n++)
{
MPI_Recv(&i, 1, MPI_INT, 0, tag, MPI_COMM_WORLD, &status);
printf("My Month is %d and rank is %d\n",i,my_rank);
}
//}
MPI_Barrier(MPI_COMM_WORLD);
if(my_rank==0)
{
printf("After program is complete\n");
}
/* shut down MPI */
MPI_Finalize();
return 0;
}
Below is the output:
Value of Processor is 0 Month 1
Value of Processor is 1 Month 2
Value of Processor is 2 Month 3
Value of Processor is 3 Month 4
Value of Processor is 4 Month 5
Value of Processor is 0 Month 6
Value of Processor is 1 Month 7
Value of Processor is 2 Month 8
Value of Processor is 3 Month 9
Value of Processor is 4 Month 10
Value of Processor is 0 Month 11
My Month is 2 and rank is 1
My Month is 7 and rank is 1
My Month is 3 and rank is 2
My Month is 8 and rank is 2
Value of Processor is 1 Month 12
My Month is 1 and rank is 0
My Month is 6 and rank is 0
My Month is 11 and rank is 0
My Month is 12 and rank is 1
My Month is 4 and rank is 3
My Month is 9 and rank is 3
My Month is 5 and rank is 4
My Month is 10 and rank is 4
答案 0 :(得分:1)
首先:你违反了MPI的一个基本规则,你必须将一个发送与一个接收匹配。
在您的示例运行中,您运行5个处理器(排名),您可以看到排名0发送3条消息到排名0和1和2消息到剩余排名。但是,每个职位帖子13都会收到。因此,他们自然会卡在等待从未发送过的消息。请记住,MPI_Recv
周围的代码由所有等级执行。所以总共会收到5 * 13的收据。
如果轮到接收,你可以通过在循环内部进行过滤来解决这个问题。但这取决于你是否真的事先知道了等级0将发送多少消息 - 你可能需要更复杂的机制。
第二名:您排名0会向自己发送阻止消息(不先发布非阻塞接收)。这可能已经导致僵局。请记住,MPI_Send
永远不会保证在匹配的收到发布之前返回,即使有时可能会在实践中。
第三:该循环for(int n=0;n<=p;n++)
运行13次。你肯定不希望这样,即使你运行12次也不正确。
最后:对于特定示例,首选解决方案是将月份保存在数组中,并使用MPI_Scatterv
将其分布在所有进程中。