版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

_yixiao

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 1216.5
散金: 15
帖子: 135
在线: 13小时
虫号: 843338
注册: 2009-09-08
专业: 钢铁冶金

[交流] 【求助】传热计算的OpenMP并行计算已有1人参与

请大家帮帮忙
在温度场的计算过程中，时间步长上的迭代过程是相互关联，相互影响的。一次迭代需要使用上一次迭代的运算结果，故迭代之间不宜实现并行化。而在一次迭代内部，温度场的计算仅需要前一次时刻的计算结果，各个计算过程之间没有影响，可以进行并行计算。但我进行以下的并行计算后结果同串行相差很大，且计算时间反而增加，不知道问题出在哪里。
以下即是我的并行程序，用到的OpenMP指令是section指令：
（其中，程序中的h_w_f，Water_Tem_in，Mould_dt，Mould_dy等为全局变量，AA、BB、C_mould，float Density_mould等为局部变量，q1，Q2,q3，Q4为调用函数。）
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://muchong.com/bbs/viewthread.php?tid=2410411

回复此楼

» 猜你喜欢

售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有5人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有3人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有4人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

VASP并行计算出错已经有9人回复
远程登录linux系统如何实现fluent并行计算？已经有6人回复
并行计算如何实现k点平均分配？已经有5人回复
gaussian 03Revison-E01版在Window 下不能并行计算啊! 已经有9人回复
intel visual fortran 多核并行计算已经有14人回复
【求助】单机多核UDF怎么并行计算已经有15人回复
中科院深圳先进技术研究院2011年校园招聘开始了，敬请关注,待遇已补上已经有79人回复
【求助】并行计算问题已经有13人回复
【求助】castep并行计算问题已经有15人回复

信仰我所必须相信的

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

frankliu

禁虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mze04532(金币+2):鼓励一下，欢迎交流~ 2010-09-28 19:38:41

本帖内容被屏蔽

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