0

在下面的代码中,在探测时,未获得正确的到达数。使用标准 MPI 函数也测试了相同的功能,并获得了正确答案。为什么 Boost 版本不能产生正确的结果?

增强版:

#include <iostream>
#include <boost/mpi.hpp>

using namespace boost;
using namespace boost::mpi;

int main()
{ 
    environment env;
    communicator world;

    if (world.rank() == 0)
    {
        int a[70];

        auto req = world.isend(0, 0, a);
        //req.wait(); // does not make any difference on the result.

        optional<status> stat;
        while (!stat)
        {
            stat = world.iprobe(0, 0);

            if (stat)
            {
                optional<int> count = (*stat).count<int>();

                if (count)
                {
                    std::cout << *count << std::endl; // output: 2, expected: 70.
                }
            }
        }
    }

    return 0;
}

标准版本:

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

int main()
{
    MPI_Init(NULL, NULL);

    int rank;
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    if (rank == 0)
    {
        int a[70];            

        MPI_Request req;
        MPI_Status stat;

        MPI_Isend(a, 70, MPI_INT, 0, 0, MPI_COMM_WORLD, &req);
        //MPI_Wait(&req, &stat);

        int flag;
        MPI_Iprobe(0, 0, MPI_COMM_WORLD, &flag, &stat);
        if (flag)
        {
            int count;
            MPI_Get_count(&stat, MPI_INT, &count);
            std::cout << count << std::endl; // output: 70.
        }
    }

    MPI_Finalize();

    return 0;
}

编辑:使用isend(dest, tag, values, n)而不是isend(dest, tag, values)给出正确答案n数组中元素的数量在哪里。

4

1 回答 1

1

你的升压版其实送的不是70 int,而是一单int [70]。对于 Boost,此类型不是 MPI 数据类型,因此正确(*stat).count<decltype(a)>();返回一个空的可选项。

现在文档有点误导:

类型 T 必须具有关联的数据类型,即is_mpi_datatype<T>必须派生mpl::true_。在类型T与传输的类型不匹配的情况下,此例程将返回一个空的optional<int>.

我似乎相反,在T与传输类型不匹配的情况下,您将得到一个虚假的结果或一个空的optional<int>. 如果它不是 mpi 数据类型,则会得到一个空的optional<int>.

你得到 2 的原因是 Boost.MPI 为每个非 MPI 数据类型消息发送两条消息。一个包含序列化缓冲区的大小和实际消息。您的探针会处理大小消息,其中包含 1 size_t,其大小与 2 相同int

不幸的是,Boost.MPI 充满了与消息实际传输方式不同的此类微妙问题和错误。

于 2017-02-18T12:04:36.437 回答