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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 500金币，帮忙找错，UDF并行
161/2返回列表12下一页
查看: 2256  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭