24小时热门版块排行榜

返回列表

_yixiao

铁虫 (小有名气)

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

[交流] 【求助】传热计算的OpenMP并行计算

请大家帮帮忙
在温度场的计算过程中，时间步长上的迭代过程是相互关联，相互影响的。一次迭代需要使用上一次迭代的运算结果，故迭代之间不宜实现并行化。而在一次迭代内部，温度场的计算仅需要前一次时刻的计算结果，各个计算过程之间没有影响，可以进行并行计算。但我进行以下的并行计算后结果同串行相差很大，且计算时间反而增加。
以下即是我的并行程序，用到的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));
}

}
}

}

}
}

回复此楼

» 猜你喜欢

基金啊基金已经有3人回复
科研孤儿太难了已经有15人回复
让我中一个面上吧！已经有10人回复
言之凿凿，明天官网会通知放榜？已经有17人回复
时间戳又变了已经有17人回复
今天系统多次维护，明天很可能放榜！已经有11人回复
时间戳变了，能看出什么问题？已经有23人回复
只有每年这种时候来逛逛小木虫已经有22人回复
放榜前的不淡定已经有12人回复
重要消息，中午系统在维护已经有12人回复

信仰我所必须相信的

1楼 2010-09-17 16:35:16

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

相关版块跳转我要订阅楼主 _yixiao 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 让我中一个面上吧！ +8	大萍1987 2026-08-20	10/500	2026-08-20 21:48 by 大萍1987
[基金申请] 言之凿凿，明天官网会通知放榜？ +15	医学老男孩 2026-08-20	17/850	2026-08-20 20:35 by CaoLG
[基金申请] 时间戳变了，能看出什么问题？ +18	基诺咪客 2026-08-17	23/1150	2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21	yaoyewhu2008 2026-08-20	22/1100	2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3	snowwithsea 2026-08-19	12/600	2026-08-20 15:56 by lfy8008
[基金申请] 重要消息，中午系统在维护 +11	yuleib84 2026-08-18	12/600	2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗？20260819 +14	kkkl_v 2026-08-19	15/750	2026-08-20 08:21 by sunzitan
[基金申请] filecode，4个jtjc了 +11	ziyangfang 2026-08-19	13/650	2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5	瞬息宇宙 2026-08-15	5/250	2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4	yuleib84 2026-08-19	5/250	2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗？ +14	布布和一二 2026-08-19	15/750	2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4	奶牛小黑 2026-08-15	10/500	2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6	wangzilk 2026-08-18	8/400	2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖？ +6	CrisMessi 2026-08-18	6/300	2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5	Shxjjxjkx 2026-08-18	5/250	2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14	archvillain 2026-08-15	26/1300	2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护祝所有人高中 +8	gjjjzhong 2026-08-18	9/450	2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4	且听虎啸 2026-08-17	5/250	2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘金属平磨液，抛光液研发工程师 +3	小天0311 2026-08-14	3/150	2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7	Tide man 2026-08-14	8/400	2026-08-15 01:11 by arzu_hma