关于利用OPENMP并行的一个问题
在并行计算中如果我想要在运行了50步时交换t1和t2的数值,在100步的时候交换t2和t3的数值,150步的时候交换t1和t2的数值,在200步的时候交换t2和t3的数值.依次类推,直到循环完为止。以下是代码program main
implicit none
real(kind=8) :: a,b,c,e0,e1,e2,t1,t2,t3,x,r1,r2,de,p
integer :: i
OPEN(31,FILE='Ea.CIR80',STATUS='unknown')
OPEN(32,FILE='Eb.CIR80',STATUS='unknown')
OPEN(33,FILE='Gc.CIR80',STATUS='unknown')
t1 = 300.d0
t2 = 310.d0
t3 = 320.d0
e0 = 10000
call random_seed ()
!$OMP PARALLEL SECTIONS
!$OMP SECTION
do i = 1,10000
call random_number (x)
r1 = dble(x)
r2 = dble(x) * 2
e1 = t1 * e0
e2 = e1 * r2
de = e2 - e1
if(de.le.0) then
p = 1.d0
else
p = dexp(-de)
end if
if(p.le.r1) then
e2 = e1
else
e1 = e2
write(*,*) i
write(31,*) i,e2
end if
end do
close(31)
!$OMP SECTION
do i = 1,10000
call random_number (x)
r1 = dble(x)
r2 = dble(x) * 2
e1 = t2 * e0
e2 = e1 * r2
de = e2 - e1
if(de.le.0) then
p = 1.d0
else
p = dexp(-de)
end if
if(p.le.r1) then
e2 = e1
else
e1 = e2
write(*,*) ' ',i
write(32,*) i,e2
end if
end do
close(32)
!$OMP SECTION
do i = 1,10000
call random_number (x)
r1 = dble(x)
r2 = dble(x) * 2
e1 = t3 * e0
e2 = e1 * r2
de = e2 - e1
if(de.le.0) then
p = 1.d0
else
p = dexp(-de)
end if
if(p.le.r1) then
e2 = e1
else
e1 = e2
write(*,*) ' ',i
write(33,*) i,e2
end if
end do
close(33)
!$OMP END PARALLEL SECTIONS
end program main
求大佬指点下迷津 随机数子程序不能用于并行域,可能会出问题,三个部分得到相同的值。logical b
b=.false.
do i=1,10000
if(mod(i,50)==0) then
b=.not.b
if(b) then
交换t1,t2
else
交换t2,t3
end if
end if
end do li913 发表于 2020-4-22 10:30
随机数子程序不能用于并行域,可能会出问题,三个部分得到相同的值。logical b
...
谢谢大佬的解答。
那如果我想在并行域产生随机数应该怎么做呢? 先生成随机数,再并行。 li913 发表于 2020-4-22 14:35
先生成随机数,再并行。
谢谢大神。最后在问一下如果有t1,t2,t3,t4一共四个量,也按照50步依次交换的话该怎么写呢? integer n
n=0
do i=1,10000
if(mod(i,50)==0) then
n=mod(n+1,3)
if(n==0) then
交换t1,t2
elseif(n==...)
交换...
end if
end if
end do li913 发表于 2020-4-22 21:55
integer n
n=0
do i=1,10000
哇,大佬牛B,这样就可以想加多少量就加多少了
页:
[1]