Параллельное программирование openmp (Умножение матрицы на вектор)

#include "stdafx.h"
#include <iostream>
#include <omp.h>
#include <thread>
#include <ctime>
using namespace std;
 
void show(int **mas, int n, int m) {
 
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < m; j++) {
            cout << mas[i][j] << " ";
        }
        cout << endl;
    }
}
 
 
 
void main() {
    setlocale(LC_ALL, "Russian");
    int m, n;
    cout << "enter n=";
    cin >> n;
    cout << endl;
 
    cout << "enter m=";
    cin >> m;
    cout << endl;
    int **mas;
    mas = new int*[n];
    for (int i = 0; i < n; i++)
        mas[i] = new int[m];
 
    int *vector = new int[m];
 
    int *result = new int[m];
 
 
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < m; j++) {
            mas[i][j] = 0 + rand() % 10;
        }
    }
 
    for (int i = 0; i < m; i++)
    {
        vector[i] = 0 + rand() % 10;
        cout << vector[i] << " ";
    }
 
    cout << endl << endl;
 
    show(mas, n, m);
 
    unsigned int start_time = clock();
    for (int i = 0; i < n; i++)
    {
        result[i] = 0;
        for (int j = 0; j < m; j++)
        {
            cout << mas[i][j] << endl;
            result[i] = result[i] + mas[i][j] * vector[j];
        }
    }
 
    unsigned int end_time = clock();
    unsigned int search_time = end_time - start_time;
 
    cout << "TIME="<<search_time << endl;
 
    for (int i = 0; i < n; i++)
    {
        cout << result[i] << " ";
    }
    
    cout << endl << endl;
 
 
    int *result1 = new int[m];
 
 
    unsigned int start_time1 = clock();
 
    omp_set_num_threads(8);
#pragma omp parallel shared(mas,vector, result1)
    {
#pragma omp for schedule (dynamic,1)
            for (int i = 0; i < n; i++)
            {
                result1[i] = 0;
                for (int j = 0; j < m; j++)
                {
                    printf("ID POTOKA=%d\n", this_thread::get_id());
                    printf("Этот поток вычисляет %d строчку\n", i);
                    result1[i] = result1[i] + mas[i][j] * vector[j];
                }
                printf("result1[%d]=%d\n", i, result1[i]);
            }
    }
    
    unsigned int end_time1 = clock();
    unsigned int search_time1 = end_time1 - start_time1;
    
    cout << "ВРЕМЯ ПАРАЛЛЕЛИ=" << search_time1 << endl;
    
    
    
    system("pause");
 
 
 
 
}

#include <iostream>
#include <omp.h>
#include <thread>
#include <time.h>
 
 
 
 
void show(int **mat, int n, int m) {
 
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < m; ++j) {
            std::cout << mat[i][j] << " ";
        }
        std::cout << '\n';
    }
    std::cout << "\n\n";
}
 
 
void show(int* arr, int size) {
    for (int i = 0; i < size; ++i) {
        std::cout << arr[i] << " ";
    }
    std::cout << "\n\n";
}
 
 
 
 
int main() {
    setlocale(LC_ALL, "Russian");
 
    int m, n;
    int** matrix;
    int* array_;
    int* result;
 
 
    std::cout << "enter n = ";
    std::cin >> n;
 
    std::cout << "enter m = ";
    std::cin >> m;
    std::cout << '\n';
 
    matrix = new int*[n];
    array_ = new int[m];
    result = new int[m];
 
    for (int i = 0; i < n; ++i) {
        matrix[i] = new int[m];
    }
 
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < m; ++j) {
            matrix[i][j] = rand() % 10;
        }
    }
    //show(matrix, n, m);
 
 
    for (int i = 0; i < m; ++i) {
        array_[i] = rand() % 10;
    }
    //show(array_, m);
 
 
    clock_t timeBegin = clock();
 
    for (int i = 0; i < n; ++i)
    {
        result[i] = 0;
        for (int j = 0; j < m; ++j)
        {
            result[i] += matrix[i][j] * array_[j];
        }
    }
    //show(result, n);
 
    std::cout << "Total ticks 1 = " << clock() - timeBegin << '\n';
 
 
 
 
    omp_set_num_threads(8);
    timeBegin = clock();
 
#pragma omp parallel shared(matrix, array_, result)
    {
#pragma omp for schedule (dynamic)
        for (int i = 0; i < n; ++i) {
            result[i] = 0;
            for (int j = 0; j < m; ++j)
            {
                result[i] += matrix[i][j] * array_[j];
            }
        }
    }
 
    std::cout << "Total ticks 2 = " << clock() - timeBegin << '\n';
 
 
 
    delete[] matrix;
    delete[] array_;
    delete[] result;
 
    system("pause");
    return 0;
}
 
 
 
/* output
enter n = 1000
enter m = 2000
 
Total ticks 1 = 12
Total ticks 2 = 8
 
*/

    for (int i = 0; i < n; ++i) {
        delete[] matrix[i];
    }
    delete[] matrix;