24小时热门版块排行榜

返回列表

sghudng

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 137.9
帖子: 85
在线: 23.9小时
虫号: 879159
注册: 2009-10-21

[交流] 【求助】哪位大侠懂并行UDF的麻烦您进来帮个忙，谢谢！已有1人参与

实在是搞不懂写UDF的时候怎么分配计算，自己照着帮助里的例子套了一天，编译都出错，可是我的非稳态问题如果不并行计算，实在是太慢了，求求各位大哥了……
程序如下：
#include "udf.h"
DEFINE_PROFILE(velocity,t,i)
{
Domain *domain=Get_Domain(1);
Thread *thread_s,*mix_thread;
double flux,flux_tot,velocity;
face_t f;
flux_tot=0.0;
mix_thread=Lookup_Thread(domain,5);
thread_s=THREAD_SUB_THREAD(mix_thread,1);
begin_f_loop(f,mix_thread)
{
      flux=F_FLUX(f,thread_s);
      flux_tot+=flux;
}
end_f_loop(f,mix_thread)
      velocity=2*flux_tot/22.319;
begin_f_loop(f,t)
{
      F_PROFILE(f,t,i)=fabs(velocity);
}
end_f_loop(f,t)
}

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有51人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复

» 本主题相关价值贴推荐，对您同样有帮助:

各位电力大侠给个建议啊，谢谢啦已经有6人回复
怎样降低量子点的荧光光谱半高宽，谢谢大家已经有15人回复
GPC问题，谢谢啦已经有10人回复
一个“懂你想法”的音乐软件。。。豆瓣FM客户端。。。【珍藏】【已搜索无重复】已经有23人回复
求助哪位大侠给分析下涂料分层的问题已经有25人回复
请教一下哪位大侠给提供一些制备柴油的方法不是生物柴油已经有15人回复
麻烦哪位高手帮忙翻译一下，后半段实在是搞不懂，我急着做PPT的，谢谢了。已经有1人回复
哪位大侠用过汽巴784光引发剂，进来帮个忙，谢谢！已经有16人回复
曳力函数的的UDF 已经有15人回复
udf 重力方向变化已经有17人回复
【求助】单机多核fluent并行udf 已经有13人回复
【求助】求助啊那位精通keil51的大侠帮个忙啊，指点迷津！！已经有7人回复
【求助】哪位用过自适应时间步长的麻烦进来帮个忙已经有5人回复
【求助】fluent计算中，udf编译通过，初始化出错已经有11人回复

1楼 2010-07-19 17:04:19

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sghudng

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 137.9
帖子: 85
在线: 23.9小时
虫号: 879159
注册: 2009-10-21

我也知道没人愿意完全帮人写程序，还是把自己照着帮助套着写的贴上来，希望大哥们指点一下其中的错误之处，希望有人能回贴……
#include "udf.h"
DEFINE_PROFILE(velocity,t,i)
{
#if !RP_HOST
Domain *domain=Get_Domain(1);
Thread *thread_s,*mix_thread;
face_t f;
#else
real velocity,flux,flux_tot=0.0;
#endif
#if !RP_HOST
mix_thread=Lookup_Thread(domain,5);
thread_s=THREAD_SUB_THREAD(mix_thread,1);

#if RP_NODE
begin_f_loop(f,mix_thread)
if PRINCIPAL_FACE_P(f,mix_thread)
{
      flux=F_FLUX(f,thread_s);
      flux_tot+=flux;
}
end_f_loop(f,mix_thread)
flux_tot=PRF_GRSUM1(flux_tot);
      velocity=2*flux_tot/22.319;
#endif
#endif
#if RP_NODE
begin_f_loop(f,t)
if PRINCIPAL_FACE_P(f,t)
{
      F_PROFILE(f,t,i)=velocity;
}
end_f_loop(f,t)
#endif
}