为什么不处理0从自身接收数据?

时间:2017-02-25 20:06:32

标签: c mpi blocking nonblocking

我刚刚开始学习MPI,我得到问题的原因并不清楚。

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


int main(int argc, char *argv[])
{
    int myRank, size;
    double myLeftBound, myRightBound;
    MPI_Status status1;
    MPI_Status status2;
    MPI_Request request1;
    MPI_Request request2;   

    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &myRank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);   


    if (myRank == 0)
    {
        int counter;

        for (counter = size - 1; counter >= 0; counter--)
        {
            myLeftBound = 0;
            myRightBound = 1;
            MPI_Isend(&myLeftBound, 1, MPI_DOUBLE, counter, 2 * counter, MPI_COMM_WORLD, &request1);
            MPI_Isend(&myRightBound, 1, MPI_DOUBLE, counter, 2 * counter + 1, MPI_COMM_WORLD, &request2);   
        }
    }

    MPI_Recv(&myLeftBound, 1, MPI_DOUBLE, 0, myRank * 2, MPI_COMM_WORLD, &status1);
    MPI_Recv(&myRightBound, 1, MPI_DOUBLE, 0, myRank * 2 + 1, MPI_COMM_WORLD, &status2);
    printf ("I received my left boundary %f\n, I'm %d\n", myLeftBound, myRank);
    MPI_Finalize();
    return 0;
}

其中一个输出正在跟随(进程数为4)

  

我收到了我的左边界0.000000

     

,我&#39; 1

     

我收到了我的左边界0.000000

     

,我&#39; 2

     

我收到了我的左边界0.000000

     

,我&#39; 3

     

在第77行的文件src / mpid / ch3 / src / ch3u_buffer.c中断言失败:FALSE   memcpy参数内存范围重叠,dst_ = 0x7ffcc26963d0 src_ = 0x7ffcc26963d0 len_ = 8

     

内部ABORT - 流程0

流程0出了什么问题?

1 个答案:

答案 0 :(得分:0)

您无法接收发送的同一缓冲区。从技术上讲,在完成非阻塞发送操作之前,不得修改发送缓冲区的任何部分。您还应该完成发送操作,这意味着您应该将其存储在多个请求对象中。

无论如何,只需使用MPI_Bcast并在您对MPI更有经验时留下非阻止操作。出于某种原因,它们经常被初学者过度使用和错误使用。