24小时热门版块排行榜    

查看: 965  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +38 AI8RaGaPaSCR 2026-08-07 40/2000 2026-08-13 23:45 by 9kfQ6ZJVaWOC
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +10 医学老男孩 2026-08-13 28/1400 2026-08-13 22:44 by 医学老男孩
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
请填处理意见