24小时热门版块排行榜    

查看: 449  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见