版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

_yixiao

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1216.5
散金: 15
帖子: 135
在線: 13小時
蟲號: 843338
注冊: 2009-09-08
專業(yè): 鋼鐵冶金

[交流] 【求助】傳熱計算的OpenMP并行計算已有1人參與

請大家?guī)蛶兔?br /> 在溫度場的計算過程中，時間步長上的迭代過程是相互關(guān)聯(lián)，相互影響的。一次迭代需要使用上一次迭代的運算結(jié)果，故迭代之間不宜實現(xiàn)并行化。而在一次迭代內(nèi)部，溫度場的計算僅需要前一次時刻的計算結(jié)果，各個計算過程之間沒有影響，可以進行并行計算。但我進行以下的并行計算后結(jié)果同串行相差很大，且計算時間反而增加，不知道問題出在哪里。
以下即是我的并行程序，用到的OpenMP指令是section指令：
（其中，程序中的h_w_f，Water_Tem_in，Mould_dt，Mould_dy等為全局變量，AA、BB、C_mould，float Density_mould等為局部變量，q1，Q2,q3，Q4為調(diào)用函數(shù)。）
for(k=0;k<100;k++)
{

            for(i=0;i<=xx;i++)
            {
                     for(j=0;j<=yy;j++)
                     {
                              thermal_1[j]=thermal_0[j];
                     }
            }

#pragma omp parallel sections
            {
#pragma omp section
                     {

                              for(i=0;i<=xx;i++)
                              {
                                    for(j=0;j<=yA;j++)
                                    {
                              if (i==xK && j==0)//K
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i+1][j]+2*BB*thermal_1[j+1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_f*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xJ && j==0)//J
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i-1][j]+2*BB*thermal_1[j+1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_f*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }

                              //////////[side]
                              else if ((i>=xK && i<=xJ) && j==0)//KJ
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_f*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xJ && (j>=0 && j<=yA))//JP
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i-1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if ((i>=xK && i<=xJ) && j==yA)//PL
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j-1])
                                             +thermal_1[j]*(1-2*AA-2*BB)+2*Q4(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xK && (j>=0 && j<=yA))//LK
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if ((i>=xK && i<=xJ) && (j<=yA) )// fix
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1.0-2*AA-2*BB));
                              }

                                    }

                     }
            }

pragma omp section
                     {
                              for(i=0;i<=xx;i++)
                              {
                                    for(j=yB;j<=yy;j++)
                                    {
                              if (i==xJ && j==yy)//E
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i-1][j]+2*BB*thermal_1[j-1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xK && j==yy)//D
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i+1][j]+2*BB*thermal_1[j-1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }

                              ////side
                              else if ((i>=xK && i<=xJ) && j==yB)//MN
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB)+2*Q2(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xJ && (j>=yB && j<=yy))//NE
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i-1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if ((i>=xK && i<=xJ) && j==yy)//ED
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j-1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_w_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dy*C_mould*Density_mould));
                              }
                              else if (i==xK && (j>=yB && j<=yy))//DM
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB));
                              }

                              else if((i>=xK && i<=xJ) && (j>=yB))// loose
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1.0-2*AA-2*BB));
                              }

                                    }
                              }
                     }

#pragma omp section
                     {
                              for(i=0;i<=xK;i++)
                              {
                                    for(j=yA;j<=yB;j++)
                                    {
                              if (i==0 && j==yA)//A
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i+1][j]+2*BB*thermal_1[j+1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_n_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==0 && j==yB)//B
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i+1][j]+2*BB*thermal_1[j-1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_n_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==xK && j==yA)//L
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +Q4(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould)+q1(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==xK && j==yB)//M
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +Q2(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould)+q1(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
///side
                              else if ((i>=0 && i<=xK) && j==yA)//AL
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if (i==xK && (j>=yA && j<=yB))//LM
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i-1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB)+2*q1(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if ((i>=0 && i<=xK) && j==yB)//MB
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j-1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if (i==0 && (j>=yA && j<=yB))//BA
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_n_l*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }

                              else if (i<=xK  && (j>=yA && j<=yB))// left
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1.0-2*AA-2*BB));
                              }

                                    }
                              }
                     }

#pragma omp section
                     {
                              for(i=xJ;i<=xx;i++)
                              {
                                    for(j=yA;j<=yB;j++)
                                    {
                              if (i==xJ && j==yA)//P
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +Q4(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould)+q3(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==xx && j==yA)//H
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i-1][j]+2*BB*thermal_1[j+1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_n_r*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==xJ && j==yB)//N
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB)
                                             +Q2(i-xK)*Mould_dt/(Mould_dy*C_mould*Density_mould)+q3(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if (i==xx && j==yB)//G
                              {
                                    thermal_0[j] = float(2*AA*thermal_1[i-1][j]+2*BB*thermal_1[j-1]+thermal_1[j]*(1-2*AA-2*BB)
                                             +2*h_n_r*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
////side
                              else if ((i>=xJ && i<=xx) && j==yA)//PH
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j+1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if (i==xx && (j>=yA && j<=yB))//HG
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i-1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB)+2*h_n_r*(Water_Tem_in - thermal_1[j])*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }
                              else if ((i>=xJ && i<=xx) && j==yB)//GN
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+2*BB*(thermal_1[j-1])+thermal_1[j]*(1-2*AA-2*BB));
                              }
                              else if (i==xJ && (j>=yA && j<=yB))//NP
                              {
                                    thermal_0[j] = float(2*AA*(thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])
                                             +thermal_1[j]*(1-2*AA-2*BB)+2*q3(j)*Mould_dt/(Mould_dx*C_mould*Density_mould));
                              }

                              else if ( i>=xJ && (j>=yA && j<=yB))// right
                              {
                                    thermal_0[j] = float(AA*(thermal_1[i-1][j]+thermal_1[i+1][j])+BB*(thermal_1[j-1]+thermal_1[j+1])+thermal_1[j]*(1.0-2*AA-2*BB));
                              }


                                    }
                              }

                     }

            }
}

本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=2410411

回復(fù)此樓

» 猜你喜歡

326求調(diào)劑已經(jīng)有8人回復(fù)
333求調(diào)劑已經(jīng)有7人回復(fù)
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)學(xué)碩調(diào)劑，一志愿985 已經(jīng)有7人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有21人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有15人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有14人回復(fù)
326求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有13人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

VASP并行計算出錯已經(jīng)有9人回復(fù)
遠程登錄linux系統(tǒng)如何實現(xiàn)fluent并行計算？已經(jīng)有6人回復(fù)
并行計算如何實現(xiàn)k點平均分配？已經(jīng)有5人回復(fù)
gaussian 03Revison-E01版在Window 下不能并行計算啊! 已經(jīng)有9人回復(fù)
intel visual fortran 多核并行計算已經(jīng)有14人回復(fù)
【求助】單機多核UDF怎么并行計算已經(jīng)有15人回復(fù)
中科院深圳先進技術(shù)研究院2011年校園招聘開始了，敬請關(guān)注,待遇已補上已經(jīng)有79人回復(fù)
【求助】并行計算問題已經(jīng)有13人回復(fù)
【求助】castep并行計算問題已經(jīng)有15人回復(fù)

信仰我所必須相信的

1樓 2010-09-17 16:41:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

frankliu

禁蟲 (正式寫手)

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖交流
mze04532(金幣+2):鼓勵一下，歡迎交流~ 2010-09-28 19:38:41

本帖內(nèi)容被屏蔽

2樓2010-09-28 13:45:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 _yixiao 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿南大，0703化學(xué)，分數(shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 306求調(diào)劑 +4	chuanzhu川燭 2026-03-18	4/200	2026-03-21 08:25 by laoshidan
[考研] 332求調(diào)劑 +4	ydfyh 2026-03-17	4/200	2026-03-21 02:20 by JourneyLucky
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-17	8/400	2026-03-21 02:12 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4	葵梓衛(wèi)隊 2026-03-18	6/300	2026-03-20 23:02 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 材料學(xué)求調(diào)劑 +4	Stella_Yao 2026-03-20	4/200	2026-03-20 20:28 by ms629
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3	123阿標 2026-03-17	3/150	2026-03-20 15:58 by 飛行琦
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 085601材料工程專碩求調(diào)劑 +10	慕寒mio 2026-03-16	10/500	2026-03-19 15:26 by 丁丁*
[考研] 一志愿福大288有機化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

_yixiao

[交流] 【求助】傳熱計算的OpenMP并行計算 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

frankliu

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[交流] 【求助】傳熱計算的OpenMP并行計算已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: