При распараллеливании программы и запуске её в Linux не обрабатываются первые два элемента первой строки матрицы

[Novichok2016]

ura_111, спасибо Вам! Я поняла. Не слишком ли я много занимаю Вашего времени? Уже не удобно к Вам обращаться, но если можете посмотрите, пожалуйста, еще программу...
Тема: "Реализация матричных операций с использованием коллективных взаимодействий процессов".

Задание: Дана действительная матрица A(n, n). Преобразовать эту матрицу, поэлементно вычитая значение последней строки из значений всех её строк, кроме последней. Запуск программы в ОС Linux. К сожалению, не могу прикрепить скрин, т.к. нет под рукой Linux...

Проблема в том, что при распараллеливании программы и запуске её в Linux не обрабатываются первые два элемента первой строки матрицы. В последовательном варианте и запуске в Windows (без MPI) - программа работает!

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<cmath>
#include "mpi.h"
#include<iomanip>

using namespace std;
const int N = 3;

int main(int argc, char *argv[])
{
      setlocale(0, "");
      int i, j, rank=0, size = 1, chunk, mas2[N][N], sum = 0;
      double mas[N][N], mas1[N][N];

      MPI_Status status;
      MPI_Init(&argc, &argv);
      MPI_Comm_rank(MPI_COMM_WORLD, &rank);
      MPI_Comm_size(MPI_COMM_WORLD, &size);

      if (rank == 0)
     {
        for (i = 0; i < N; i++)
       {
         for (j = 0; j < N; j++)
         {
           mas[i][j] = rand() % 10;
          }
        }
        chunk = N / size;
       }

       MPI_Bcast(&chunk, 1, MPI_INT, 0, MPI_COMM_WORLD);

       cout << endl;
       cout << "Первоначальная матрица: " << endl;
       cout << endl;
       for (i = 0; i < chunk; i++)
      {
        for (j = 0; j < chunk; j++)
       {
         mas1[i][j] = rand() % 10;
         cout << setw(3) << mas1[i][j] << " ";
       }
       cout << endl << endl;
      }

      MPI_Scatter(mas, chunk, MPI_DOUBLE, mas1, chunk, MPI_DOUBLE, 0,     MPI_COMM_WORLD);

      for (int i = 0; i < chunk - 1; ++i)
     {
       for (int j = 0; j < chunk; ++j)
      {
        mas2[i][j] = mas1[i][j] - mas1[N - 1][j];
        mas2[N - 1][j] = mas1[N - 1][j];
      }
     }

      MPI_Gather(mas1, chunk, MPI_DOUBLE, mas, chunk, MPI_DOUBLE, 0, MPI_COMM_WORLD);

      if (rank == 0)
     {
       cout << endl;
       cout << "Преобразованная матрица: " << endl << endl;
       for (i = 0; i < N; i++)
      {
        for (j = 0; j < N; j++)
        {
         cout << setw(3) << mas2[i][j] << " ";
        }
         cout << endl << endl;
        }
     }

    MPI_Finalize();
    system("pause");
    return 0;
}

[ura_111]

Пока я не могу. Выложи программу в новой ветке (на всеобщее обозрение). Может кто-нибудь раньше просмотрет.

У меня нет люнекса.

[Novichok2016]

Сообщение от ura_111

Пока я не могу. Выложи программу в новой ветке (на всеобщее обозрение). Может кто-нибудь раньше просмотрет.

У меня нет люнекса.

Я выкладывала уже, просмотров много, а ответов нет. А в Windows ещё... Необходимо возвести матрицу во вторую степень. Последовательный код я написала, он работает....

#include "stdafx.h"
#include <stdio.h>
#include <iostream>
#include <iomanip>
#include <ctime>

using namespace std;

int main()
{
	setlocale(0, "");
	int const n = 3;
	int a[n][n], c[n][n];
	cout << endl;
	cout << "Матрица А: " << endl << endl;
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			a[i][j] = rand() % 10 + 1;
			cout << setw(3) << a[i][j] << "      ";
		}
		cout << endl;
	}
	
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			c[i][j] = 0;
			for (int k = 0; k < n; ++k)
			{
				c[i][j] = c[i][j] + a[k][j] * a[i][k];
			}
		}
	}
	cout << endl;
	
	cout << "Матрица А*А: " << endl << endl;
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			cout << setw(3) << c[i][j] << "      ";
		}
		cout << endl;
	}
	cout << endl;
	
	system("pause");
	return 0;
}

Но нужно решить задачу параллельно с использованием метода
сдваивания. По описанию метода обработку двумерного массива нужно свести к одномерному виду... Например, если бы по заданию требовалось организовать параллельное выполнение умножения матрицы на вектор, по последовательный код ...

	
	
for (i = 0; i < m; i++)
{
   c[i] = 0;
   for (j = 0; j < n; j++)
   {
     c[i] += A[i][j]*b[j]
   }
}

... заменился бы на функцию потока...

void SerialResultCalculation(double* pMatrix, double* pVector, double* pResult, int Size, int N)
{
   int i, j; // переменные для организации цикла
   for (i=0; i<Size; i++) 
   {
     for (j=0; j<N; j++)
     pResult[i] += pMatrix[i*N+j]*pVector[j];
   }
}

Вопрос: не могу написать функцию потока для условия своей задачи (возведения матрицы в "квадрат")

[ura_111]

Мне нужно твоя помощь. Ты мне расскажи все что знаешь и что не получается. Только не здесь на ютубе (у меня канал на нём). Я создал новое видео, ссылка вот:

https://www.youtube.com/watch?v=G-C2dn-Wg1k&edit=vd

Под ним, в качестве комментариев, и будем общаться. Создай себе профиль на ютубе (под именем Novichok2016).
Нужно синхронизовать время выхода на связь. Вместе будем решать твою задачу.

[Alex11223]

В чем прикол общаться в комментах ютуба? Они ж неудобные.

Не проще уж тогда например скайп или вк взять? (из самого популярного)

[ura_111]

Никакого прикола нет, просто первое что пришло в голову.
Мне все равно где общаться.

Хотя для скайпа у меня нет оборудования (веб, гарнитуры). А вк - фамилию не хочется разглашать.

[Alex11223]

В скайпе текстовый чат есть.
Который не грузится 4 секунды при каждом обновлении страницы и не скрывает сообщения

[Smitt&Wesson]

#include <thread>

Старенький он уже, но гад, работает.

[Croessmah]

Кто старенький?
thread? Всего-то пять лет )))

[Novichok2016]

Сообщение от ura_111

Мне нужно твоя помощь. Ты мне расскажи все что знаешь и что не получается. Только не здесь на ютубе (у меня канал на нём). Я создал новое видео, ссылка вот:

https://www.youtube.com/watch?v=G-C2dn-Wg1k&edit=vd

Под ним, в качестве комментариев, и будем общаться. Создай себе профиль на ютубе.
Нужно синхронизовать время выхода на связь. Вместе будем решать твою задачу.

Вам какое время удобно? Я сейчас могу