24小时热门版块排行榜

返回列表

iZeeyo

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 29
在线: 21.9小时
虫号: 3929670
注册: 2015-06-17
专业: 流体力学

[求助] 在fortran程序加上openmp并行后耗时比不加并行更慢

在此段程序中，串行计算用时15.5s，运用openmp用两个线程计算，用时17.3s，运用openmp用4个线程计算，用时55s！！！！运行的结果是一样的，但是这个耗时也太不对劲了，希望有经验的大佬指点一下！！！感谢！
程序略长，但是结构简单，用了三个部分，在计算函数部分使用了并行。

module com！全局变量
implicit none
integer a,ia
parameter (a=10000,ia=3)
integer summ(a,a),sumn(a,a),su(ia,3000,3000)
end

program xunhuan！主程序
use com
implicit none
include 'omp_lib.h'

integer i,j,ic
real(kind=8) time_begin,time_end
i=a
j=a

call omp_set_num_threads(2)
call cpu_time(time_begin)
call loop(i,j)
call cpu_time(time_end)

open(10,file="result.txt"

write(10,*) time_end - time_begin
do i=1,500
do j=1,500
write(10,*) i,j,summ(i,j), sumn(i,j),summ(i,j)+sumn(i,j)
do ic=1,ia
   write(10,*) su(ic,i,j)
end do
end do
end do
close(10)
end

subroutine loop(ii,jj)!计算函数
use com
implicit none
integer i1,j1,ii,jj,aa1,bb1,aa2,bb2,cc1,cc2,dd1,dd2,ib

!$omp parallel private(aa1,bb1,aa2,bb2,cc1,cc2,dd1,dd2) shared(summ,sumn,su,ii,jj)
!$omp do

do i1=1,ii
do j1=1,jj
   if(abs(mod(i1,2))==1)then
   aa1=3*i1-j1
   bb1=i1
   else if(abs(mod(i1,2))==0) then
      aa1=i1
      bb1=3*i1-j1
   end if
   cc1=min(aa1,bb1)
   if(abs(i1-25)>=20)then
      aa2=abs(i1-20)*2-j1
      bb2=abs(j1-50)+i1
   else if(abs(i1-25)<20)then
      aa2=abs(i1-10)*2-j1
      bb2=abs(j1-70)+i1
   end if
   cc2=floor(sqrt(1.0*aa2*bb2/abs(aa2-bb2)))

   dd1=cc1+i1
   dd2=cc2+j1


summ(i1,j1)=floor(1.0*sqrt(abs(0.1*(dd1*(j1+1)-i1*dd2)*min(cc1,cc2))/22))+floor(1.0*sqrt(abs(0.01*cc1*cc2)))
sumn(i1,j1)=floor(1.0*sqrt(abs(0.1*((i1+1)*dd2-dd1*j1)*max(cc1,cc2))/22))+floor(1.0*sqrt(abs(0.01*cc1*cc2)))
if(i1<=3000.and.j1<=3000)then
do ib=1,ia
   su(ib,i1,j1)=floor(sqrt(1.0*abs(summ(i1,j1)*cc1)))+floor(sqrt(1.0*abs(sumn(i1,j1)*cc2)))+ib
end do
end if

end do
end do

!$omp end do
!$omp end parallel

end

回复此楼