如何解决OPENMP F90 / 95嵌套DO循环-与串行实现相比有待改进的问题
|| 我已经进行了一些搜索，但是找不到与我的问题有关的任何内容（对不起，如果我的问题很冗长！）。无论如何，正如标题所指出的那样，我很难对我的代码的串行实现进行任何改进。我需要并行化的代码段如下（这是带有OpenMP的Fortran90）：do n=1,lm
do m=1,jm
do l=1,im
sum_u = 0
sum_v = 0
sum_t = 0
do k=1,lm
!$omp parallel do reduction (+:sum_u,sum_v,sum_t)
do j=1,jm
do i=1,im
exp_smoother=exp(-(abs(i-l)/hzscl)-(abs(j-m)/hzscl)-(abs(k-n)/vscl))
sum_u = sum_u + u_p(i,j,k) * exp_smoother
sum_v = sum_v + v_p(i,k) * exp_smoother
sum_t = sum_t + t_p(i,k) * exp_smoother
sum_u_pert(l,m,n) = sum_u
sum_v_pert(l,n) = sum_v
sum_t_pert(l,n) = sum_t
end do
end do
end do
end do
end do
end do
解决方法
放置omp编译指示的明显位置是在外部循环中。
对于每个（l，m，n），您都在计算扰动变量和指数平滑器之间的卷积。每个（l，m，n）计算完全独立于其他计算，因此您可以将其放在最外面的循环上。因此，例如最简单的事情
!$omp parallel do private(n,m,l,i,j,k,exp_smoother) shared(sum_u_pert,sum_v_pert,sum_t_pert,u_p,v_p,t_p),default(none)
do n=1,lm
do m=1,jm
do l=1,im
do k=1,lm
do j=1,jm
do i=1,im
exp_smoother=exp(-(abs(i-l)/hzscl)-(abs(j-m)/hzscl)-(abs(k-n)/vscl))
sum_u_pert(l,n) = sum_u_pert(l,n) + u_p(i,k) * exp_smoother
sum_v_pert(l,n) = sum_v_pert(l,n) + v_p(i,k) * exp_smoother
sum_t_pert(l,n) = sum_t_pert(l,n) + t_p(i,k) * exp_smoother
end do
end do
end do
end do
end do
end do
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。