2

要使用 MPI 运行我的软件，我需要创建一个简单的方法来返回MPI::COMM_WORLD.

所以在我的课堂上，我们有：

#include <mpi.h>
class Parallel{      
    public:

        MPI::Comm getCommunicator(){
            return MPI::COMM_WORLD;
        }
    protected:

    int iproc;
};

int main(int argc, char *argv[]){

    Parallel* parallel;
    MPI::Init(argc, argv);
    int my_rank;
    my_rank = parallel->getCommunicator().Get_rank();
    MPI::Finalize();
    return 0;
}

我应该如何实现该getCommunicator()方法才能返回MPI::COMM_WORLD？当我尝试编译上述内容时，出现以下错误：

invalid abstract return type for member function 'MPI::Comm Parallel ::getCommunicator()
4

1 回答 1

2

Zulan 是对的，从 MPI 3 中删除了 C++ 绑定，因此真的不应该使用它们编写新代码。

它们被删除的原因是它们没有经过特别深思熟虑或维护，而且它们绝对不是很惯用的，所以如果你使用它们会遇到很多奇怪的情况。Boost::MPI要好得多，但不幸的是只涵盖 MPI 1。

如果出于维护现有代码的目的，您必须拥有它，那么问题就如这里所述- 至少在 OpenMPI 中，MPI::Comm 被定义为纯虚拟类，因此您不能返回该类型的对象因为没有办法创造一个；你只能返回一个子类型。（我假设它是通过这种方式完成的，因此您可以将内部和内部通信器作为子类型）。

处理这种情况的经典方法是返回对对象的引用，而不是对象本身，并让编译器处理向上转换：

#include <mpi.h>
#include <iostream>

class Parallel{      
    public:
        MPI::Comm &getCommunicator(){
            return MPI::COMM_WORLD;
        }
    protected:
        int iproc;
};

int main(int argc, char *argv[]){

    Parallel* parallel;
    MPI::Init(argc, argv);
    int my_rank, size;

    my_rank = parallel->getCommunicator().Get_rank();
    size = parallel->getCommunicator().Get_size();
    std::cout << my_rank << "/" << size << std::endl;
    MPI::Finalize();
    return 0;
}

给

$ mpic++ -o foo foo.cpp
$ mpiexec -np 4 ./foo
1/4
2/4
3/4
0/4
于 2016-07-15T15:35:17.463 回答