Гергель В.П. Общий курс Теория и практика параллельных вычислений Лекция 16 Методы разработки параллельных программ для многопроцессорных систем с общей. - презентация

1 Гергель В.П. Общий курс Теория и практика параллельных вычислений Лекция 16 Методы разработки параллельных программ для многопроцессорных систем с общей памятью (стандарт OpenMP) – 2 Нижегородский Государственный Университет им. Н.И. Лобачевского

2 Параллельные Гергель В.П Содержание Директивы OpenMPДирективы OpenMP –Синхронизация Директивы master, critical, barrier, atomic, flush, ordered, threadprivate –Управление областью видимости данных Параметры директив shared, private, firstprivate, lastprivate, copyin, default, reduction –Совместимость директив и их параметров Библиотека функций OpenMPБиблиотека функций OpenMP –Функции для контроля/запроса параметров среды исполнения –Функции синхронизации Переменные среды Реализации OpenMPРеализации OpenMP

3 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива master определяет фрагмент кода, который должен быть выполнен только основным потоком; все остальные потоки пропускают данный фрагмент кода (завершение директивы по умолчанию не синхронизируется) #pragma omp master newline structured_block

4 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива critical определяет фрагмент кода, который должен выполняться только одним потоком в каждый текущий момент времени (критическая секция) основным потоком; все остальные потоки пропускают данный фрагмент кода (завершение директивы по умолчанию не синхронизируется) #pragma omp critical [ name ] newline structured_block

5 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива critical (пример) #include main() { int x; x = 0; #pragma omp parallel shared(x) { #pragma omp critical x = x + 1; } /* end of parallel section */ }

6 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива barrier – определяет точку синхронизации, которую должны достигнуть все процессы для продолжения вычислений (директива должны быть вложена в блок) #pragma omp barrier newline

7 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива atomic – определяет переменную, доступ к которой (чтение/запись) должна быть выполнена как неделимая операция #pragma omp atomic newline statement_expression Возможный формат записи выражения x binop = expr, x++, ++x, x--, --x x должна быть скалярной переменной expr не должно ссылаться на x binop должна быть неперегруженной операцией вида +, -, *, /, &, ^, |, >>,

8 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива flush – определяет точку синхронизации, в которой системой должно быть обеспечено единое для всех процессов состояние памяти (т.е. если потоком какое-либо значение извлекалось из памяти для модификации, измененное значение обязательно должно быть записано в общую память) #pragma omp flush (list) newline Если указан список list, то восстанавливаются только указанные переменные Директива flush неявным образом присутствует в директивах barrier, critical, ordered, parallel, for, sections, single

9 Параллельные Гергель В.П Директивы OpenMP Синхронизация… Директива ordered – указывает фрагмент кода параллельного цикла, который должен выполняться точно в таком же порядке, как и при последовательном выполнении #pragma omp ordered newline structured_block В каждый момент времени в блоке ordered может находиться только один поток На одной итерации цикла может быть только одна директива ordered и эта директива может выполниться только однократно Цикл, в котором имеется директива ordered, должен иметь параметр ordered

10 Параллельные Гергель В.П Директивы OpenMP Синхронизация Директива threadprivate – используется для создания поточных копий для глобальных переменных программы; созданные копии не видимы между потоками, но существуют во все время выполнения программы #pragma omp threadprivate (list)

11 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Общие (разделяемые между потоками, shared) переменные – static, переменные с областью видимости в пределах файла Локальные (private) данные потоков - переменные циклы Управление областью видимости обеспечивается при помощи параметров (clause) директив private, firstprivate, lastprivate, shared, default, reduction, copyin которые определяют, какие соотношения существуют между переменными последовательных и параллельных фрагментов выполняемой программы

12 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Параметр shared определяет список переменных, которые будут общими для всех потоков параллельной области; правильность использования таких переменных должна обеспечиваться программистом shared (list) Параметр private определяет список переменных, которые будут локальными для каждого потока; переменные создаются в момент формирования потоков параллельной области; начальное значение переменных является неопределенным private (list)

13 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Параметр firstprivate позволяет создать локальные переменные потоков, которые перед использованием инициализируются значениями исходных переменных firstprivate (list) Параметр lastprivate позволяет создать локальные переменные потоков, значения которых запоминаются в исходных переменных после завершения параллельной области (используются значения потока, выполнившего последнюю итерацию цикла или последнюю секцию) lastprivate (list)

14 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Параметр copyin позволяет выполнить инициализацию переменных директивы threadprivate copyin (list) Параметр default устанавливает область видимости переменных по умолчанию default (shared | none)

15 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Параметр reduction определяет список переменных, для которых выполняется операция редукции; перед выполнением параллельной области для каждого потока создаются копии этих переменных, потоки формируют значения в своих локальных переменных и при завершении параллельной области на всеми локальными значениями выполняются необходимые операции редукции, результаты которых запоминаются в исходных (глобальных) переменных reduction (operator: list)

16 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных… Параметр reduction (правила записи) Возможный формат записи выражения x = x op expr x = expr op x x binop = expr x++, ++x, x--, --x x должна быть скалярной переменной expr не должно ссылаться на x op (operator) должна быть неперегруженной операцией вида +, -, *, /, &, ^, |, &&, || binop должна быть неперегруженной операцией вида +, -, *, /, &, ^, |

17 Параллельные Гергель В.П Директивы OpenMP Управление областью видимости данных Параметр reduction (пример) #include main () { /* vector dot product */ int i, n, chunk; float a[100], b[100], result; /* Some initializations */ n = 100; chunk = 10; result = 0.0; for (i=0; i < n; i++) { a[i] = i * 1.0; b[i] = i * 2.0; } #pragma omp parallel for \ default(shared) private(i) \ schedule(static,chunk) \ reduction(+:result) for (i=0; i < n; i++) result = result + (a[i] * b[i]); printf("Final result= %f\n",result); }

18 Параллельные Гергель В.П Директивы OpenMP Совместимость директив и их параметров

19 Параллельные Гергель В.П void omp_set_num_threads(int num_threads) Позволяет назначить максимальное число потоков для использования в следующей параллельной области (если это число разрешено менять динамически). Вызывается из последовательной области программы. int omp_get_max_threads(void) Возвращает максимальное число потоков. int omp_get_num_threads(void) Возвращает фактическое число потоков в параллельной области программы. Библиотека функций OpenMP Функции для контроля/запроса параметров среды исполнения…

20 Параллельные Гергель В.П Функции для контроля/запроса параметров среды исполнения int omp_get_thread_num(void) Возвращает номер потока. int omp_get_num_procs(void) Возвращает число процессоров, доступных приложению. int omp_in_parallel(void) Возвращает.TRUE., если вызвана из параллельной области программы. void omp_set_dynamic(int dynamic) int omp_get_dynamic(void) Устанавливает/запрашивает состояние флага, разрешающего динамически изменять число потоков. void omp_get_nested(int nested) int omp_set_nested(void) Устанавливает/запрашивает состояние флага, разрешающего вложенный параллелизм. Функции для контроля/запроса параметров среды исполнения Библиотека функций OpenMP

21 Параллельные Гергель В.П Функции синхронизации… В качестве замков используются общие переменные типа omp_lock_t или omp_nestlock_t. Данные переменные должны использоваться только как параметры примитивов синхронизации. void omp_init_lock(omp_lock_t *lock)void omp_nest_init_lock(omp_nest_lock_t *lock) Инициализирует замок, связанный с переменной lock. void omp_destroy_lock(omp_lock_t *lock) void omp_destroy_nest__lock(omp_nest_lock_t *lock) Удаляет замок, связанный с переменной lock. Библиотека функций OpenMP

22 Параллельные Гергель В.П Функции синхронизации void omp_set_lock(omp_lock_t *lock) void omp_set_nest__lock(omp_nest_lock_t *lock) Заставляет вызвавший поток дождаться освобождения замка, а затем захватывает его. void omp_unset_lock(omp_lock_t *lock) void omp_unset_nest__lock(omp_nest_lock_t *lock) Освобождает замок, если он был захвачен потоком ранее. void omp_test_lock(omp_lock_t *lock) void omp_test_nest__lock(omp_nest_lock_t *lock) Пробует захватить указанный замок. Если это невозможно, возвращает.FALSE. Библиотека функций OpenMP

23 Параллельные Гергель В.П Переменные среды OMP_SCHEDULE Определяет способ распределения итераций в цикле, если в директиве DO использована клауза SCHEDULE(RUNTIME). OMP_NUM_THREADS Определяет число нитей для исполнения параллельных областей приложения. OMP_DYNAMIC Разрешает или запрещает динамическое изменение числа нитей. OMP_NESTED Разрешает или запрещает вложенный параллелизм. Компилятор с поддержкой OpenMP определяет макрос "_OPENMP", который может использоваться для условной компиляции отдельных блоков, характерных для параллельной версии программы

24 Параллельные Гергель В.П Пример #include #define THREADNUMS 2 main () { /* вычисление числа */ long StepNums = 10000; double step, x, pi, sum=0.0; int i; step = 1.0/StepNums; omp_set_num_threads(THREADNUMS); #pragma omp parallel for reduction(+:sun) private(i,x) for (i=0; i

25 Параллельные Гергель В.П Сравнение технологий MPI и OpenMP для систем с общей памятью… Сравнение времени выполнения последовательного варианта программы для задачи матричного умножения с вариантами OpenMP и MPI для 2- процессорного сервера

26 Параллельные Гергель В.П Сравнение технологий MPI и OpenMP для систем с общей памятью… Ускорение матричного умножения при использовании параллельных вычислений

27 Параллельные Гергель В.П Сравнение технологий MPI и OpenMP для систем с общей памятью… Сравнение времени выполнения последовательного варианта программы для задачи матричного умножения с вариантами OpenMP и MPI для 4- процессорного сервера

28 Параллельные Гергель В.П Сравнение технологий MPI и OpenMP для систем с общей памятью Ускорение матричного умножения при использовании параллельных вычислений

29 Параллельные Гергель В.П Комбинированная (MPI+OpenMP) технология программирования для систем с общей памятью… Сравнение времени выполнения MPI и MPI+OpenMP вариантов программы для двух 2-процессорных серверов

30 Параллельные Гергель В.П Комбинированная (MPI+OpenMP) технология программирования для систем с общей памятью Ускорение матричного умножения при использовании смешанного MPI+OpenMP варианта параллельной программы

31 Параллельные Гергель В.П Реализации OpenMP 1) Silicon Graphics. Fortran 77/90 (IRIX), планируется поддержка OpenMP для C/C++. 2) Compaq/DEC. DIGITAL Fortran 3) Kuck & Associates (KAI). Fortran, C/C++ (Unix, Windows) 4) Portland Group (PGI). Fortran и C/C++ для Windows NT, Linux, Solaris (x86). 5) OdinMP. OpenMP-препроцессор для языка С, генерация программы в стандарте POSIX threads.OdinMP 6) Sun. Планируется поддержка OpenMPSun 7) Pacific-Sierra Research предлагает распараллеливающие препроцессоры VAST/Parallel для Fortran и С, которые обеспечивают автоматическое распознавание параллелизма в программах и выполнение необходимой трансформации программ путем добавления соответствующих директив OpenMP.Pacific-Sierra Research 8) Intel. Fortran, C/C++ (Unix, Windows)

32 Параллельные Гергель В.П Информационные ресурсы Что такое OpenMP OpenMP C/C++ specification v Introduction to OpenMP - MAIN.html MAIN.html Chandra, R., Menon, R., Dagum, L., Kohr, D., Maydan, D., McDonald, J. (2000). Parallel Programming in OpenMP. Morgan Kaufmann Publishers.

33 Параллельные Гергель В.П Вопросы для обсуждения Методы синхронизации обработки данных в OpenMP Программирование с использованием библиотеки функций OpenMP

34 Параллельные Гергель В.П Задания для самостоятельной работы Разработка параллельных методов для задач линейной алгебры при использовании интерфейса OpenMP

35 Параллельные Гергель В.П Заключение Методы синхронизации обработки данных Библиотека функций OpenMP Переменные среды окружения