在fortran下进行openmp并行计算编程

最新推荐文章于 2025-06-18 18:18:37 发布

转载

最新推荐文章于 2025-06-18 18:18:37 发布 · 1.9k 阅读

6 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/China3S/p/3500132.html

博主分享了在Fortran中使用OpenMP进行并行计算的经验，强调了临界与归约、条件并行以及负载平衡的重要性。通过示例代码解释了如何处理并行中的公共变量问题，以及如何通过条件编译和调度策略来优化并行效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

在fortran下进行openmp并行计算编程

最近写水动力的程序，体系太大，必须用并行才能算的动，无奈只好找了并行编程的资料学习了。我想我没有必要在博客里开一个什么并行编程的教程之类，因为网上到处都是，我就随手记点重要的笔记吧。这里主要是openmp的~

1 临界与归约
在涉及到openmp的并行时，最需要注意的就是被并行的区域中的公共变量，对于需要reduce的变量，尤其要注意，比如这段代码：

program main
implicit none
include 'omp_lib.h'
integer N,M,i
real(kind=8) t
N=20000
t=0.0
!$OMP PARALLEL DO
do i=1,N
t=t+float(i);
M=OMP_get_num_threads()
enddo
write(*, "('t =  ', F20.5, ' running on ', I3, ' threads.')") t,M
pause
stop
end

串行代码可以很容易的得到正确结果：

t = 200010000.00000 running on 1 threads.
不幸的是，如果是并行的话，可能每次都得到一个不同的结