5.1.18 OpenMP

OpenMP (Open Multi-Processing) — это интерфейс прикладного программирования (API), который поддерживает многоплатформенное программирование с общей памятью для многопроцессорной обработки в C/C++ и Fortran на многих архитектурах, включая Unix и Microsoft Windows. Он состоит из набора директив компилятора, процедур библиотеки и переменных среды, которые влияют на поведение во время выполнения.

GNU Fortran стремится к совместимости с OpenMP Application Program Interface v4.5.

Для включения обработки директивы OpenMP !$omp в исходном коде свободного формата; директив c$omp, *$omp и !$omp в фиксированном формате; условных sentinels !$ в свободном формате; и sentinels c$, *$ и !$ в фиксированном формате, необходимо вызвать gfortran с -fopenmp. Это также обеспечивает автоматическую линковку библиотеки GNU Offloading и Multi Processing Runtime Library libgomp.

Библиотечные процедуры OpenMP для Fortran предоставляются как в форме модуля Fortran 90 с именем omp_lib, так и в форме файла Fortran include с именем omp_lib.h.

Пример паралелизованного цикла взят из Приложения A.1 к OpenMP Application Program Interface v2.5:

SUBROUTINE A1(N, A, B)
  INTEGER I, N
  REAL B(N), A(N)
!$OMP PARALLEL DO !I is private by default
  DO I=2,N
    B(I) = (A(I) + A(I-1)) / 2.0
  ENDDO
!$OMP END PARALLEL DO
END SUBROUTINE A1

Пожалуйста, обратите внимание:

• -fopenmp подразумевает -frecursive, т.е. все локальные массивы будут выделены в стеке. При переносе существующего кода в OpenMP это может привести к неожиданным результатам, особенно к ошибкам сегментации, если размер стека ограничен.
• На системах, основанных на glibc, приложения с включенным OpenMP не могут быть статически слинкованы из-за ограничений в реализации pthreads. Возможно получение работоспособного решения, если добавить -Wl,--whole-archive -lpthread -Wl,--no-whole-archive в командную строку. Однако это не поддерживается gcc и поэтому не рекомендуется.