关于利用OPENMP并行的一个问题

白手002 · 发表于 2020-4-21 22:08:04

在并行计算中如果我想要在运行了50步时交换t1和t2的数值，在100步的时候交换t2和t3的数值，150步的时候交换t1和t2的数值，在200步的时候交换t2和t3的数值.依次类推，直到循环完为止。以下是代码

[Fortran] 纯文本查看 复制代码

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

白手002 · 发表于 2020-4-21 22:08:05

求大佬指点下迷津

li913 · 发表于 2020-4-22 10:30:59

随机数子程序不能用于并行域，可能会出问题，三个部分得到相同的值。

[Fortran] 纯文本查看 复制代码

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

白手002 · 发表于 2020-4-22 12:19:56

li913 发表于 2020-4-22 10:30
随机数子程序不能用于并行域，可能会出问题，三个部分得到相同的值。[mw_shl_code=fortran,true]logical b
...

谢谢大佬的解答。
那如果我想在并行域产生随机数应该怎么做呢？

li913 · 发表于 2020-4-22 14:35:58

先生成随机数，再并行。

白手002 · 发表于 2020-4-22 21:08:19

li913 发表于 2020-4-22 14:35
先生成随机数，再并行。

谢谢大神。最后在问一下如果有t1,t2,t3,t4一共四个量，也按照50步依次交换的话该怎么写呢？

li913 · 发表于 2020-4-22 21:55: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

白手002 · 发表于 2020-4-23 21:10:51

li913 发表于 2020-4-22 21:55
integer n
n=0
do i=1,10000

哇，大佬牛B，这样就可以想加多少量就加多少了

[并行] 关于利用OPENMP并行的一个问题