Форум программистов, компьютерный форум CyberForum.ru

[OMP] Как распараллелить for - C++

Восстановить пароль Регистрация
 
Kill100
 Аватар для Kill100
359 / 248 / 33
Регистрация: 11.12.2010
Сообщений: 1,068
Завершенные тесты: 1
05.12.2012, 16:03     [OMP] Как распараллелить for #1
Надо написать вычисление интеграла и распараллелить его.
Я написал левыми прямоугольниками. Но проблема в том что #pragma omp parallel for принимает только for с типом int а у меня double
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
#include "stdafx.h"
#include "omp.h"
#include "conio.h"
#include "string.h"
#include "math.h"
#include <iostream>
#include <vector>
#include <queue>
#include <stack>
#include <time.h>
using namespace std;
 
//сама F(x)
double Fx(double x)
{
    return (0.5 * x + x * log(x));
}
//Численный метод прямоугольников
//H - количество фигур на которое разбиваем
double Integral(double a, double b, double H)
{
    double P = 0;//ответ
    double h = (b - a) / H;//шаг
    omp_set_dynamic(1);//Динамичное распределение потоков
//само рзспараллливание
#pragma omp parallel for reduction(+: P)
    for (double x = a; x < b; x += h)
    {
        P += ((Fx(x) + Fx(x + h)) / 2);
    }
    //printf("%d\n",omp_get_num_threads());
    return P * h;
}
int _tmain(int argc, _TCHAR* argv[])
{
    time_t t1, t2;
    t1 = clock();
    printf("%f\n",Integral(1,2,99999990));//ответ
    t2 = clock();
    double duration = ((double)t2-(double)t1)/double(CLOCKS_PER_SEC);
    printf("%f\n",duration);//время выполнения
    system("pause");
    return 0;
}
Добавлено через 18 часов 53 минуты
Всем спасибо за внимание
В итоге получил.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
 //Расчет интеграла, распараллеливание for + редукция
    double Integral(double a, double b, int H)
    {
        //ответ
        double P = 0;
        //шаг
        double h = (b - a) / H;
        omp_set_dynamic(1);
        omp_set_num_threads(8);
        //параллельный блок
#pragma omp parallel  
        {
#pragma omp for reduction(+: P)
            for (int dd=0; dd < H; dd ++)
            {
                P += ((Fx((dd*h+a)) + Fx((dd+1)*h+a)) / 2);
            }
            //вывод количества испльзуемых потоков
#pragma omp critical
            printf("%d\n",omp_get_num_threads());
        }
        return P * h;
    }
Можно ли еще как ни будь улучшить?
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
05.12.2012, 16:03     [OMP] Как распараллелить for
Посмотрите здесь:

C++ Помогите распараллелить цикл!!!
Проблемы с распараллеливанием(omp) простой программы C++
C++ OpenMP, Visual Studio. выдает ошибку: fatal error C1083: Не удается открыть файл включение: omp.h: No such file or directory
Можно ли это распараллелить? C++
Распараллеливание потоков для нахождения суммы матрицы <omp.h> C++
Расспаралеливание omp.h C++
MPI. Как распараллелить алгоритм нахождение нормы матрицы C++
Как распараллелить цикл while? C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 09:22. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru