Задача на распараллеливание с подпрограммой

#include <iostream>
#include <clocale>
#include <cstdlib>
using namespace std;
void Input(int* M, int n);  
void Print(int* M, int n);   
void RandomFill(int* M, int n);   
void Work(int* M1, int* M2, int* res, int n);   
 
 
int main()
{
    setlocale(LC_ALL, "Rus");
    int n;
    cout << "Размер векторов?";
    cin >> n;
    int* A = new int[n];
    int* B = new int[n];
    int* C = new int[n];
 
    //Input(A,n);
    RandomFill(A, n);
    cout << "Первый массив\n";
    Print(A, n);
 
    RandomFill(B, n);
    cout << "Второй массив\n";
    Print(A, n);
 
    Work(A, B, C, n);
    cout << "Сумма\n";
    Print(C, n);
    delete[] A;
    delete[] B;
    delete[] C;
    system("pause");
    return 0;
}
 
void Input(int* M, int n)
{
    for (int i = 0; i < n; i++)
    {
        cout << "M[" << i << "]=?";
        cin >> M[i];
    }
}
 
void RandomFill(int* M, int n)
{
    for (int i = 0; i < n; i++)
    {
        M[i] = rand() % 100;
    }
}
 
void Print(int* M, int n)
{
    for (int i = 0; i < n; i++)
    {
        cout << M[i] << " ";
    }
    cout << endl;
}
 
void Work(int* M1, int* M2, int* res, int n) 
{
    for (int i = 0; i < n; i++)
    {
        res[i] = M1[i] + M2[i];
    }
}

#define VECLEN 100
float a[VECLEN], b[VECLEN];
float s;
 
float dotprod ()
{
int i;
 
s = 0;
#pragma omp for reduction(+:s)  reduction используется только если нужна общая сумма
  for (i=0; i < VECLEN; i++)
    {
    s = s + (a[i]*b[i]);    здесь формулы согласно задаче
    }
return(s);  если ответ: массив, то передавайте через ссылку или указатель
}
 
void main()
{
double t1, t2, sum;
for (i=0; i < VECLEN; i++)  заполнение массивов
 
t1 = omp_get_wtime();   начальный замер времени
#pragma omp parallel
sum = dotprod();
t2 = omp_get_wtime();
 
cout<<"Sum = "<<sum<<endl;
cout<<"time: "<<t2 - t1<<endl;
 
system("pause");
}

s += a[i] * b[i];