Реализация матричных операций с использованием коллективных взаимодействий процессов

[Novichok2016]

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

Ребят, подскажите, пожалуйста, где ошибка? Не обрабатывает первые два элемента первой строки матрицы. Запуск программы в ОС Linux. К сожалению, не могу прикрепить скрин, т.к. нет под рукой Linux

#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;
}