Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.83/6: Рейтинг темы: голосов - 6, средняя оценка - 4.83
0 / 0 / 0
Регистрация: 20.12.2018
Сообщений: 1
1

Параллельные вычисления: произведение двух векторов (OpenMP)

20.12.2018, 23:54. Просмотров 1153. Ответов 2
Метки нет (Все метки)

Помогите с написанием программы вычисления произведения двух векторов с использованием параллельного алгоритма OpenMP)

Добавлено через 4 минуты
КАК ЭТО РАСПАРАЛЛЕЛИТЬ ?





C++ (Qt)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
//Scalar.cpp
#include <stdlib.h>
#include <iostream>
#include <vector>
#include <cmath>
 
using namespace std;
 
/**
This function returns the scalar product of two vectors "a" and "b"
*/
double scalar_product(vector<double> a, vector<double> b)
{
    //In C++, you should declare every variable before you use it. So, you declare product and initialize it to 0.
    double product = 0;
    //Here you check whether the two vectors are of equal size. If they are not then the vectors cannot be multiplied for scalar product.
    if(a.size()!=b.size()){
        cout << "Vectors are not of the same size and hence the scalar product cannot be calculated" << endl;
        return -1;  //Note: This -1 is not the answer, but just a number indicating that the product is not possible. Some pair of vectors might actually have a -1, but in that case you will not see the error above.
    }
 
    //you loop through the vectors. As bobo also pointed you do not need two loops.
    #pragma omp parallel for
    for (int i = 0; i < a.size(); i++)
    {
        product = product + a[i]*b[i];
    }
 
    //finally you return the product
    return product;
}
 
 
 //This is your main function that will be executed before anything else.
int main() {
    //you declare two vectors "veca" and "vecb" of length 2 each
    vector<double> veca(2);
    vector<double> vecb(2);
 
    //put some random values into the vectors
    veca[0] = 1.5;
    veca[1] = .7;
    vecb[0] = 1.0;
    vecb[1] = .7;
 
    //This is important! You called the function you just defined above with the two parameters as "veca" and "vecb". I hope this cout is simple!
    cout << scalar_product(veca,vecb) << endl;
}
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
20.12.2018, 23:54
Ответы с готовыми решениями:

Параллельные вычисления с использованием параллельных секций в OpenMP
Доброго времени суток, ребята! Ребят, нужно вычислить интеграл на отрезке для функции. Нужно...

Перевести программу с паскаля на фортран, используя Параллельные вычисления и OpenMP
Есть несложная программа на паскале, надо ее перевести в фортран с помощью параллельных вычислений...

Параллельные сортировки с OpenMP
Привет, читателю . Надеюсь не поленитесь прочесь до конца, кто понимает =) Вообщем суть такова, у...

OPENmp параллельные потоки QT
#include &lt;iostream&gt; #include &quot;omp.h&quot; using namespace std; int main() { #pragma omp...

2
2112 / 1653 / 824
Регистрация: 21.12.2010
Сообщений: 2,926
Записей в блоге: 12
21.12.2018, 10:26 2
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
#include <iostream>
#include <vector>
#include <chrono>
#include <random>
#include <omp.h>
using namespace std;
 
 
double scalar_product(vector<double> const& a, vector<double> const& b)
{
    double product = 0.;
 
    auto tp1 = std::chrono::high_resolution_clock::now();
 
    omp_set_num_threads(8);
#pragma omp parallel for reduction(+:product)
    for (int i = 0; i < a.size(); ++i)
    {
//#pragma omp atomic
        product += (a[i] * b[i]);
    }
 
    auto tp2 = std::chrono::high_resolution_clock::now();
    std::cout << "dur: " << std::chrono::duration_cast<std::chrono::milliseconds>(tp2 - tp1).count() << std::endl; // время выполнения цикла
 
    return product;
}
 
 
int main()
{
    vector<double> va ( 2'000'000 );
    vector<double> vb ( va.size() );
 
    std::mt19937 eng{ std::random_device{}() };
    std::uniform_real_distribution<> dst{-10., 10.};
    omp_set_num_threads(8);
#pragma omp parallel for 
    for (int i = 0; i < va.size(); ++i) // заполнение векторов рандомными double из интервала -10, 10
    {
        // va[i] = 1.; vb[i] = 1.; // for testing
 
        va[i] = dst(eng); vb[i] = dst(eng);
    }
 
    cout << "prod: " << scalar_product(va, vb) << endl;
}
0
647 / 269 / 92
Регистрация: 04.07.2014
Сообщений: 757
24.12.2018, 13:13 3
Цитата Сообщение от igorrr37 Посмотреть сообщение
C++
1
2
3
4
5
6
7
8
9
  std::mt19937 eng{ std::random_device{}() };
  std::uniform_real_distribution<> dst{-10., 10.};
  omp_set_num_threads(8);
#pragma omp parallel for 
  for (int i = 0; i < va.size(); ++i) // заполнение векторов рандомными double из интервала -10, 10
  {
    // va[i] = 1.; vb[i] = 1.; // for testing
va[i] = dst(eng); vb[i] = dst(eng);
  }
Так у тебя объекты dst и eng общие для нескольких потоков!!!

C++
1
2
3
4
5
6
7
8
9
10
  #pragma omp parallel
  {
    std::mt19937 eng{ std::random_device{}() };
    std::uniform_real_distribution<> dst{-10., 10.};
    #pragma omp for 
    for (int i = 0; i < va.size(); ++i) // заполнение векторов рандомными double из интервала -10, 10
    {
      va[i] = dst(eng); vb[i] = dst(eng);
    }
  }
1
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
24.12.2018, 13:13

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь.

Вычисления скалярного произведение N-мерных векторов
Помогите написать алгоритм вычисления скалярного произведение N-мерных векторов

Скалярное произведение двух векторов
найти скалярное произведение 2 векторов из n элементов введенных с клавиатуры

Скалярное произведение двух векторов
Даны два вектора (одномерных массива), содержащих n вещественных элементов. Найти скалярное ...

Найти скалярное произведение двух векторов
Доброго дня суток. Нужна помощь в решении этой задачи: Найти скалярное произведение a и b У меня...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
3
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.