将向量分散到可变 MPI 中

Scatter a vector into a variable MPI

本文关键字：MPI 向量 分散      更新时间：2023-10-16

我正在尝试从秩 = 0 的进程向所有其他进程 （1,2,3,4） 发送一个来自向量n_pointsBuffer的值。问题是只有进程 0 获得其值，而其他进程则没有。谁能向我解释一下我是否可能要做什么，如果是这样，如何做？这是我第一次使用 MPI。

#include <mpi.h>
#include <iostream>
#include <cstdlib>
#include <cmath>
using namespace std;
int main(int argc, char* argv[]) {
MPI::Init(argc, argv);
int num_procs = MPI::COMM_WORLD.Get_size();
int rank = MPI::COMM_WORLD.Get_rank();
srand(getpid());
int n_points;
if (rank == 0) {
    int n_pointsBuffer[] = { 1000000, 1203100, 1231230, 1231000, 1312322 };
    MPI::COMM_WORLD.Scatter(n_pointsBuffer, 1, MPI::INT, &n_points, 1,
            MPI::INT, 0);
}
cout << "Rank = " << rank << ", n_points = " << n_points << "n";
double sum = 0;
for (int i = 0; i < n_points; i++) {
    double x = rand() / ((double) (RAND_MAX));
    double f = 1.0 / (1.0 + x * x);
    sum += f;
}
double avg_sum = 0;
MPI::COMM_WORLD.Reduce(&sum, &avg_sum, 1, MPI::DOUBLE, MPI::SUM, 0);
if (rank == 0) {
    double pi = 4.0 * (avg_sum / num_procs / ((double) (n_points)));
    cout << "Pi is approx " << pi << " with error " << pi - M_PI << ".n";
}
MPI::Finalize();
return 0;
}

这是因为每个人都必须打电话给Scatter()。不仅仅是根。引用另一个答案：

MPI_Scatter函数包含发送和接收逻辑。这 根进程（此处指定为 0）发送数据，并且所有 接收者收到;每个参与者都必须调用例行程序。 分散是 MPI 集合操作的一个例子，其中所有任务 在沟通者中必须调用相同的例程。广播、屏障、 缩减操作和收集操作是其他示例。

请注意，每个集体行动都是一样的。这里有一个使用广播的例子（你也可以在同一个仓库中找到一些其他程序，这些是我写的，用来测试我的由转储机器组成的 MPI 集群）。

相关文章：