Fortran OpenMP代码：每次迭代仅执行一次任务及嵌套循环并行疑问

嘿，作为Fortran OpenMP新手碰到这个线程利用率的问题太正常了！咱们先搞清楚为啥会出现这种情况，再给你几个实用的解决方案。

问题根源 #

你现在只给内层循环加了并行指令，相当于每次外层i迭代开始时，线程组要等上一轮i对应的内层j循环全跑完，才能重新被调用处理下一个i的任务。本质是外层循环完全串行，内层循环的并行是"单次性"的，线程没法同时处理不同i的迭代，自然会出现等待，浪费了CPU资源。

解决方案 #

根据你的变量var的特性，有两种常用的优化方式：

1. 合并两层循环并行（推荐优先尝试）

如果每个外层i对应的var是独立的，而且内层循环的计算（比如数组访问）没有跨i的依赖，那直接用COLLAPSE(2)把两层循环合并并行化，让线程同时处理不同i和j的迭代，彻底消除外层串行的等待。

示例代码：

!$OMP PARALLEL DEFAULT(NONE) PRIVATE(i,j,var) SHARED(array, N, M)
!$OMP DO COLLAPSE(2)
do i = 1, N  ! 规模较小的外层循环
    var = some_function(i)  ! 每个i的var独立计算
    do j = 1, M  ! 规模最大的内层循环
        ! 你的核心计算逻辑，比如：
        array(i,j) = var * sin(real(j, kind=8))
    end do
end do
!$OMP END DO
!$OMP END PARALLEL

• 关键说明：COLLAPSE(2)会把两层循环展开成一个大的迭代空间，线程均匀分配任务；PRIVATE(i,j,var)确保每个线程有自己的循环变量和var副本，避免数据竞争。
• 编译注意：需要用支持OpenMP 3.0及以上的编译器（现在主流的GCC、Intel Fortran都支持），编译时加OpenMP选项（比如GCC用-fopenmp）。

2. 外层循环并行化+内层串行/嵌套并行

如果var的计算本身比较耗时，或者你不想合并循环，那可以把并行区域提到外层循环，让不同线程同时处理不同的i迭代，这样下一个i的任务不用等上一个i的j循环全跑完。

示例代码：

!$OMP PARALLEL DEFAULT(NONE) PRIVATE(i,j,var) SHARED(array, N, M)
!$OMP DO
do i = 1, N
    var = some_function(i)
    ! 这里内层循环可以选择串行（如果单个i的j循环已经足够大，能占满线程）
    do j = 1, M
        array(i,j) = var * cos(real(j, kind=8))
    end do
    ! 要是还想内层并行，需要开启嵌套并行（编译时加-fopenmp-nested，GCC为例）
    !$OMP DO
    do j = 1, M
        array(i,j) = var * cos(real(j, kind=8))
    end do
    !$OMP END DO
end do
!$OMP END DO
!$OMP END PARALLEL

• 关键说明：外层!$OMP DO会把i的迭代分配给不同线程，线程可以同时处理各自的i和对应的j循环，不会互相等待；嵌套并行适合内层循环规模极大的场景，但会增加线程管理开销，一般优先用第一种方案。

重要注意事项 #

• 数据依赖检查：一定要确保不同i的计算之间没有冲突，比如不要让多个线程同时写入同一个数组位置。
• 变量作用域：var是外层i循环定义的，必须设为PRIVATE，否则会出现多个线程共用同一个var的错误。
• 性能测试：两种方案都可以测试一下，看哪种更适合你的实际代码（比如合并循环的缓存友好性可能更好）。

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