0

После использование двух и более подряд функций MPI_Bcast ломается параллельность (далее выполняет код только 0 поток). Пример кода:

#include <iostream>
#include <mpi.h> 
#include <vector>
using namespace std;

int main(int argc, char** argv)
{
    int rank, p;
    int n = 0;
    vector<int> matr;
    vector<int> vec;
    vector<int> res;
    MPI_Comm comm = MPI_COMM_WORLD;
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(comm, &rank);
    MPI_Comm_size(comm, &p);
    if (rank == 0)
    {
        cin >> n;
        matr = vector<int>(n*n,2);
        vec = vector<int>(n,4);
        res = vector<int>(n,4);
    }

    MPI_Barrier(comm);

    MPI_Bcast(&n, 1, MPI_INT, 0, comm);
    MPI_Bcast(&matr, (n*n), MPI_INT, 0, comm);
    MPI_Bcast(&vec, n, MPI_INT, 0, comm);
    MPI_Bcast(&res, n, MPI_INT, 0, comm);

    cout << n;
    MPI_Finalize();
}

cout << n выполняет только 0 поток. После чего программа завершается. Если же попробовать вывести элементы какого-нибудь вектора, программа и вовсе крашится с ошибкой (см. скришнот). В чем проблема?

#include <iostream>
#include <mpi.h> 
#include <vector>
using namespace std;

int main(int argc, char** argv)
{
    int rank, p;
    int n = 0;
    vector<int> matr;
    vector<int> vec;
    vector<int> res;
    MPI_Comm comm = MPI_COMM_WORLD;
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(comm, &rank);
    MPI_Comm_size(comm, &p);
    if (rank == 0)
    {
        cin >> n;
        matr = vector<int>(n*n,2);
        vec = vector<int>(n,4);
        res = vector<int>(n,4);
    }

    MPI_Barrier(comm);

    MPI_Bcast(&n, 1, MPI_INT, 0, comm);
    MPI_Bcast(&matr, (n*n), MPI_INT, 0, comm);
    MPI_Bcast(&vec, n, MPI_INT, 0, comm);
    MPI_Bcast(&res, n, MPI_INT, 0, comm);


    for (int i = 0; i < vec.size(); i++)
        cout << vec[i] << " ";

    MPI_Finalize();
}

введите сюда описание изображения

5
  • ЗдесьMPI_Bcast(&matr, (n*n), MPI_INT, 0, comm); вы передали указатель на объект вектор (суть указатель на структуру из нескольких указателей), а не на буфер данных, подходящего размера (типа int*). В дочерних процессах, место под данные вовсе не выделено. Дальше справитесь?
    – Chorkov
    30 мая 2021 в 20:19
  • А как правильно? При использовании matr.data() выбрасывает исключение null buffer. 30 мая 2021 в 20:33
  • Даже если так, почему тогда при cout << n результат выводит только 0 поток? 30 мая 2021 в 20:34
  • Потому что UB. Выучи темы указателей, классов, контейнеров в C++, и только после этого используй их. 30 мая 2021 в 20:49
  • 1
    Падение на matr.data() , потому что вы не выделили память ни в одном процессе кроме главного. (для пустого вектора, data возвращает нулевой указатель).
    – Chorkov
    30 мая 2021 в 21:00

0

Ваш ответ

By clicking “Отправить ответ”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Посмотрите другие вопросы с метками или задайте свой вопрос.