MPI:多个重叠的通信器

MPI:多个重叠的通信器,mpi,communicator,Mpi,Communicator,我想创建MPI通讯器,将秩为0的进程链接到其他每个进程。假设n是进程的总数。然后,秩为0的进程应该有n-1个通信器,而其他每个进程都有一个通信器。这可能吗?如果可能,为什么我不能使用下面的程序来实现这一点 使用mpic++编译器编译下面的代码会在我的计算机上终止,而不会出现警告和错误。但是,当我使用3个或更多进程(mpiexec-n3)运行生成的程序时,它永远不会终止 我可能误解了MPI中的沟通者的概念。也许有人能帮助我理解为什么下面的程序会被卡住,还有什么更好的方法来创建这些沟通者?谢谢 #i

我想创建MPI通讯器,将秩为0的进程链接到其他每个进程。假设n是进程的总数。然后,秩为0的进程应该有n-1个通信器,而其他每个进程都有一个通信器。这可能吗?如果可能,为什么我不能使用下面的程序来实现这一点

使用mpic++编译器编译下面的代码会在我的计算机上终止,而不会出现警告和错误。但是,当我使用3个或更多进程(mpiexec-n3)运行生成的程序时,它永远不会终止

我可能误解了MPI中的沟通者的概念。也许有人能帮助我理解为什么下面的程序会被卡住,还有什么更好的方法来创建这些沟通者?谢谢

#include <iostream>
#include <vector>
#include <thread>
#include <chrono>
#include "mpi.h"
void FinalizeMPI();
void InitMPI(int argc, char** argv);
int main(int argc, char** argv) {
    InitMPI(argc, argv);

    int rank,comm_size;
    MPI_Comm_rank(MPI_COMM_WORLD,&rank);
    MPI_Comm_size(MPI_COMM_WORLD,&comm_size);

    if(comm_size<2) {
        FinalizeMPI();
        return 0;
    }

    MPI_Group GroupAll;
    MPI_Comm_group(MPI_COMM_WORLD, &GroupAll);

    if(rank==0) {
        std::vector<MPI_Group> myGroups(comm_size-1);
        std::vector<MPI_Comm> myComms(comm_size-1);
        for(int k=1;k<comm_size;++k) {
           int ranks[2]{0, k};
           MPI_Group_incl(GroupAll, 2, ranks, &myGroups[k-1]);
           int err = MPI_Comm_create(MPI_COMM_WORLD, myGroups[k-1], &myComms[k-1]);
           std::cout << "Error: " << err << std::endl;
        }
    } else {
        MPI_Group myGroup;
        MPI_Comm myComm;
        int ranks[2]{0,rank};
        MPI_Group_incl(GroupAll, 2, ranks, &myGroup);
        int err = MPI_Comm_create(MPI_COMM_WORLD, myGroup, &myComm);
        std::cout << "Error: " << err << std::endl;
    }
    std::cout << "Communicators created: " << rank << std::endl;
    std::this_thread::sleep_for(std::chrono::milliseconds(1000));

    FinalizeMPI();
    return 0;
}

void FinalizeMPI() {
   int flag;
   MPI_Finalized(&flag);
   if(!flag)
      MPI_Finalize();
}

void InitMPI(int argc, char** argv) {
   int flag;
   MPI_Initialized(&flag);
   if(!flag) {
      int provided_Support;
      MPI_Init_thread(&argc, &argv, MPI_THREAD_MULTIPLE, &provided_Support);
      if(provided_Support!=MPI_THREAD_MULTIPLE) {
          exit(0);
      }
   }
}
#包括
#包括
#包括
#包括
#包括“mpi.h”
void-empi();
void InitMPI(int argc,char**argv);
int main(int argc,字符**argv){
InitMPI(argc,argv);
整数秩,通信大小;
MPI通信等级(MPI通信世界和等级);
MPI_通信大小(MPI_通信世界和通信大小);

如果(comm_size
MPI_comm_create
是初始通信器上的一个集合操作(
MPI_comm_WORLD
),则必须在所有进程上调用它

解决此问题的最简单方法是按您的方式使用
MPI\u Comm\u create\u group
。它与
MPI\u Comm\u create
类似,只是它是整个组的集合