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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 500金币，帮忙找错，UDF并行
161/2返回列表12下一页
查看: 2303  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”：确已离世，正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭