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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 500金币，帮忙找错，UDF并行
161/2返回列表12下一页
查看: 1729  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

DJ-Fish

铁虫 (初入文坛)

[求助] 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)

}
}
回复此楼

» 收录本帖的淘帖专辑推荐

弹药库

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-01-19 19:23:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DJ-Fish

铁虫 (初入文坛)
以下是小妹修改的其中一个版本

#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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DJ-Fish

铁虫 (初入文坛)
修改后的UDF 在加上这一句：node_to_host_real_1(flux);之后 报错 ：浮点溢出
赞一下 回复此楼
3楼2013-01-19 19:25:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DJ-Fish

铁虫 (初入文坛)
小妹愚钝，望高手帮忙指点，不胜感激！！！
赞一下 回复此楼
4楼2013-01-19 19:26:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kylafree

至尊木虫 (知名作家)

【答案】应助回帖

★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢参与！ 2013-01-24 19:44:44
对于你的这个代码，并行与串行的代码应该是一样的。你不妨直接把在单核中可用的代码不作任何改动拿到多核中试试
赞一下 回复此楼
5楼2013-01-20 19:43:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

GXFGMB
★ ★
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:51
6楼2013-01-20 19:55:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duzhong
★ ★
xiegangmai: 金币-2, 专业版块，请勿灌水 2013-01-24 19:44:57
7楼2013-01-20 20:21:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangdian

木虫 (正式写手)
上面的只是一个建议，需要试试看，同时建议把你的问题适当描述一下。。。。
赞一下 回复此楼
9楼2013-01-21 10:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gearboy78

禁虫 (小有名气)
本帖内容被屏蔽
10楼2013-01-22 11:03:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 DJ-Fish 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 328，0703考生求调剂，一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
[硕博家园] 木虫好像不热闹了，是不是？ +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
关闭
请填处理意见
关闭