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

MPI. полный чайник. разбор ошибок - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Потоки не выполняются параллельно http://www.cyberforum.ru/cpp-beginners/thread1763181.html
Почему потоки выполняются один за другим? using namespace std; void proc1(const int num) { int max_i = 50000; int i = 0; while(i<max_i)
C++ В чем ошибка в коде? #include <cstdlib> #include <iostream> using namespace std; int main(int argc, char *argv) { char rs='y'; if (rs == 'y') { http://www.cyberforum.ru/cpp-beginners/thread1763163.html
C++ Ширина бинарного дерева (нерекурсивно)
нерекурсивный метод определения ширины дерева с помощью стека. pop- извлечение из стека вершины дерева, push- занесение вершины, st1, st2 - стеки. ширина почему то выходит больше , т.е где то ошибка Struct node { uzel * info; node * next;};// стек, куда помещаем корень дерева Struct uzel {int info; uzel* next}; // функция опред. ширины
Функции max и min C++
Доброго времени суток) Как записать функции max и min? Сколько лазил, но понять не могу, как сформировать их по этому примеру. Подскажите пожалуйста) Заранее спасибо)
C++ Вычислить значение выражения с помощью процедур http://www.cyberforum.ru/cpp-beginners/thread1763134.html
Доброго времени суток) Что это такое и есть ли оно? Мне нужно создать прогу по формуле, но сколько не ищу, понять что такое процедуры в С++ никак не могу. Можете как-то подсобить или привести пример? Заранее спасибо!
C++ Создать класс Complex и перегрузить операторы Создайте класс с именем Complex для выполнения арифметических операций с комплексными числами. Напишите программу для тестирования вашего класса. Комплексные числа имеют вид: realPart + imaginaryPart*i, где i равно квадратному корню из -1. Используйте переменные с плавающей точкой для представления закрытых данных класса. Определите конструктор, который дает возможность инициализировать... подробнее

Показать сообщение отдельно
Маринчик
 Аватар для Маринчик
0 / 0 / 0
Регистрация: 12.10.2012
Сообщений: 123
16.06.2016, 15:25  [ТС]     MPI. полный чайник. разбор ошибок
Цитата Сообщение от AlexVRud Посмотреть сообщение
Теперь в каждом из потоков можно вычислять максимальную сумму элементов в чётных строках матриц
благодаря вашей помощи, получается вот так
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
// lab3!.cpp: определяет точку входа для консольного приложения.
 
#include "stdafx.h"
#include "mpi.h"
#include <stdlib.h>
#include <iostream>
using namespace std;
 
int _tmain(int argc, char* argv[])
{   int M,N;
int size, rank,i,j;
int lengthM=0, lengthN=0;/*, locallength=0;*/
double *A, *localA;;
/*int *displs, *counts;*/
int localM=0, localN=0;
 
double sum, localsum;
double maxsum;
 
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &size);//Определение количества процессов в параллельной программе 
MPI_Comm_rank(MPI_COMM_WORLD, &rank);//определение ранга процесса
 
 
if (!rank)
{    //ввод размеров матрицы
    cout<<"N = "; cin>>N;
    cout<<"vvedite chetnoe M = "; cin>>M;
    if (M%2==0)
    {
        A=new double[M*N];//динамическая матрица
        for ( i = 0; i < M; i++)  
         for ( j = 0; j < N; j++) 
          {A[i*N+j]=rand() % 100;///задаем
          }
 
        for ( i = 0; i < M; i++) // вывод 
        {for ( j = 0; j < N; j++) 
         printf("%7.2f",A[i*N+j]);
         printf("\n");
        }
    }
    else {printf("Ne sootvetstvie zadaniyu.Vvedte dannie zanovo"); return 0;}
}
//операция рассылки
 
MPI_Bcast(&M, 1, MPI_INT, 0, MPI_COMM_WORLD);
MPI_Bcast(&N, 1, MPI_INT, 0, MPI_COMM_WORLD);
 
localM=M/size;//+ (rank<(M%size));
localN=N;//Делим на полосы, а не квадратиками
localA = new double[localM*localN]; 
 
MPI_Scatter(A,localM*localN,MPI_DOUBLE,localA,localM*localN,MPI_DOUBLE,0, MPI_COMM_WORLD);
 
for (i=0; i<localM; i++)
 {for (j=0;j<localN; j++)
  cout <<localA[i*localN+j]<<" ";
  cout <<"\n "<<endl; 
 }
 
 for (i=0; i<localM; i+=2)//показываем на четность строк.. 0, 2, 4 ...
   {sum=0;
   { for (j=0;j<localN; j++)
     {// if (i%2==0)
      sum+=localA[i*localN+j];// сумма элементов четных строк 
   } }
 }   maxsum=sum;
   
 
MPI_Reduce (&maxsum,&localsum, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);//редукция
        if (!rank)
            printf("\n Maksimalnaya summa elementov chetnih  strok: %7.2f \n",localsum);
 
delete[]A;
delete[]localA;
/*delete[]displs;
delete[]counts;*/
MPI_Finalize();
 
system("pause");
return 0;
 
}
только вот не пойму,или мне уже кажется или и правда почему-то у меня по четным строкам не всегда максимальное из сумм вытягивает, если вручную просматривать..

Добавлено через 14 часов 51 минуту
А для Scatterv -как должно быть?? что не правильно делаю??

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
// lab3!.cpp: определяет точку входа для консольного приложения.
 
#include "stdafx.h"
#include "mpi.h"
#include <stdlib.h>
#include <iostream>
using namespace std;
 
int _tmain(int argc, char* argv[])
{   int M,N;
int size, rank,i,j,x,y;
int lengthM=0, lengthN=0;/*, locallength=0;*/
double *A, *localA;;
int *displs, *counts;
int localM=0, localN=0;
 
double sum, localsum;
double maxsum;
 
MPI_Init(&argc, &argv);
MPI_Comm_size(MPI_COMM_WORLD, &size);//Определение количества процессов в параллельной программе 
MPI_Comm_rank(MPI_COMM_WORLD, &rank);//определение ранга процесса
 
 
if (!rank)
{    //ввод размеров матрицы
    cout<<"N = "; cin>>N;
    cout<<"vvedite chetnoe M = "; cin>>M;
    if (M%2==0)
    {
        A=new double[M*N];//динамическая матрица
        for ( i = 0; i < M; i++)  
         for ( j = 0; j < N; j++) 
          {A[i*N+j]=rand() % 100;///задаем
          }
 
        for ( i = 0; i < M; i++) // вывод 
        {for ( j = 0; j < N; j++) 
         printf("%7.2f",A[i*N+j]);
         printf("\n");
        }
    }
    else {printf("Ne sootvetstvie zadaniyu.Vvedte dannie zanovo"); return 0;}
}
//операция рассылки
 
MPI_Bcast(&M, 1, MPI_INT, 0, MPI_COMM_WORLD);
MPI_Bcast(&N, 1, MPI_INT, 0, MPI_COMM_WORLD);
 
localM=M/size;//+ (rank<(M%size));
localN=N;//Делим на полосы, а не квадратиками
localA = new double[localM*localN]; 
counts= new int[size];
displs = new int [size];
//MPI_Scatter(A,localM*localN,MPI_DOUBLE,localA,localM*localN,MPI_DOUBLE,0, MPI_COMM_WORLD);
 
 
 
for (i=0; i<localM; i++)
 {for (j=0;j<localN; j++)
  cout <<localA[i*localN+j]<<" ";
  cout <<"\n "<<endl; 
 }
for (i=0; i<size; i++)
 for (j=0;j<size; j++)
   counts[i*size+j]=M/size+(i<(1%M/size));
{
   for (x=0; x<i; x++)
     for (y=0;y<j; y++)
      displs[x*j+y]=i/size+(x<(1%i/size));
}
MPI_Scatterv(A, counts, displs, MPI_DOUBLE, localA, localM*localN, MPI_DOUBLE, 0, MPI_COMM_WORLD);
 
 
 
for (i=0; i<localM; i+=2)//показываем на четность строк.. 0, 2, 4 ...
   {sum=0;
   { for (j=0;j<localN; j++)
     {// if (i%2==0)
      sum+=localA[i*localN+j];// сумма элементов четных строк 
   } }
 }   maxsum=sum;
   
 
MPI_Reduce (&maxsum,&localsum, 1, MPI_DOUBLE, MPI_MAX, 0, MPI_COMM_WORLD);//редукция
        if (!rank)
            printf("\n Maksimalnaya summa elementov chetnih  strok: %7.2f \n",localsum);
 
delete[]A;
delete[]localA;
delete[]displs;
delete[]counts;
MPI_Finalize();
 
system("pause");
return 0;
 
}
 
Текущее время: 03:01. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru