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

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

Восстановить пароль Регистрация
 
White Luna
 Аватар для White Luna
32 / 26 / 2
Регистрация: 08.09.2010
Сообщений: 402
22.04.2013, 20:06     mpi распаралеливанние процессов когда потоков меньше чем сторк у матрицы #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
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++ Объясните тупому почему задержка потока по Sleep() больше чем ... с вызывает задержки и других потоков?
C++ Одновременное выполнение действий без использования потоков/процессов
[MPI] Как отправить кусок матрицы nxm C++
Понятие синхронизации потоков и процессов. Семафоры и мьютексы C++
Произведение на 1 меньше, чем должно быть C++
Перенести в массив те элементы из массивов, которые меньше чем (max_a+max_b)/2 C++
MPI. Как распараллелить алгоритм нахождение нормы матрицы C++
Записать условие, которое является истинной, когда только одно из чисел А, В и С меньше 50 C++

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

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

Текущее время: 19:26. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru