Понимание MPI_Type_vector - PullRequest
Новая
       2

Понимание MPI_Type_vector

1 голос
/ 29 марта 2019

Я ожидал, что следующий код выдаст 1,0,2,0, но он выдаст 1,2,0,0.Почему он не дополняет полученные данные? 

// tmp.cpp
#include <iostream>
#include <mpi.h>

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

    int rk;
    MPI_Comm_rank( MPI_COMM_WORLD, &rk );

    if ( rk == 0 )
    {
        double b[2] = {1.0, 2.0};
        MPI_Send( b, 2, MPI_DOUBLE, 1, 0, MPI_COMM_WORLD );
    }
    if ( rk == 1 )
    {
        double b[4] = {0.0, 0.0, 0.0, 0.0};

        MPI_Datatype padded;
        MPI_Type_vector( 1, 1, 2, MPI_DOUBLE, &padded );
        MPI_Type_commit( &padded );

        MPI_Recv( b, 2, padded, 0, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE );

        std::cout << b[0] << ',' << b[1] << ','
                  << b[2] << ',' << b[3] << '\n';

        MPI_Type_free( &padded );
    }

    MPI_Finalize();
    return 0;
}

Однако, если я изменяю принимаемую часть на 

        MPI_Datatype padded;
        MPI_Type_vector( 2, 1, 2, MPI_DOUBLE, &padded );
        MPI_Type_commit( &padded );

        MPI_Recv( b, 1, padded, 0, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE );

, он выдает 1,0,2,0.

MPI: OpenMPI 3.1.2

GCC: 7.2.1

Сборка: mpicxx tmp.cpp

Запуск: mpiexec -n 2 ./a.out

1 Ответ

0 голосов
/ 29 марта 2019

Видимо, это правильное поведение. Я нашел объяснение здесь , в частности, слайды 12-18.

...