Форум программистов, компьютерный форум, киберфорум
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Kill100
406 / 272 / 37
Регистрация: 11.12.2010
Сообщений: 1,157
Завершенные тесты: 1
#1

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

05.12.2012, 16:03. Просмотров 939. Ответов 0
Метки нет (Все метки)

Надо написать вычисление интеграла и распараллелить его.
Я написал левыми прямоугольниками. Но проблема в том что #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;
    }
Можно ли еще как ни будь улучшить?
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
05.12.2012, 16:03
Здравствуйте! Я подобрал для вас темы с ответами на вопрос [OMP] Как распараллелить for (C++):

Как правильно задать количество потоков в #pragma omp for - C++
Пробовал сделать так,но когда питаюсь вывести номер потока с помощью cout &lt;&lt; omp_get_thread_num() по идее должны быть разные номера,но...

Как распараллелить цикл while? - C++
Доброго времени суток. При распараллелировании применять технологии Openmp. Было бы здорово, если бы ответ на вопрос, что в заглавии,...

Как распараллелить обработку массива данных? - C++
Есть массив данных. Есть функция, которая обрабатывает элемент этого массива. Надо распараллелить обработку. Для этого хочу в цикле...

MPI. Как распараллелить алгоритм нахождение нормы матрицы - C++
собственно задача описана в заголовке. разобрал все функции MPI, но как применить к этому алгоритму - хз. везде в примерах легкотня одна. ...

Не могу разобраться с OpenMP! Как правильно распараллелить формулу? - C++
Только начинаю изучать ПП. Помогите разобраться как распараллелить формулу. Данные зависимы. Не могли бы показать на данном примере как...

OMP parallel for - C++
#pragma omp parallel for for (int j = 0; j &lt; whatFind.size(); j++) { string temp = &quot;&quot;; int t = Find(whatFind, 0,...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.12.2012, 16:03
Привет! Вот еще темы с ответами:

Распараллеливание через omp.h - C++
// labaOMP.cpp: определяет точку входа для консольного приложения. // #include &quot;stdafx.h&quot; #include &lt;iostream&gt; #include &lt;omp.h&gt; ...

Проблемы с распараллеливанием(omp) простой программы - C++
Нужно реализовать и распараллелить простенький алгоритм: Найти минимальное число, большее заданного N, которое делится нацело на все...

Почему не работает #pragma omp task? - C++
Сижу быструю сортировку распараллеливаю. Использую подручные средства, а именно Visual Studio 2013 Ultimate. Программа никак не хочет...

Распараллеливание потоков для нахождения суммы матрицы <omp.h> - C++
Доброго дня. Сегодня слушал занимательную лекцию об распараллеливание потоков и получил задание: есть матрица 10,10 инициализированная...


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

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

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