24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2205  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程专硕324分,一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 318求调剂,计算材料方向 +7 吸喵有害笙命 2026-04-01 8/400 2026-04-01 23:54 by 义文wang
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
[考研] 261求调剂 +3 明仔· 2026-04-01 3/150 2026-04-01 20:52 by cq2548
[考研] 0817化工学硕调剂 +11 努力上岸中! 2026-03-31 11/550 2026-04-01 20:30 by 赖春艳
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 调剂 +3 好好读书。 2026-04-01 3/150 2026-04-01 17:06 by zhouyuwinner
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +8 YYYYYNNNNN 2026-03-26 10/500 2026-04-01 14:39 by hexingyi
[论文投稿] chinese chemical letters英文版投稿求助 120+4 Yishengeryi 2026-03-30 5/250 2026-04-01 14:11 by 陆小果画大饼
[考研] 0805求调剂 +6 是水分 2026-03-31 6/300 2026-04-01 09:05 by oooqiao
[考研] 353求调剂 +3 江上枫_26 2026-03-28 3/150 2026-03-31 15:53 by jp9609
[考研] 085600,专业课化工原理,320分求调剂 +6 大馋小子 2026-03-29 6/300 2026-03-31 10:03 by 氯化亚硝酰
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 求调剂 +10 家佳佳佳佳佳 2026-03-29 10/500 2026-03-30 18:34 by 544594351
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 070300化学求调剂 +4 起个名咋这么难 2026-03-27 4/200 2026-03-27 21:39 by 83503孙老师
[考研] 一志愿上海理工能源动力(085800)310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
信息提示
请填处理意见