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

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

Войти
Регистрация
Восстановить пароль
 
Joe_Hopkins
0 / 0 / 0
Регистрация: 31.05.2013
Сообщений: 2
#1

Задача с openMP - C++

20.12.2013, 14:29. Просмотров 291. Ответов 0
Метки нет (Все метки)

Здравствуйте, я "написал" программу для решения интеграла по методу трапеций, но когда указываю количество нитей(num_threads) больше 1, то ответ становится не правильный, хотя для 1 нити все решает правильно..
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
#include "stdafx.h"
#include <stdio.h>
#include <math.h>
#include <locale.h>
#include <iostream>
#include <omp.h>
#include <ctime>
#include <cstdlib>
 
#define n 100000000
 
double Integral(double x)
{
    double F;
    F = exp(x); // Сюда задается подиннтегральная функция
    return F;
}
 
int main()
{
    omp_set_dynamic(0);
    omp_set_num_threads(1);
    setlocale(LC_ALL, "Rus");
    double a, b, h, S=0, x, eps;
 
    while (true)
    {
        puts("Введите диапазон [a;b] интегрирования.\n");
        printf("a = ");
        scanf_s("%lf", &a);
        printf("b = ");
        scanf_s("%lf", &b);
        if (a>b)
            printf("Неверно заданы пределы интегрирования!\n");
        else
            break;
    }
    h = (b - a) / n;
    clock_t time = clock();
    //S = (Integral(a) + Integral(b))/2;
    
    int i; double pos;
 
#pragma omp parallel shared(S,h,x,a,b) private(i,pos)
 
    {
#pragma omp for
 
        for (i = 1; i <= n; i++)
        {
            //pos = a + i*h;
            S += Integral(i*h + a) + Integral((i + 1)*h + a);
        }
    }
        S = S*h/2;
    
    time = clock() - time;
    printf("Интеграл f(x) равен = %lf ", S);
    printf("\nВремя работы равно %f\n", (double)time / CLOCKS_PER_SEC);
    //printf("\nНажмите любую клавишу для выхода из программы");
    system("pause");
    return 0;
}
Может кто знает в чем может быть проблема. Заранее спасибо
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.12.2013, 14:29
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Задача с openMP (C++):

OpenMP. Время выполнения программы больше чем без OpenMP - C++
Сегодня первый раз сел за OpenMP. Читаю на сайте майкрософта как работает этот API. Так вот там сказано:&quot;Директива #pragma omp for...

OpenMP - C++
есть код в общем виде... #pragma omp parallel for for (i = 0; i &lt; n; ++i) { temp = B; // обращаемся for (j = 0; j &lt; k; ++j)...

OpenMP - C++
#pragma omp parallel private(i) shared(j) { #pragma omp for for (j = 0; j &lt; n; j++) { for (i = 0; i &lt; m;...

OpenMP цикл for - C++
Здравствуйте. Хочу что бы каждый поток вывел по 2 случайных числа, используя цикл for. Вот программа: #include &lt;iostream&gt; using...

Вопросы по OpenMP - C++
Всем доброго времени суток! 1. Заметил, что потоки создаваемые в цикле for так и висят в диспетчере задач Windows до завершения работы...

Многопоточность. OpenMP - C++
Друзья! Есть небольшая задача. У меня есть программка, которая выводит просто какой-либо текст сообразно количеству потоков на компьютере....

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

Цикл for и OpenMP - C++
Недавно распаралелил свой цикл for, а теперь решил посмотреть на сколько я выигрываю от этого. Но вот неожиданность, все многопоточные...

Магия OpenMP - C++
Есть такой код: Expr = ... // vector&lt;pair&lt;size_t, double&gt;&gt; #ifdef _OPENMP #pragma omp parallel num_threads( 2==(size_level-ptr) ? 1...

OpenMP и SIMD - C++
Добрый день. Не могу исправить код, чтобы использовалось 4 ядра, а не 1. В функцию proizv нужно добавить распараллеливание вычислений с...

Программирование на OpenMP - C++
Такой вопрос, у меня есть прога, которая работает ~сутки, не буду вдаваться что и как она делает, мне посоветовали для быстроты...


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

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

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