Распараллеливание алгоритма Беллмана-Форда

#include <iostream>
#include <math.h>
#include <ctime>
#include <omp.h>
 
#ifndef _OPENMP
static_assert(false, "openmp support required");
#endif
 
using namespace std;
 
const int INF = INFINITY;
const int Vmax = 10000;
const int Emax = Vmax * (Vmax - 1) / 2;
 
struct edges {
    int u, v, w;
};
 
int i, j, start_peak;
int peak = 0;
int e = 0;
edges edge[Emax];
int d[Vmax];
 
int num_threads;
 
//главная функция
int main()
{
    srand(time(0));
    setlocale(LC_ALL, "Rus");
    int weight;
 
    cout << "Количество вершин > "; cin >> peak;
 
 
    //Заполнение графа
    for (i = 0; i < peak; i++) {
        for (j = 0; j < peak; j++)
        {
            //cout << "Вес " << i + 1 << "->" << j + 1 << " > "; 
            weight = rand() % 100; //cout << w << endl;
            //cin >> w;
            if (weight != 0)
            {
                edge[e].v = i;
                edge[e].u = j;
                edge[e].w = weight;
                e++;
            }
        }
    }
 
    cout << "Стартовая вершина > "; cin >> start_peak;
 
    cout << "Кол-во потоков: "; cin >> num_threads; cout << endl;
 
    for (int i = 0; i < peak; i++) {
        d[i] = INF;
        d[start_peak] = 0;
    }
 
    //num_threads = 3;
 
    int start_t = omp_get_wtime();
    omp_set_dynamic(0);
    omp_set_num_threads(num_threads);
    #pragma omp parallel
    {
    #pragma omp for schedule (static)   
        for (int i = 0; i < peak - 1; i++) {
            for (int j = 0; j < e; j++) {
                if (d[edge[j].v] + edge[j].w < d[edge[j].u]) {
                    d[edge[j].u] = d[edge[j].v] + edge[j].w;
                }
            }
        }
        int end_t = omp_get_wtime();
        cout << endl << "Время расчета: " << 1000 * (end_t - start_t) << " мс"; cout << endl;
 
        //Вывод кратчайших путей
        //cout << "Список кратчайших путей:";
        //for (i = 0; i < n; i++) if (d[i] == inf)
            //cout << endl << start_peak << "->" << i + 1 << "=" << "Not";
        //else cout << endl << start_peak << "->" << i + 1 << "=" << d[i];
    }
 
    system("pause");
    return 0;
}