Linux Fortran OpenMP：从任务调用的内部子例程访问全局变量是否合法？

问题背景 #

从OpenMP Task调用的内部子例程中访问程序全局变量是否合法有效？使用ifort 2021.7.0 20220726编译无报错，但运行结果随编译选项变化出现异常：预期全局变量g的输出值始终为42，实际却得到随机值，调整编译选项会出现不同的异常结果。实际项目中存在数十个并行阶段仅读取、不修改的全局变量，子例程调用层级较深，逐个传参成本极高。

示例代码 #

program test1
  implicit none
  integer :: i, j, g
  g = 42
  !$OMP PARALLEL DEFAULT(SHARED)
     !$OMP SINGLE
        i = 0
        j = 1
        do while (j < 60)
           i = i + 1
           !$OMP TASK DEFAULT(SHARED) FIRSTPRIVATE(i,j)
              call sub(i,j)
           !$OMP END TASK
           j = j + j
        end do
     !$OMP END SINGLE
  !$OMP END PARALLEL
  stop

contains

  subroutine sub(i,j)
    implicit none
    integer i,j
    !$OMP CRITICAL(unit6)
       write(6,*) i,j,g
    !$OMP END CRITICAL(unit6)
  end subroutine sub

end program test1

编译命令 #

ifort -o test1 test1.f90 -qopenmp -warn all -check all

原因分析 #

从OpenMP标准层面来说，全局变量（包括Fortran内部子例程访问的宿主变量）在DEFAULT(SHARED)的并行区域和任务中属于共享变量，访问本身是合法的。出现异常的核心原因是编译器优化导致的内存可见性问题：

• 全局变量g仅在并行区域外赋值，并行区域内没有修改操作，编译器可能会将其视为常量进行折叠，或者将值缓存到寄存器中；
• OpenMP Task是延迟调度执行的，当任务最终被其他线程执行时，可能无法正确读取到g的全局内存值，而是读取到错误的缓存值或未初始化的内存区域，导致输出随机值。

解决方案 #

针对这类只读全局变量，无需逐个传参，可通过以下方式解决：

• 显式声明变量的共享属性：在!$OMP TASK指令中显式指定SHARED(g)，强制编译器将g视为任务需要访问的共享变量，避免不正确的优化。修改后的任务指令如下：

  !$OMP TASK DEFAULT(SHARED) FIRSTPRIVATE(i,j) SHARED(g)
     call sub(i,j)
  !$OMP END TASK
  

• 禁用特定优化：如果有多个全局变量，可通过编译选项禁用可能导致该问题的优化（如-O0），但这种方法会影响整体性能，不推荐在生产环境使用；
• 使用模块封装全局变量：将只读全局变量放入模块中，在子例程中通过USE语句引用，这种方式能让编译器更清晰地识别变量的作用域和访问方式，减少优化误判。

内容的提问来源于stack exchange，提问作者Peter McGavin