Ubuntu下Fortran多线程编程指南

Ubuntu下Fortran多线程编程指南

一 环境准备与编译

1. 安装编译器与运行时
   1. 安装 GNU Fortran 编译器：sudo apt update && sudo apt install gfortran
   2. 可选，安装 OpenMP 运行时与开发包：sudo apt install libomp-dev
2. 编译与运行
   1. 启用 OpenMP：在 gfortran 命令中加入 -fopenmp
   2. 示例：gfortran -fopenmp -O3 -o omp_demo omp_demo.f90 && ./omp_demo
3. 版本与特性
   1. GNU Fortran（gfortran）自版本 4.9 起原生支持 OpenMP，适用于共享内存多线程并行。

二 OpenMP快速上手

1. 常用指令与子句
   1. 并行区域：!$omp parallel [private(…)] [shared(…)]
   2. 工作共享：!$omp do [schedule(…)] [reduction(…)]
   3. 同步：!$omp barrier, !$omp critical, !$omp single
   4. 运行时库：使用 use omp_lib 调用如 omp_get_thread_num()、omp_get_num_threads()、omp_set_num_threads()
2. 最小示例
   1. 说明：并行打印线程号，演示 private/shared 与库函数用法
   2. 代码：

      program omp_hellouse omp_libimplicit noneinteger :: i, nn = 8call omp_set_num_threads(4)!$omp parallel private(i) shared(n)i = omp_get_thread_num()print '("Hello from thread ", i0, " of ", i0)', i, omp_get_num_threads()!$omp end parallel!$omp parallel do schedule(static,1) reduction(+:n)do i = 1, 100 n = n + iend do!$omp end parallel doprint '("Reduction sum = ", i0)', nend program omp_hello

   3. 编译运行：gfortran -fopenmp -O2 -o omp_hello omp_hello.f90 && ./omp_hello
3. 要点
   1. 循环迭代变量通常声明为 private
   2. 存在数据竞争时，使用 reduction、critical 或原子操作
   3. 默认调度为 schedule(static)，可根据负载特性选择 dynamic/guided。

三 线程控制与性能调优

1. 环境变量
   1. 设置线程数：OMP_NUM_THREADS=8
   2. 线程绑定：OMP_PROC_BIND=true（提升缓存局部性）
   3. 调度策略：OMP_SCHEDULE=“dynamic,4”
2. 运行时 API
   1. 在程序内设置线程数：call omp_set_num_threads(N)
   2. 获取线程信息：omp_get_thread_num()、omp_get_num_threads()
3. 调优建议
   1. 线程数尽量与物理核心数匹配；超线程场景下可尝试 2×核心数 做实验
   2. 避免线程间频繁共享写，优先使用 reduction 或局部变量
   3. 合理设置 schedule 与 chunk size，I/O 密集段谨慎并行
   4. 编译时开启优化：-O3 -march=native，并链接数学库（如 -lm）以确保浮点性能。

四 常见问题与排查

1. 未启用并行
   1. 现象：输出始终为 1 个线程
   2. 排查：确认编译命令含 -fopenmp；确认运行环境未强制设 OMP_NUM_THREADS=1
2. 数据竞争与结果不确定
   1. 现象：结果偶发不一致
   2. 排查：将循环迭代变量设为 private；对累加/归约使用 reduction；必要时加 critical/atomic
3. 打印交错混乱
   1. 现象：线程输出混杂
   2. 排查：在打印段加 !$omp critical 或使用缓冲/日志文件
4. 性能不升反降
   1. 排查：线程过多导致上下文切换；调度不合适；存在锁竞争；内存带宽瓶颈
5. 兼容性
   1. 确认 gfortran ≥ 4.9；旧版本不支持 OpenMP 或支持不完整。

五 进阶并行模型

1. 混合编程 OpenMP + MPI
   1. 适用：节点内多线程 + 节点间多进程
   2. 安装 MPI：sudo apt install libopenmpi-dev openmpi-bin
   3. 编译运行：mpif90 -fopenmp -O3 -o hybrid hybrid.f90；mpirun -np 2 --bind-to core ./hybrid
2. 适用场景
   1. OpenMP：共享内存、单机多核、快速并行化
   2. MPI：分布式内存、跨节点扩展、大规模并行