Вызов функций при использовании параллельного программирования

У меня просто довольно простой вопрос об использовании MPI в программе на C++. Действительно, возьмем этот пример:

#include <iostream>
#include <sstream>
#include <cblas.h>
#include <cmath>

using namespace std;

#include <mpi.h>

void multiply(double* x,double* y,int tai,double dot){
    for(int i=0; i<tai;i=i+1){
        dot=dot+x[i]*y[i];
    }
}

int main(int argc, char* argv[]) {
    const int n=32;
    int rank;
    int size;

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

    int tai=n/size;
    double* x=new double[tai];
    double* y=new double[tai];
    srand(time(0)+rank);

    for(int i=0;i<tai;i=i+1){
        x[i]=(double)rand()/RAND_MAX*10;
        y[i]=(double)rand()/RAND_MAX*10;
    }
    double dot=0;
    multiply(x,y,tai,dot);
    double ddot;
    MPI_Reduce(&dot, &ddot, 1, MPI_DOUBLE, MPI_SUM,0, MPI_COMM_WORLD);
    if (rank==0){
        cout<<"product:"<<ddot<<endl;
    }
    MPI_Finalize();
    return 0;
}

Затем я вызываю свою функцию «умножение» в своей программе для умножения двух векторов, но, к сожалению, она возвращает начальное значение «точки» (которое в основном равно 0).

Я просто хотел бы знать, есть ли что-нибудь особенное в вызове функции при параллельном программировании.

PS: Я знаю, что вся программа работает, так как она дает мне хорошие результаты, когда я напрямую умножаю свои 2 вектора в функции "main".