Форум программистов, компьютерный форум CyberForum.ru

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
White Luna
32 / 26 / 2
Регистрация: 08.09.2010
Сообщений: 402
#1

mpi распаралеливанние процессов когда потоков меньше чем сторк у матрицы - C++

22.04.2013, 20:06. Просмотров 243. Ответов 0
Метки нет (Все метки)

Задача умножение матрицы на вектор, это рабочий код когда количество поток и строк совпадают или кол-во строк меньше, а вот когда количество срок больше чем потоков с этим уже проблемы

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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
#include<iostream>
#include<mpi.h>
#include<conio.h>
 
int main(int argc, char ** argv)
{
  int mynode, totalnodes;
  int sum,startval,endval,accum;
 
  MPI_Status status;
 
  MPI_Init(&argc,&argv);
  MPI_Comm_size(MPI_COMM_WORLD, &totalnodes);
  MPI_Comm_rank(MPI_COMM_WORLD, &mynode);
 
  int  a1, a2;
  if(mynode == 0)
  {
    
   std::cout << "Vvedite razmer matrich stroku: ";  std::cin >> a1;
   std::cout << "Vvedite razmer matrich stolbech: " ;   std::cin>>a2;
   int w=0;
    int * vector6 =  new int[a2];
   int * vector =  new int[a2];
   
   int ** matr = new int *[a1];
   for(int i=0; i<a1; i++)
   {
       matr[i] = new int [a2];
   }
   for(int i=0;i<a2;i=i+1)
    {
   std::cout << "Vvedite element vectora" << i <<": "; std::cin>>vector[i];
   }
    for(int i1=0;i1<a1;i1=i1+1)
    {
         for(int i2=0;i2<a2;i2=i2+1)
        {
        std::cout << "Vvedite element matrich" << i1<<","<< i2 <<": "; std::cin>>matr[i1][i2];
        }
    }
 
std::cout<<"Vector "<<std::endl;
    for(int i=0;i<a2;i=i+1)
    {
std::cout<<vector[i]<<" ";
    }
    std::cout<<" "<<std::endl;
    std::cout<<"Matrich "<<std::endl;
     for(int i=0;i<a1;i=i+1)
    {
         for(int j=0;j<a2;j=j+1)
        {
            std::cout<<matr[i][j]<<" ";
         }
         std::cout<<" "<<std::endl;
     }
 
     
          
         for(int y=1; y<totalnodes; y++)
         {
               std::cout<<"Matr str "<<std::endl;
         MPI_Send(&a2,1,MPI_INT,y,1,MPI_COMM_WORLD); 
         MPI_Send(vector,a2,MPI_INT,y,2,MPI_COMM_WORLD); // отсылка данных
         MPI_Send(matr[y],a2,MPI_INT,y,3,MPI_COMM_WORLD);
         int  w2;
         MPI_Recv(&w2,1,MPI_INT,y,4,MPI_COMM_WORLD, &status);
         vector6[y]=w2;
         }
          
      
 
            for (int i = 0; i<a2 ;i++ ) 
            {
                std::cout<<matr[0][i]<<" ";
                w+=matr[0][i] * vector[i];
    std::cout<<" "<<std::endl;
            }
     vector6[0]=w;
std::cout<<"Umnogen strok v glav="<<w<<std::endl;
    
std::cout<<"Summa=";
for(int i=0; i<a1; i++)
{
std::cout<<vector6[i]<<" "<<std::endl;
}
 
 
  delete [] vector;
  delete [] matr;
  }
 
 
  if(mynode!=0)
  {
    
      int b1;
       MPI_Recv(&b1,1,MPI_INT,0,1,MPI_COMM_WORLD, &status);
      int * vector1 =  new int[b1];
          MPI_Recv(vector1,b1,MPI_INT,0,2,MPI_COMM_WORLD, &status);
      
       int * stroka_matr = new int[b1];
      MPI_Recv(stroka_matr,b1,MPI_INT,0,3,MPI_COMM_WORLD, &status);
      
 
  std::cout<<"Vector "<<std::endl;
    for(int i=0;i<b1;i=i+1)
    {
std::cout<<vector1[i]<<" ";
    }
    std::cout<<" "<<std::endl;
 
 
 
 std::cout<<"Matr str "<<std::endl;
 for(int i=0;i<b1;i=i+1)
    {
    std::cout<<stroka_matr[i]<<" ";
    }
std::cout<<" "<<std::endl;
 
/////////////////////////cчет
int w1=0;
for(int i=0; i<b1; i++)
{
w1+=stroka_matr[i] * vector1[i];
}
 
std::cout<<"Umnogen strok="<<w1<<std::endl;
 
 MPI_Send(&w1,1,MPI_INT,0,4,MPI_COMM_WORLD); 
      
  }
 
 
 
  MPI_Finalize();
  getch();
  return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
22.04.2013, 20:06     mpi распаралеливанние процессов когда потоков меньше чем сторк у матрицы
Посмотрите здесь:

Понятие синхронизации потоков и процессов. Семафоры и мьютексы - C++
Хотелось бы обсудить и уточнить некоторые понятия: 1. Корректно ли следующее утверждение. Синхронизация - это обеспечение...

Одновременное выполнение действий без использования потоков/процессов - C++
Нужно ожидать ввода команды со стандартного потока и в то же время в начале каждого часа вызывать некую функцию ( func назовем). ...

Записать условие, которое является истинной, когда только одно из чисел А, В и С меньше 50 - C++
Записать условие, которое является истинной, когда только одно из чисел А, В и С меньше 50. Записал вот так: if (a &lt; 50 || b &lt;...

[MPI] Как отправить кусок матрицы nxm - C++
Пишу программу умножения матрицы на вектор. Но она падает :(. Методом добавления вывода этапов пришёл к выводу что падает на принятии...

MPI. Как распараллелить алгоритм нахождение нормы матрицы - C++
собственно задача описана в заголовке. разобрал все функции MPI, но как применить к этому алгоритму - хз. везде в примерах легкотня одна. ...

Произведение на 1 меньше, чем должно быть - C++
Помогите пожалуйста с решением проблемы: если ввести k_usd = 0.13, а uah = 1000, то ответ будет 129 вместо 130, и эта ошибка расчетов...

Вывести информацию о холодильнике с ценой меньше, чем Z - C++
#include &lt;iostream&gt; #include &lt;cstring&gt; using namespace std; struct Fridge { string name; string company; int kol; ...

Объясните тупому почему задержка потока по Sleep() больше чем ... с вызывает задержки и других потоков? - C++
Какое отношение они имеют к данному потоку?Ну спит он себе и спит, пусть он выполняется раз в 5 секунд, допустим, чаще не надо-другие-то...

Выдает ошибку когда задействую поток.Подскажите в чем ошибка? - C++
Задача: Задан двумерный массив A из N x N элементов, расположенный в локальной памяти каждой нити, то есть массив A разделен на m...

Перенести в массив те элементы из массивов, которые меньше чем (max_a+max_b)/2 - C++
Заданы массивы а1, а2, … , аn и b1, b2, … , bm Перенести в массив С те элементы из массивов, которые меньше чем (max_a+max_b)/2


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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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