Параллельное умножение матриц с указанием повторов

#include "stdlib.h"
#include "time.h"
#include "iostream"
#include "conio.h"
#include "ctime"
#include <omp.h>
#include <Windows.h> //Библиотека для установки русского языка в консоли
 
using namespace std;
double S;
 
int main()
{
    SetConsoleCP(1251); //Позволяет выводить на экран русские символы
    SetConsoleOutputCP(1251);
    double start, end, times;
    int Q, N;
    double **aMatrix, **bMatrix, **Result;
    cout << "Введите размер матрицы" << endl;
    cin >> N;
    cout << "Введите Q" << endl;
    cin >> Q;
    int i, j, k;
    double a = 0.5, b = -0.5;
    //Ввод матрицы A
    aMatrix = new double*[N];
    for (i = 0; i < N; i++)
    {
        aMatrix[i] = new double[N];
        for (j = 0; j < N; j++)
            aMatrix[i][j] = rand()*(b - a) / RAND_MAX + a;
    }
    //Ввод матрицы B
    bMatrix = new double*[N];
    for (i = 0; i < N; i++)
    {
        bMatrix[i] = new double[N];
        for (j = 0; j < N; j++)
            bMatrix[i][j] = rand()*(b - a) / RAND_MAX + a;
    }
    cout << "Вывод матрицы А" << endl;
    //Вывод матрицы А
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
        {
            cout << aMatrix[i][j] << " ";
        }
        cout << endl;
    }
    cout << "Вывод матрицы B" << endl;
    //Вывод матрицы B
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
        {
            cout << bMatrix[i][j] << " ";
        }
        cout << endl;
    }
    srand(time(NULL));
    cout << "Вывод результата" << endl;
    omp_set_num_threads(4); //указание числа потоков
 
    //Умножение матрицы
 
    //start = clock();
    
    start = omp_get_wtime();
    Result = new double*[N];
    
    for (int l = 0; l < Q; l++)
    {
            for (i = 0; i < N; i++)
            {
                Result[i] = new double[N];
                
                for (j = 0; j < N; j++)
                {
                    //Result[i][j] = 0;
                    S = 0;
                    #pragma omp for reduction (+:S)
                    for (k = 0; k < N; k++)
                    {
                        //Result[i][j] += aMatrix[i][k] * bMatrix[k][j];
                        S += aMatrix[i][k] * bMatrix[k][j];
                    }
                    cout << S << "  ";
                    //cout << Result[i][j] << "  ";
                }
                cout << "\n";
            }
        //end = clock();
        end = omp_get_wtime();
        times = (end - start);
    }
    //cout << "Время=" << times / CLK_TCK << "секунды" << "\n";
    cout << "Время=" << times << "секунды" << "\n";
}