出于某种原因必须两次发送MPI消息

时间:2012-10-13 20:53:21

标签: c++ c multiprocessing mpi

我一直致力于一个MPI项目,其中所有从站都将数据发送回主站。出于某种原因,如果我连续连续发送2次,主设备将仅接收数据。这很奇怪,我认为这会引起一些其他奇怪的问题。知道是什么会导致这个吗?我认为第一个发送是发送某种垃圾数据或其他东西。发送是完全相同的代码行。

编辑:以下代码......

if (numProcs > 0)
    MPI_Barrier( MPI_COMM_WORLD ) ; //only wait if there are other processes to wait for

if (rank != 0)
{
    MPI_Send(handArray, 10, MPI_DOUBLE, 0, TAG_HAND, MPI_COMM_WORLD);
    MPI_Send(handArray, 10, MPI_DOUBLE, 0, TAG_HAND, MPI_COMM_WORLD);
}
//8. After the main loop the master process receives and sums together the hand counts array
//   from each slave process.
else
{
    int activeProcs = numProcs - 1;
    getHandsFromSlaves( activeProcs, handArray );

然后主人继续打印一些数据......

这是getHands FromSlaves方法。请注意我也尝试过使用阻塞调用以及同样的问题。

void getHandsFromSlaves( int& activeCount, double handTotals[10] ){

static MPI_Request request;
static int msgBuff, recvFlag;
static double handBuff[10];
MPI_Status status;

while (activeCount > 0)
{
    if( request )
    {
        // Already listening for a message

        // Test to see if message has been received
        MPI_Test( &request, &recvFlag, &status );
        //cout << "TAG: " << status.MPI_TAG << " SOURCE: "<< status.MPI_SOURCE  << " ERROR: " << status.MPI_ERROR << endl;
        if( recvFlag )
        {
            // Message received
            if( status.MPI_TAG == TAG_HAND )
            {
                cout << "Hand Received!" << endl;

                for(int m = 0; m < 10; ++m)
                {
                    handTotals[m] += handBuff[m];
                }

                activeCount--;
            }
            else
            {
                //error report... what happened?
                cout << "TAG: " << status.MPI_TAG << " SOURCE: "<< status.MPI_SOURCE  << " ERROR: " << status.MPI_ERROR << endl;
            }

            // Reset the request handle
            request = 0;
        }
    }

    if( !request && activeCount > 0 )
        // Start listening again
        MPI_Irecv(&handBuff, 10, MPI_DOUBLE, MPI_ANY_SOURCE, MPI_ANY_TAG, MPI_COMM_WORLD, &request);
}
}

1 个答案:

答案 0 :(得分:1)

好吧,您可能正在尝试处理太多邮件,因为在输入request例程时您的getHandsFromSlaves()变量未定义。从输入开始,request几乎肯定不为零,即使您尚未发布MPI_Test,您也会立即尝试Irecv发送消息。

事实上,这里发布的代码摘录有很多非常奇怪的事情。为什么局部变量static?为什么要在MPI_Test()而不是MPI_Wait()上实施自己的busywait?如果你在接收之间没有做任何有用的事情,你为什么要使用非阻塞接收?事实上,如果你只是总结所有的数组,为什么你要做个别的点对点接收,而不是做MPI_Reduce()

以下更短的代码看起来像你上面要做的那样:

#include <stdio.h>
#include <mpi.h>


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

    int rank, numProcs;

    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &numProcs);
    double handArray[10];
    double handTotal[10];

    for (int i=0; i<10; i++)
        handArray[i] = rank + i;

    if (rank == 0)  // Since apparently rank 0 doesn't do anything
    {
        for (int i=0; i<10; i++)
            handArray[i] = 0;
    }

    MPI_Reduce(handArray, handTotal, 10, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);

    if (rank == 0) {
        printf("Hand Totals= \n");
        for (int i=0; i<10; i++)
            printf(" %lf ", handTotal[i]);
        printf("\n");
    }

    MPI_Finalize();
}