24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2203  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

18772958660

新虫 (初入文坛)

[求助] 怎样把用Fortran编写的程序改写为并行程序 已有2人参与

有没懂Fortran并行编程的,目前我面临的一个问题:用Fortran 编写的程序在服务器上跑了几天就不出结果,想把它改成并行的用多核计算,但没学过并行的,所以很头疼,需要一位大神帮帮我。或者是优化一下程序的算法,下面一段程序,用单核计算跑的太久,怎样改进一下算法:
subroutine initial_condition()
  use global
  implicit none
    real :: r1, p1 = 0.5
        do i = 1,n
                call random_number(r1)
                if(r1 < p1)then
                        s(i) = 1.0
                else
                        s(i) = -1.0
                end if
                sum1 = sum1 + s(i)
        end do
        return
end

subroutine mcs()
  use global
  implicit none
    integer :: ii
    real :: r2, r3, sum2 = 0.0
    real :: ws(n)
        do ii = 1,n
                sum2 = 0.0
                call random_number(r2)
                i = n*r2 + 1
                sum1 = sum1 - s(i)
                do j = 1,n
                        if(ad(i,j) == 1)then         
                                sum2 = sum2 + s(j)
                        end if
                end do
                if(sum2<0)then
                        ws(i) = (1-(1-2*q)*s(i)*(-1.0))/2
                else if(sum2>0)then
                        ws(i) = (1-(1-2*q)*s(i)*(1.0))/2
                else
                        ws(i) = 0.5
                end if
                call random_number(r3)
                if(r3 <= ws(i))then
                        s(i) = -s(i)  
                else
                        s(i) = s(i)
                end if
                sum1 = sum1 + s(i)
        end do
        return
end

!**************************Ising dynamics*******************************
subroutine dynamics()
use global
implicit none
integer :: mm, t = 0
real :: sum3 = 0.0, sum4 = 0.0, sum5 = 0.0
real :: m(time)
do mm = 1,n1
q = q + 0.005
call initial_condition()
do while(t < time)
        call mcs()       
        t = t + 1
        m(t) = abs(sum1/n)
        if(t>200000)then
        sum3 = sum3 + m(t)**4
        sum4 = sum4 + m(t)**2
        sum5 = sum5 + m(t)
        end if
end do
MNq(mm) = sum5/(time-200000)
UNq(mm) = 1 - (sum3/(time-200000))/(3*((sum4/(time-200000))**2))
t = 0
sum1 = 0.0
sum3 = 0.0
sum4 = 0.0
sum5 = 0.0
end do
return
end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongwu22

新虫 (初入文坛)

【答案】应助回帖

主函数是什么样的啊?循环内有函数调用,并行的难度比较大啊,而且需要知道一些变量的值,才能确定函数的热点,也就是程序运行时哪最耗时间
有事没事常来玩啊~~
3楼2017-04-25 13:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

yongwu22

新虫 (初入文坛)

【答案】应助回帖


jjdg: 金币+1, 感谢参与 2017-04-25 18:25:18
你打算在什么样的硬件环境跑程序啊,是共享内存还是分布内存啊?
有事没事常来玩啊~~
2楼2017-04-25 13:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

糖爷妖猴

新虫 (初入文坛)

4楼2018-06-07 09:32:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085801 总分275 本科新能源 求调剂 +10 bradoner 2026-04-01 11/550 2026-04-01 22:57 by 来看流星雨10
[考研] 286分调剂 +18 Faune 2026-03-30 20/1000 2026-04-01 21:37 by 495374996
[考研] 085601材料工程找调剂 +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[考研] 337求调剂 +9 《树》 2026-03-29 9/450 2026-04-01 18:05 by 记事本2026
[考研] 调剂 +5 好好读书。 2026-03-28 7/350 2026-04-01 15:32 by 王亮_大连医科大
[考研] 085410人工智能 初试316分 求调剂 +3 残星拂曙 2026-03-31 3/150 2026-04-01 11:09 by 小熊raider
[考研] 一志愿武理材料工程302调剂环化或化工 +15 Doleres 2026-03-31 16/800 2026-04-01 09:49 by lfj11
[考研] 318求调剂 +3 笃行致远. 2026-03-31 3/150 2026-03-31 20:27 by 求调剂zz
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 272求调剂,接受跨专业调剂! +3 闲鱼卢 2026-03-31 3/150 2026-03-31 13:00 by 替代品000
[考研] 22408 359分调剂 +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
[考研] 本科211总分289,08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 抱歉 +4 田洪有 2026-03-30 4/200 2026-03-30 21:26 by mumin1990
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 297求调剂 +17 田洪有 2026-03-26 18/900 2026-03-30 18:32 by nothing投稿中
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
信息提示
请填处理意见