矢量的 MPI 发送/接收

Question

大家好，我正在开展一个涉及使用 MPI 库的大学项目，不幸的是我无法分享整个代码，但我希望有人能够给我一些指导。 我需要将一个数组从 proc X 发送到 proc 0，但是正如我在 google 上读到的那样，发送动态创建的数组的唯一方法是找出大小，将大小发送到 Proc 0，然后才发送数组（我无法预先确定它的大小），所以这就是我在 Proc X 上所做的：

vector<int> tosend = aglib.ExtractEvaluationFunctions();
int tosend_size = tosend.size();
MPI_Send(&tosend_size, 1, MPI_INT, 0, 666, MPI_COMM_WORLD);
MPI_Send(&tosend, tosend_size, MPI_INT, 0, 777, MPI_COMM_WORLD);

这就是 Proc 0 上发生的情况（我无法共享相同的缓冲区，因为要发送的向量每次都是在 Proc1 中本地创建的）：

 vector<int> results_and_rank;
 int results_rank_size;
 MPI_Recv(&results_rank_size, 1, MPI_INT, MPI_ANY_SOURCE, 666, MPI_COMM_WORLD, &status);
 MPI_Recv(&results_and_rank, results_rank_size, MPI_INT, MPI_ANY_SOURCE, 777, MPI_COMM_WORLD, &status); 
 cout << "size of results_and_rank is :"<< results_and_rank.size()endl;
 cout<< "last element of array is :"<< results_and_rank.back()<<endl;

我认为通信工作正常，因为我能够读取接收到的向量的大小，该向量与我发送的向量相同，但是每当我尝试访问数组的元素时，我的代码就会崩溃

results_and_rank

，从而在最后一个时崩溃打印。

顺便说一下，我需要使用阻塞通信来实现我的项目。

我错过了什么吗？感谢您的时间和帮助。

Answer 1

您想要向向量发送/接收data，因此您需要向量内部的数据，而不是向量本身。

在这里，您将

std::vector

本身的地址传递给

MPI_Recv()

，这会损坏

std::vector

。

MPI_Recv(&results_and_rank, results_rank_size, MPI_INT, MPI_ANY_SOURCE, 777, MPI_COMM_WORLD, &status);

正确做法：

发件人：

MPI_Send(&tosend[0], tosend_size, MPI_INT, 0, 777, MPI_COMM_WORLD);

接收者：

results_and_rank.resize(results_rank_size);
MPI_Recv(&results_and_rank[0], results_rank_size, MPI_INT, MPI_ANY_SOURCE, 777, MPI_COMM_WORLD, &status);