MPI 发送和接收死锁答案

【问题标题】：MPI send and receive deadlockMPI 发送和接收死锁
【发布时间】：2016-06-15 17:15:50
【问题描述】：

我对 MPI 很陌生，我只是在编写一个基本的发送和接收模块，在该模块中我向 n 个处理器发送 12 个月，每个月接收并打印它的值。所以我能够正确发送值并且也能够接收所有这些值，但是我的程序被卡住了，即它最后没有打印“程序完成后”。你能帮忙吗？

#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

【问题讨论】：

标签： c mpi

【解决方案1】：

首先：你违反了 MPI 的一个基本规则，你必须匹配一个发送和一个接收。

在您的示例运行中，您使用 5 个处理器（等级）运行，您可以看到等级 0 向等级 0 发送 3 条消息，向剩余等级发送 1 和 2 条消息。但是，每个级别的职位都有 13 名。因此，他们自然会陷入等待从未发送过的消息。请记住，围绕MPI_Recv 的循环中的代码由所有等级执行。所以总共会有 5 * 13 次接收。

如果轮到你接收，你可以通过在循环内过滤来解决这个问题。但这取决于您是否真的事先知道等级 0 将发送多少条消息 - 您可能需要更复杂的机制。

第二次： 你排名 0 向自己发送一条阻塞消息（没有先发布非阻塞接收）。这已经可能导致死锁。请记住，MPI_Send 永远不会保证在匹配的接收发布之前返回，即使在实践中有时可能会返回。

第三次： 那个循环 for(int n=0;n<=p;n++) 运行了 13 次。你肯定不希望这样，即使你运行 12 次它是不正确的。

最后： 对于具体示例，首选的解决方案是将月份保存在一个数组中，并使用MPI_Scatterv 将其分布在所有进程中。

【讨论】：