Параллельное программирование mpi

#include "stdafx.h"
#include <cstdlib>
#include <iostream>
#include <mpi.h>
#include <ctime>
#define ARRAY_SIZE 10
 
using namespace std;
 
int main(int argc, char* argv[])
{
    double time_start=0, time_end, TIME;
    int size=0, rank=0;
    int *vec1, *vec2,*vec3;
    const double size1 = ARRAY_SIZE;
    MPI_Status status;
 
    MPI_Init(&argc, &argv);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
 
    double mid = size1 / size;//кол-во элементов на кол-во процессов
    double size_for_procces = 0;
    float ost = modf(mid, &size_for_procces);
    int dop = 0;
 
 
//заполнение векторов
    if (rank)
    {
        vec1 = new int[(int)size_for_procces];
        vec2 = new int[(int)size_for_procces];
        for (int i = 0; i < size_for_procces; i++)
        {
            vec1[i] = rand() % 10 + 1;
            vec2[i] = rand() % 10 + 1;
        }
    }
    else
    {
        if (ost > 0)
            dop = round(ost*size);
 
        vec1 = new int[(int)size_for_procces + dop];
        vec2 = new int[(int)size_for_procces + dop];
        for (int i = 0; i < size_for_procces + dop; i++)
        {
            vec1[i] = rand() % 10 + 1;
            vec2[i] = rand() % 10 + 1;
            cout << vec1[i]<< ' '<<' ';
        }
        cout << endl;
        for (int i = 0; i < size_for_procces + dop; i++)
        {
            cout << vec2[i] << ' '<<' ';
        }
    }
    //измерение времени
    if (rank == 0)
        time_start = MPI_Wtime();
 
    //поиск суммы векторов
    vec3= new int[(int)size_for_procces + dop];
    if (rank == 0)
    {
        for (int i = 0; i < size_for_procces + dop; i++)
            vec3[i] = vec1[i] + vec2[i];
    }
    else
    {
        for (int i = 0; i < size_for_procces; i++)
        {
            vec3[i] = vec1[i] + vec2[i];
 
            MPI_Send(&vec3[i], 1, MPI_INT, 0, 2, MPI_COMM_WORLD);
        }
    }
    cout << endl;
    if (rank == 0)
    {
        int *n_vec = new int[size];
        n_vec[0] = vec3[1];
 
        for (int i = 1; i < size; i++)
            MPI_Recv(&n_vec[i], 1, MPI_INT, i, 2, MPI_COMM_WORLD, &status);
 
        for (int i = 1; i < size; i++)
        {
            vec3[i] = vec1[i] + vec2[i];
        //out << vec3[i] << ' ';
        }
        time_end = MPI_Wtime();
        TIME = (time_end - time_start);
        for (int i = 0; i < size_for_procces; i++)
        cout << vec3[i] << ' '<<' ';
        cout << "\n TIME: " << TIME * 1000 << " milliseconds" << endl;
        delete[] n_vec;
    }
    delete[] vec1;
    delete[] vec2;
    delete[] vec3;
    MPI_Finalize();
    return 0;
}