24小时热门版块排行榜    

查看: 927  |  回复: 2

luckyhuman1

铜虫 (小有名气)

[求助] 请大伙帮我看看这段UDF,串行没问题,并行就出错了

下面的UDF,我是想把一个出口边界的温度值,增加一个常数,然后给入口边界,这段在串行计算可以计算正确,但是放到并行里面,加边界就出错了
#include "udf.h"
DEFINE_PROFILE(tem_profile,t,i)
{
real T_Total=0.0;
real T_average_indutor_out=0.0;
real current_time=CURRENT_TIME;
real n=0;
real T_WALL;
real T_inductorinlet=0.0;
real T_inlet;
real x[ND_ND];
real m=0.0;

#if !RP_HOST
face_t f;
cell_t c;
Thread *tc;
Thread *t0;
Domain *domain=Get_Domain(1);
#endif

#if !RP_NODE
FILE *fout;
FILE *fout123;
#endif

#if !RP_NODE
fout=fopen("inductor_outlet_aveT.txt","a+";
fout123=fopen("inductor_inlet.txt","a+";
Message("\n  openfile \n";
#endif

#if !RP_HOST
tc=Lookup_Thread(domain,16);
begin_f_loop_int(f,tc)
{
  T_WALL=F_T(f,tc);
  T_Total+=T_WALL;
  n++;
}
end_f_loop_int(f,tc)

# if RP_NODE
T_Total=PRF_GRSUM1(T_Total);
n=PRF_GRSUM1(n);
#endif

#endif

node_to_host_real_2(T_Total,n);

#if !RP_NODE
T_average_indutor_out=T_Total/n;
fprintf(fout," %lg  %lg \n ",current_time,T_average_indutor_out);
fclose(fout);
#endif

#if !RP_HOST
begin_f_loop_int(f,t)
{
   F_PROFILE(f,t,i)=T_average_indutor_out+8;
   T_inductorinlet+=F_T(f,t);
   m++;
}
end_f_loop_int(f,t)

# if RP_NODE
T_Total=PRF_GRSUM1(T_inductorinlet);
m=PRF_GRSUM1(m);
#endif
#endif

node_to_host_real_2(T_inductorinlet,m);

#if !RP_NODE
T_inlet=T_inductorinlet/m;
fprintf(fout123," %lg  %lg \n ",current_time,T_inlet);
fclose(fout123);
#endif
}
回复此楼

» 猜你喜欢

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

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

luckyhuman1

铜虫 (小有名气)

xiegangmai: , 建议分享解决方法,金币奖励! 2012-05-20 22:46:02
已经解决了,谢谢关注
2楼2012-05-19 16:37:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

missyou8326

金虫 (著名写手)

解决办法是什么,给大家参考一下啊
3楼2012-05-20 21:15:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luckyhuman1 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 买卖文章的刷屏了! +3 能丫 2026-06-13 3/150 2026-06-14 09:20 by bobvan
[考研] *O_o^ +3 jk31rb6vyz 2026-06-12 3/150 2026-06-14 09:08 by 8bmv3h4hvh
[硕博家园] T^T +7 m584r2z4j3 2026-06-11 8/400 2026-06-14 08:34 by 8bmv3h4hvh
[公派出国] T^T +6 q1zg6x4lf3 2026-06-12 6/300 2026-06-14 08:28 by 8bmv3h4hvh
[考博] →_→ +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-14 08:28 by 8bmv3h4hvh
[考研] O(∩_∩)O +4 q1zg6x4lf3 2026-06-12 5/250 2026-06-14 08:23 by 8bmv3h4hvh
[考研] ╭(╯ε╰)╮ +5 m584r2z4j3 2026-06-11 5/250 2026-06-14 08:03 by 8bmv3h4hvh
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +3 q1zg6x4lf3 2026-06-12 3/150 2026-06-13 19:07 by 1js5n6upxm
[公派出国] +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 19:05 by 1js5n6upxm
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[考研] o>_<o +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 18:22 by lkmi35jdzt
[考博] T^T +4 jk31rb6vyz 2026-06-12 4/200 2026-06-13 12:09 by 80j16ujcu8
[基金申请] 面上项目没有好文章就没希望了吗? +19 1234567wang 2026-06-07 21/1050 2026-06-12 21:24 by littlezl
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[基金申请] 前几天时间戳更新了 +8 vito刘 2026-06-11 13/650 2026-06-12 14:17 by BlakeReary
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见