| 查看: 2141 | 回复: 15 | |||
[求助]
500金币，帮忙找错，UDF并行
|
|
单机下能正常使用，单机多核并行报错： 单机UDF如下 /*UDF功能：从出口取出颗粒质量流量给到入口*/ #include "udf.h" #include "mem.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include "metric.h" #define UDF_FILENAME "hxout_s_flux.dat"/*每次开始新的算例要修改名字否则会接着原来的文件写入数据*/ real last_ts1=-1; DEFINE_PROFILE(inlet_1,t,index) { Domain *domain; Thread *t_s,*t_mix; face_t f; double flux,flux_tot,veloc; real ro_s; real A[ND_ND]; real area; real current_time; current_time=RP_Get_Real("flow-time"  ;/*取得当前的计算时间*/if(last_ts1!=current_time) { last_ts1=current_time; domain=Get_Domain(1);/*1：混合相，2：主相，3：第二相*/ t_mix=Lookup_Thread(domain,5);/*5：出口边界ID*/ t_s=THREAD_SUB_THREAD(t_mix,1);/*0：主相，1：第二相，以此类推*/ begin_f_loop(f,t_mix) { flux=F_FLUX(f,t_s); flux_tot+=flux; } end_f_loop(f,t_mix) begin_f_loop(f,t) { ro_s=F_R(f,t); F_AREA(A,f,t); area+=NV_MAG(A); } end_f_loop(f,t) veloc=flux_tot/ro_s/area/0.3/2.0; begin_f_loop(f,t) { F_PROFILE(f,t,index)=veloc; } end_f_loop(f,t) } } |
回复此楼» 收录本帖的淘帖专辑推荐
 弹药库 |
» 猜你喜欢
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有97人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Fluent中编写UDF之后无法初始化。。。
已经有5人回复
- 【活动】关于papershow活动最终奖项颁发的投票
已经有4人回复
|
以下是小妹修改的其中一个版本 #include "udf.h" #include "mem.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include "metric.h" #define UDF_FILENAME "hxout_s_flux4.0.dat" real last_ts1=-1; real veloc; DEFINE_PROFILE(outlet_pressure,t,index) { real flux; #if !RP_HOST face_t f; Thread *t_s; #endif #if !RP_HOST t_s=THREAD_SUB_THREAD(t,1);/*0：主相，1：第二相，以此类推*/ begin_f_loop(f,t) { F_PROFILE(f,t,index)=-4000.0; flux+=F_FLUX(f,t_s); } end_f_loop(f,t) #endif #if RP_NODE flux=PRF_GRSUM1(flux); #endif #if RP_NODE flux=PRF_GRSUM1(flux); #endif node_to_host_real_1(flux); veloc=flux/2678.0/0.004145/0.3/2.0; } DEFINE_PROFILE(inlet,thread_in,index) { face_t inlet_face; #if !RP_HOST real time=RP_Get_Real("flow-time" ; begin_f_loop(inlet_face,thread_in) { F_PROFILE(inlet_face,thread_in,index)=veloc; } end_f_loop(inlet_face,thread_in) #endif /*printf("%f\n",inlet_yi); */ } |
2楼2013-01-19 19:24:47
3楼2013-01-19 19:25:55
4楼2013-01-19 19:26:50
5楼2013-01-20 19:43:32
★ ★
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:51
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:51
 |
6楼2013-01-20 19:55:49
★ ★
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:57
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:57
|
|
7楼2013-01-20 20:21:43
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+3, 谢谢建议！ 2013-01-24 19:45:19
感谢参与，应助指数 +1
xiegangmai: 金币+3, 谢谢建议！ 2013-01-24 19:45:19
|
首先，按照kylafree的说法试试看，少计算几步对比结果看看。 如果二者结果差别较大，建议做如下改动，然后再对比并行和串行计算结果： 在每个面循环中再加上主面判断。如（其它地方改动类似）： begin_f_loop(inlet_face,thread_in) { if PRINCIPAL_FACE_P(f,tf) F_PROFILE(inlet_face,thread_in,index)=veloc; } end_f_loop(inlet_face,thread_in) |
8楼2013-01-21 10:02:25
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
9楼2013-01-21 10:04:06
|
本帖内容被屏蔽 |
10楼2013-01-22 11:03:07
