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

Параллельное программирование: нахождение интеграла - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 9, средняя оценка - 4.56
senna
 Аватар для senna
0 / 0 / 0
Регистрация: 14.05.2012
Сообщений: 19
19.01.2013, 22:17     Параллельное программирование: нахождение интеграла #1
Нужно распараллелить метод средних прямоугольников.

Долго не мучаясь, решил воспользоваться OpenMP, вот что получилось:

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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
// intagral_omp.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include <omp.h>
#include <iostream>
#include <time.h>
#include <math.h>
#include <conio.h>
using namespace std;
#define n 2000
 
 
double f(double x)
{
    return (x*x)*4*pow(sin(x),5);
}
    
int _tmain(int argc, _TCHAR* argv[])
{
    omp_set_dynamic(0);
    omp_set_num_threads(4);
    double a=-5.0;
    double b=10.0;
    double h;
    double* x=new double[n];
    double result=0;
    int thr=0;
 
 
    x[0]=a;
    x[n]=b;
    h=fabs(a-b)/n;
    
    cout<<"Integral calculator"<<endl;
    cout<<"\n\na = "<<a<<endl;
    cout<<"\nb = "<<b<<endl;
    cout<<"\nh = "<<h<<endl;
    cout<<"\n\n Number of divisions = "<<n<<endl;
    for(int i=0;i<=n;i++)
    {
        x[i]=a+i*h;
        //cout<<x[i]<<endl;
    }
    
#pragma omp parallel for reduction(+: result) 
    for(int i=1;i<=n;i++)
    {
        result+=h*f((x[i-1]+x[i])/2);
        if(omp_get_num_threads()>=thr)
        {
        thr=omp_get_num_threads();
        }
 
    }
    cout<<"\n\nMax number of used threads: "<<thr<<endl;
    cout<<result;
    //system("pause");
    _getch();
 
 
    return 0;
}
Но проблема в том, что при вычислении используется только один поток (определил с помощью omp_get_num_threads()). Что не так??

Спасибо
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
19.01.2013, 22:17     Параллельное программирование: нахождение интеграла
Посмотрите здесь:

параллельное программирование C++
Параллельное программирование C++
C++ Параллельное программирование
Параллельное программирование в С++ C++
параллельное программирование C++
параллельное программирование C++
Параллельное программирование C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
DU
1477 / 1053 / 45
Регистрация: 05.12.2011
Сообщений: 2,279
19.01.2013, 22:35     Параллельное программирование: нахождение интеграла #2
поддержка openmp включена при сборке проекта?
senna
 Аватар для senna
0 / 0 / 0
Регистрация: 14.05.2012
Сообщений: 19
19.01.2013, 23:20  [ТС]     Параллельное программирование: нахождение интеграла #3
DU, о боже как я ступил))) конечно она выключена) спасибо, что заглянул сюда

Ну ладно, может кому нибудь мой код пригодится.
Yandex
Объявления
19.01.2013, 23:20     Параллельное программирование: нахождение интеграла
Ответ Создать тему
Опции темы

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