Fortran Coder

标题: 关于利用OPENMP并行的一个问题 [打印本页]

---

作者: 白手002    时间: 2020-4-21 22:08
标题: 关于利用OPENMP并行的一个问题
在并行计算中如果我想要在运行了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
求大佬指点下迷津

---

作者: li913    时间: 2020-4-22 10:30
随机数子程序不能用于并行域，可能会出问题，三个部分得到相同的值。

[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

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

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

---

作者: li913    时间: 2020-4-22 14:35
先生成随机数，再并行。

---

作者: 白手002    时间: 2020-4-22 21:08

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

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

---

作者: li913    时间: 2020-4-22 21:55
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

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

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

---

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)

Powered by Discuz! X3.2