24小时热门版块排行榜    

查看: 451  |  回复: 2

sghudng

铜虫 (小有名气)

[交流] 【求助】哪位大侠懂并行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)            
}
回复此楼

» 猜你喜欢

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

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

sghudng

铜虫 (小有名气)

我也知道没人愿意完全帮人写程序,还是把自己照着帮助套着写的贴上来,希望大哥们指点一下其中的错误之处,希望有人能回贴……
#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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sghudng

铜虫 (小有名气)

不并行算的时候,程序完全没问题,开并行以后,能导进去,但是程序设置的入口条件不起作用……真没人会吗?
3楼2010-07-20 17:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sghudng 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +3 ssyjh 2026-08-27 3/150 2026-08-27 09:59 by 阿拉贡aragon
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +4 a089 2026-08-26 5/250 2026-08-27 09:24 by Howard28
[基金申请] 基金不中,共勉 +7 eulota 2026-08-26 7/350 2026-08-27 09:02 by hsding2026
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见