| 查看: 450 | 回复: 2 | |||
[交流]
【求助】哪位大侠懂并行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) } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位电力大侠给个建议啊，谢谢啦
已经有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人回复
|
我也知道没人愿意完全帮人写程序，还是把自己照着帮助套着写的贴上来，希望大哥们指点一下其中的错误之处，希望有人能回贴…… #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 } |
2楼2010-07-20 10:32:41
3楼2010-07-20 17:06:00
25