Отправляет email-рассылки с помощью сервиса Sendsay
  Все выпуски  

RFpro.ru: Программирование на C / C++


Хостинг портала RFpro.ru:
Московский хостер
Профессиональный ХОСТИНГ на базе Linux x64 и Windows x64

РАССЫЛКИ ПОРТАЛА RFPRO.RU

Лучшие эксперты по данной тематике

Асмик Гаряка
Статус: Советник
Рейтинг: 11043
∙ повысить рейтинг »
Коцюрбенко Алексей aka Жерар
Статус: Советник
Рейтинг: 4067
∙ повысить рейтинг »
CradleA
Статус: Бакалавр
Рейтинг: 2150
∙ повысить рейтинг »

/ КОМПЬЮТЕРЫ И СОФТ / Программирование / C/C++

Номер выпуска:1777
Дата выхода:19.01.2013, 22:00
Администратор рассылки:Киселёва Алёна aka Verena (Академик)
Подписчиков / экспертов:96 / 67
Вопросов / ответов:1 / 1

Консультация # 187104: Здравствуйте! Прошу помощи в следующем вопросе: нужно изменить код программы под другую формулу

Код :
// OpenMP-Example.cpp
#include &q
...

Консультация # 187104:

Здравствуйте! Прошу помощи в следующем вопросе: нужно изменить код программы под другую формулу

Код :
// OpenMP-Example.cpp
#include "stdafx.h"
#include <omp.h>
#include <stdlib.h>
#include <iostream>
#include <stdio.h>
//число разбиений области интегрирования
const int N = 100000000;
//подынтегральная функция
double f (double x)
{
return log(1.0/x);
}
int _tmain(int argc, _TCHAR* argv[])
{
//включаем поддержку русского языка в консоли
setlocale(LC_ALL, "Russian");
//диапазон интегррования
double a = 0.000001, b = 1.0;
//отключение выбора количества потоков (нитей) по умолчанию
omp_set_dynamic(false);
//заданние числа используемых потоков
omp_set_num_threads(2);
//определение типа переменных блока распараллеливания
int nTheads, theadNum;
#pragma omp parallel private(nTheads, theadNum)
{
nTheads = omp_get_num_threads();
18
theadNum = omp_get_thread_num();
//вывод на экран номера и количества используемых потоков
printf("OpenMP поток №%d из %d потоков\n", theadNum, nTheads);
}
//переменные для хранения суммы, значения аргумента и счетчик цикла
double sum = 0.0, x;
int i;
//сохраняем время начала теста
double startTime = omp_get_wtime();
//начало области распараллеливания цикла переменная i -
//приватная, переменная sum обновляется всеми потоками
#pragma omp parallel for default(shared)\
private(i) schedule(static,10)\
reduction (+:sum)
for (i = 0; i < N; i++)
{
x = a + (b-a)*(i-0.5)/N;
sum += f(x);
}
//сохраняем время окончания теста
double endTime = omp_get_wtime();
//выводим сообщения об окончании работы потоками
#pragma omp parallel private(theadNum)
{
theadNum = omp_get_thread_num();
printf("Поток №%d окончил вычисления... \n", theadNum);
}
//вычисление итогового результата и печать на экран
sum = sum/N;
printf("Результат = %f (Ошибка вычисления: %f)\n", sum, 1.0 - sum);
//выводим время работы теста на экран
printf("Время вычисления %f\n", endTime - startTime);
return 0;
}

Дата отправки: 16.01.2013, 14:51
Вопрос задал: Евгений (Посетитель)
Всего ответов: 1
Страница онлайн-консультации »


Консультирует Лысков Игорь Витальевич (Старший модератор):

Здравствуйте, Евгений!
Изменения:
1) добавил #include <math.h> для использования математических функций sin(), exp()
2) изменил функцию f(x)
3) изменил пределы интегрирования a и b
4) в выводе убрал вывод погрешности, который для данной функции не соответствует действительности

Код :
#include "stdafx.h"
#include <omp.h>
#include <stdlib.h>
#include <iostream>
#include <stdio.h>
#include <math.h>
//число разбиений области интегрирования
const int N = 100000000;
//подынтегральная функция
double f (double x)
{
	return exp(-x*x+0.38)/(2.+sin(1./(1.5+x*x)));
}
int _tmain(int argc, _TCHAR* argv[])
{
	//включаем поддержку русского языка в консоли
	setlocale(LC_ALL, "Russian");
	//диапазон интегррования
	double a = 0.4, b = 1.0;
	//отключение выбора количества потоков (нитей) по умолчанию
	omp_set_dynamic(false);
	//заданние числа используемых потоков
	omp_set_num_threads(2);
	//определение типа переменных блока распараллеливания
	int nTheads, theadNum;
	#pragma omp parallel private(nTheads, theadNum)
	{
		nTheads = omp_get_num_threads();
		theadNum = omp_get_thread_num();
		//вывод на экран номера и количества используемых потоков
		printf("OpenMP поток №%d из %d потоков\n", theadNum, nTheads);
	}
	//переменные для хранения суммы, значения аргумента и счетчик цикла
	double sum = 0.0, x;
	int i;
	//сохраняем время начала теста
	double startTime = omp_get_wtime();
	//начало области распараллеливания цикла переменная i -
	//приватная, переменная sum обновляется всеми потоками
	#pragma omp parallel for default(shared)\
		private(i) schedule(static,10)\
		reduction (+:sum)
	for (i = 0; i < N; i++)
	{
		x = a + (b-a)*(i-0.5)/N;
		sum += f(x);
	}
	//сохраняем время окончания теста
	double endTime = omp_get_wtime();
	//выводим сообщения об окончании работы потоками
	#pragma omp parallel private(theadNum)
	{
		theadNum = omp_get_thread_num();
		printf("Поток №%d окончил вычисления... \n", theadNum);
	}
	//вычисление итогового результата и печать на экран
	sum = sum/N;
	printf("Результат = %f\n", sum);
	//выводим время работы теста на экран
	printf("Время вычисления %f\n", endTime - startTime);
	return 0;
}

Консультировал: Лысков Игорь Витальевич (Старший модератор)
Дата отправки: 17.01.2013, 13:11
Рейтинг ответа:

НЕ одобряю 0 одобряю!


Оценить выпуск | Задать вопрос экспертам

главная страница  |  стать участником  |  получить консультацию
техническая поддержка  |  восстановить логин/пароль

Дорогой читатель!
Команда портала RFPRO.RU благодарит Вас за то, что Вы пользуетесь нашими услугами. Вы только что прочли очередной выпуск рассылки. Мы старались. Пожалуйста, оцените его. Если совет помог Вам, если Вам понравился ответ, Вы можете поблагодарить автора - для этого в каждом ответе есть специальные ссылки. Вы можете оставить отзыв о работе портале. Нам очень важно знать Ваше мнение. Вы можете поближе познакомиться с жизнью портала, посетив наш форум, почитав журнал, который издают наши эксперты. Если у Вас есть желание помочь людям, поделиться своими знаниями, Вы можете зарегистрироваться экспертом. Заходите - у нас интересно!
МЫ РАБОТАЕМ ДЛЯ ВАС!



В избранное