Форум программистов, компьютерный форум, киберфорум
Наши страницы

C++

Войти
Регистрация
Восстановить пароль
 
Mellotron
0 / 0 / 0
Регистрация: 10.09.2012
Сообщений: 31
#1

C++ MPI выделение потоков - C++

14.01.2014, 21:45. Просмотров 853. Ответов 1
Метки нет (Все метки)

Здравствуйте, я к вам с немного нубским вопросом. Написал вот такой код:

Кликните здесь для просмотра всего текста
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
#include <iostream>
#include <vector>
#include <ctime>
#include <mpi.h>
using namespace std;
 
void FillMatrix(const int n, const int m, int** mtrx, int* ivc);
void ReceiveMatrix(const int n, const int m, int** mtrx, int* rvc);
void PrintRcv(int vec[], int n, int &rng, int operation);
void sortVc(int* vector, int count);
 
void main(int argc, char **argv)
{
    setlocale(LC_ALL,"Russian");
    srand(time(NULL));
    int size, range;
    const int n=7;
    const int m=6;
    const int msize = m*n;
 
    int scount;
 
    int** matrix = new int* [n];
    int ivc[msize];
    int* rcv = new int[msize];
    for(unsigned i = 0; i < n; i++)
        matrix[i] = new int [m];
    MPI_Init(&argc, &argv);
    MPI_Comm_size (MPI_COMM_WORLD, &size);
    MPI_Comm_rank (MPI_COMM_WORLD, &range);
    scount = msize/size;
    int* tvc = new int[scount];
    if(range==0)
    {
        FillMatrix(n, m, matrix, ivc);      
    }
    
    if((msize%size)!=0)
    {
        scount++;
    }
    MPI_Scatter(&ivc[0],scount, MPI_INT, &tvc[0], scount, MPI_INT, 0, MPI_COMM_WORLD);
    PrintRcv(tvc, scount, range, 1);
    MPI_Barrier(MPI_COMM_WORLD);
    sortVc(tvc, scount);
    PrintRcv(tvc, scount, range, 2);
    //GATHERING HERE!
    MPI_Barrier(MPI_COMM_WORLD);
    MPI_Gather(&tvc[0], scount, MPI_INT, &rcv[0], scount, MPI_INT, 0, MPI_COMM_WORLD);
    
    MPI_Finalize();
    if(range == 0)
    {
        ReceiveMatrix(n, m, matrix, rcv);
    }
}
 
void FillMatrix(const int n, const int m, int** mtrx, int* ivc)
{
    int rc=0;
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < m; j++)
        {
            mtrx[i][j] = rand() % 64;
            cout<<"["<<i<<"]["<<j<<"]="<<mtrx[i][j]<<"\t";
            ivc[rc]=mtrx[i][j];
            rc++;
        }
        cout<<endl;
    }
}
 
void PrintRcv(int vec[], int n, int &rng, int operation)
{
    if(operation==1)
        cout<<"   Scattered Thread # "<<rng<<": ";
    if(operation==2)
        cout<<"   Sorting Thread # "<<rng<<": ";
    for(int i=0;i<n;i++)
    {
        cout.width(3);
        cout<<vec[i]<<" ";
    }
    cout<<endl;
}
 
void sortVc(int* vector, int count)
{
    for(int i = 0; i < count; i++)
    {
        for(int j = 0; j < count; j++)
            if(vector[i] < vector[j])
            {
                int tmp = vector[i];
                vector[i]=vector[j];
                vector[j]=tmp;
            }
    }
}
 
void ReceiveMatrix(const int n, const int m, int** mtrx, int* rvc)
{
    int rc = 0;
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < m; j++)
        {
            mtrx[i][j] = rvc[rc];
            rc++;
            cout<<"["<<i<<"]["<<j<<"]="<<mtrx[i][j]<<"\t";
        }
        cout<<endl;
    }
}


Код сортирует элементы в строках матрицы, но работает немного некорректно. Требуется, чтобы при любом количестве потоков сортировка элементов строки всегда работала правильно, а в моем коде она работает только если выделить количество потоков, кратное (лучше равное) количеству строк в матрице, иначе оно распределяет все элементы по количеству текущих потоков, и сортировка получается в корне неправильной. Пытался переделать так, чтобы потоки заново выделялись для обработки оставшихся строк - не вышло. Собственно, прошу помощи. Заранее спасибо.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
14.01.2014, 21:45
Здравствуйте! Я подобрал для вас темы с ответами на вопрос C++ MPI выделение потоков (C++):

Компиляция MPI - C++
Доброго времени суток всем, даже не знаю где создать тему, но так как программа написана на с++ решил здесь. В общем вопрос следующий: как...

MPI Распределение задач - C++
При инициализации процесса присваиваем N некое количество прохода цикла. При выполнении процессом всех проходов N нужно уменьшить N других...

MPI C++. Построение топологии сети - C++
Всем доброго времени суток. Задача следующая: каждый узел в сети знает только своих соседей(локальную топологию). Необходимо...

Передача в MPI vector<char*> - C++
Пишу программу на MPI. И возникли сложности. Имеется массив слов и его нужно передать между процессами. Сначала у меня был vector&lt;string&gt;,...

MPI использует только одно ядро - C++
MPI_Comm_size(MPI_COMM_WORLD, &amp;num); всегда возвращает 1. Ядра 100% два. В диспетчере задач проверил - Visual Studio использует оба ядра....

Вычисление суммы чисел, от параллельных процессов MPI - C++
Задача такая, запускаются параллельно 5 процессов, каждый из них генерит рандомное число и отправляет соседнему процессу, нужно вывести...

1
gng
645 / 491 / 133
Регистрация: 08.09.2013
Сообщений: 1,326
17.01.2014, 00:05 #2
Mellotron, если я вас правильно понял, то, например, так:
Групповыми операциями, как у вас, похоже, будет сложнее, так как придется динамически менять группу и создавать коммутатор.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
int main(int argc, char **argv)
{
  int size, range;
  const int n=7;
  const int m=6;
  const int msize = m*n;
 
  int** matrix = new int* [n];
  for(unsigned i = 0; i < n; i++) matrix[i] = new int [m];
  int ivc[msize];
  int* rcv = new int[msize];
  int* tvc = new int[n];
  MPI_Init(&argc, &argv);
  MPI_Comm_size (MPI_COMM_WORLD, &size);
  MPI_Comm_rank (MPI_COMM_WORLD, &range);
  if (!range) FillMatrix(n, m, matrix, ivc);
 
  MPI_Status st= {0};
  MPI_Barrier(MPI_COMM_WORLD);
  int count= 0;
  if(range==0) for (int k= 0; k < n + size - 1; k++) {
      MPI_Recv (tvc, m, MPI_INT, MPI_ANY_SOURCE,  MPI_ANY_TAG, MPI_COMM_WORLD, &st);
      MPI_Get_count (&st, MPI_INT,&count);
      if (count == m) memcpy (rcv + m * st.MPI_TAG, tvc, m * sizeof (int));
      count=  (k < n) ? m : 0;
      MPI_Send (ivc + k * m, count, MPI_INT, st.MPI_SOURCE, k, MPI_COMM_WORLD);
  }
  else for (;;) {
      MPI_Send (tvc, count, MPI_INT, 0, st.MPI_TAG, MPI_COMM_WORLD);
      MPI_Recv (tvc, m, MPI_INT, 0, MPI_ANY_TAG, MPI_COMM_WORLD, &st);
      MPI_Get_count (&st, MPI_INT,&count);
      if (!count) break;
      PrintRcv(tvc, m, range, 1);
      sortVc(tvc, m);
      PrintRcv(tvc, m, range, 2);
  }
  MPI_Finalize();
 
  if(range == 0) ReceiveMatrix(n, m, matrix, rcv);
  return 0;
}
1
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
17.01.2014, 00:05
Привет! Вот еще темы с ответами:

Вопрос про mpi (рассылка элементов массива) - C++
Привет всем, вопрос у меня следующий: можно ли с помощью mpi организовать одновременный доступ ко всему массиву для всех...

Error LNK2019 при компиляции программы MPI - C++
Здравствуйте! Это уже становится чересчур.. все решения перепробованы. Вобщем, на VS 2012 устанавливаю MS SDK MPI. В свойствах проекта...

Параллельная программа для метода холецкого с помощью openMp и mpi - C++
Товарищи,помогите пожалуйста с параллельным программированием: надо написать параллельную программу для метода холецкого с помощью openMp...

Создание потоков - C++
Всем доброго времени суток. Как всегда, последняя надежда на этот сайт. Сразу к делу: Написать программу, которая принимает в качестве...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru