Question

在下面的代码中，在探测时，未获得正确的到达次数。使用标准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是数组中元素的数量。

Answer 1

您的Boost版本实际上并不发送70 int，而只发送一个int [70]。对于Boost，此类型不是MPI数据类型，因此正确的(*stat).count<decltype(a)>();返回空的可选项。

现在documentation有点误导：

类型T必须具有关联的数据类型，即is_mpi_datatype<T>必须派生mpl::true_。如果类型T与传输类型不匹配，则此例程将返回空optional<int>。

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

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

不幸的是，Boost.MPI充满了微妙的问题以及与实际传输消息的不同方式有关的错误。

boost :: iprobe不会返回正确的计数

1 个答案: