VPF::Распараллеливание с помощью openmp

shonny

Дата 14.11.2017, 20:03 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 1
Регистрация: 14.11.2017

Репутация: нет
Всего: нет

Здравствуйте, вообщем пытаюсь использовать Openmp, чтобы распараллелить задачу вычисления экспоненты с помощью ряда Тейлора. Программа работает, но при сравнении времени выполнения последовательного кода и с Openmp, нет прироста скорости, а даже наоборот выполняется дольше с Openmp. В чем может быть причина, не могли бы вы помочь?

Код


#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <omp.h>
int fact(int c);
double taylor(double x, int n);
 
int main(int argc, char *argv[]) {
  double x=0;
  int n=0;
  int n_max = 32;
  if (argc != 3) {
  printf("Input three parametrs!\n");
  return 1;
}
  sscanf(argv[1], "%lf", &x);
  sscanf(argv[2], "%d", &n);
 
 if((n>n_max)||(n<0)) {
 printf("Invalid parameter n, enter n from 0..10\n");
 return 1;
 }
 
 int i;
 double sum=0 ;
 double t = omp_get_wtime();
 #pragma omp parallel for private(i) reduction(+: sum)
   for (i=1; i<=n; i++) {
      double taylor = pow(x,i)/fact(i);
      sum+=taylor;
  }
 
  double diff = omp_get_wtime() - t;
  printf("e^x = %lf\n Time: \t %lf\n", ++sum, diff);
  return 0;
}
 
int fact(int c) {
   int fact=1;
   int i;
   for (i = 1; i <= c; i++ ) {
      fact *=i;
   }
   return fact;
}

./taylor 1 32
e^x = 2.718282
Time: 0.000126 - время выполнения с openmp

/staylor 1 32
e^x = 2.718282
Time: 0.000020 - время выполнения последовательного кода

xvr

Дата 14.11.2017, 20:51 (ссылка)

(голосов:1)

Загрузка ...

Эксперт

Профиль
Группа: Комодератор
Сообщений: 7046
Регистрация: 28.8.2007
Где: Дублин, Ирландия

Репутация: 60
Всего: 223

Распараллеливание цикла на 32 привара не даст. Количество итераций должно быть больше порядков на 5-6
У вас накладные расходы на создание потоков больше, чем работа, которая в них проводится

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »