Создать параллельную версию одномерного массива

    std::string product[]{"milk", "butter", "bread"};
    int price[]{100, 200, 50};

#include <iostream>
#include <cmath>
#include <ostream>
 
 
double theta(double a, double b, double epsilon) {
    return (epsilon - std::fabs(a - b) < 0) ? 0 : 1;
}
 
template<std::size_t N, typename Function>
double correlationSum(double (&system)[N], Function f, double epsilon) {
    double sum = 0;
    for (std::size_t i = 0; i < N; ++i) {
        for (std::size_t j = 0; j < N; ++j) {
            sum += f(system[i], system[j], epsilon);
        }
    }
    return sum / (N * N);
}
 
int main() {
 
    constexpr std::size_t N = 10;
    double system[N] = { 3, 4, 5, 3, 4, 5, 1, 3, 1, 7 };
    double epsilon = 0;
 
    std::cout << correlationSum(system, theta, epsilon) << std::endl;
 
    return 0;
}

#include <iostream>
#include <cmath>
#include <ostream>
#include <cstdlib>
#include <Windows.h>
#include <conio.h>
#include <stdlib.h>
#include "omp.h"
#include <ctime>
 
using namespace std;
 
double theta(double a, double b, double epsilon)
{
    return (epsilon - fabs(a - b) <= 0) ? 0 : 1;
}
 
template<size_t N, typename Function>
double CorSum(double(&Num)[N], Function f, double epsilon)
{
    int P;
    cout << "Кол-во потоков: P= "; cin >> P;
    int i, j;
    double sum = 0;
#pragma omp parallel num_threads(P) reduction(+: sum)
    {
# pragma omp for 
        for (i = 0; i <= N - 1; i++)
        {
            for (j = 0; j <= N - 1; j++)
            {
                sum += f(Num[i], Num[j], epsilon);
                //cout << sum<<endl;
            }
        }
    }
    return sum / (N * N);
}
 
int main()
{
    setlocale(0, "");
    srand(time(NULL));
    float fTimeStart1 = clock() / (float)CLOCKS_PER_SEC;
    const size_t B = 5000;
    double s[B] = { rand() % 10 };
    double epsilon = 4;
    cout << CorSum(s, theta, epsilon) << endl;
    float fTimeStop = clock() / (float)CLOCKS_PER_SEC;
    printf("Время работы %f с.\n", fTimeStop - fTimeStart1);
    cout << endl;
    system("pause");
    return 0;
}