通过并行化加速含独立do循环的Fortran代码是否可行？

可行性判断 #

• 你的场景完全可以通过并行计算提升运行速度，这种无跨迭代依赖的循环属于典型的易并行计算场景，是并行优化收益最高的场景之一。如果单次循环内的计算量足够大，并行后的理论加速比接近可用CPU核心数。
• 注意：如果单次循环内的计算量极小（仅包含几次算术运算），多核心并行的线程调度开销可能抵消收益，这种场景优先使用SIMD向量化优化即可。

仅开启Intel Fortran并行选项的适用范围 #

• Intel Fortran提供的自动并行编译选项（ifort/ifx对应参数为-qparallel，需搭配-O2及以上优化等级使用）确实可以在不修改代码的前提下自动识别无依赖循环并做并行化，但存在明确局限：
  • 仅能处理逻辑非常简单的循环，一旦循环内存在函数调用、指针操作、间接数组索引等编译器无法100%确认无依赖的逻辑，就会放弃自动并行
  • 自动并行的任务拆分策略固定，无法根据实际计算负载调整，容易出现负载不均、并行开销超过收益的问题
  • 无法对嵌套循环做针对性的并行层级选择，容易选错并行粒度反而降低性能

更稳妥的实现方案 #

• 不需要编写复杂的并行子程序，仅需添加极少的OpenMP编译指示即可实现可控的循环并行，改造成本极低，改造示例如下：

!$omp parallel do private(i)  ! 声明循环变量i为线程私有，避免访问冲突
do i=1,1000
    ! 原有计算逻辑不需要做任何修改
end do
!$omp end parallel do

• 编译时添加-qopenmp参数即可启用OpenMP并行，运行时会自动根据环境变量OMP_NUM_THREADS设置的线程数分配循环迭代到不同核心执行，比自动并行的稳定性、可控性高得多，是Fortran循环并行的主流实现方案。

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