向量的 MPI 发送/接收

Question

大家好，我正在从事一个涉及使用 MPI 库的大学项目，不幸的是我无法分享整个代码，但我希望有人能够给我一些指示。 我需要将一个数组从 proc X 发送到 proc 0，但是正如我在 google 上读到的，发送一个 dinamycally 创建的数组的唯一方法是找出大小，将大小发送到 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

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

在这里，您将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);