版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huanghuang36

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 138.6
帖子: 92
在线: 43小时
虫号: 2050668
注册: 2012-10-09
性别: MM
专业: 机械工程

[求助] UDF 改为并行计算

我的UDF 在单机下能正常使用，但是并行不是需要改一下语句嘛，本人愚钝，照着教程改了改但是一直不对，希望高人能帮忙指点一下，不胜感激啊
如下是我的UDF：主要使用的就是CG_MOTION

#include
#include "udf.h"
#include "math.h"
real r_movefluid;
DEFINE_CG_MOTION(moveface, dt, cg_vel, cg_omega, time, dtime)
{
      Thread *t;
face_t f;
      real x0,y0,r;
real v0;
      real v;
real omega=M_PI;
      real theta=omega*time;
real x[ND_ND];
      t = DT_THREAD(dt);
begin_f_loop(f,t)
      {
            F_CENTROID(x,f,t);
            x0=x[0];
            y0=x[1];
            r=sqrt(pow(x0,2)+pow(y0,2));
      }
end_f_loop(f,t)
      v0=omega*r;
      r_movefluid=r;
if (time>=0 && time<=M_PI/(180*omega))
         v=0;
else if (time<=8*M_PI/(180*omega))
         v=12.96*(theta-(M_PI/180))*omega/(7*pow(M_PI,2));
else if (time<=10*M_PI/(180*omega))
         v=0.072*omega/M_PI;
else if (time<=29*M_PI/(180*omega))
         v=(0.072*omega/M_PI)-12.96*(theta-M_PI/18)*omega/(19*M_PI*M_PI);
else if (time<=31*M_PI/(180*omega))
         v=0;

      cg_vel[0]=v*cos(theta)-v0*sin(theta);
      cg_vel[1]=v*sin(theta)+v0*cos(theta);
      cg_vel[2]=0;
      cg_omega[0]=0;
      cg_omega[1]=0;
      cg_omega[2]=omega;
}
DEFINE_CG_MOTION(rotate, dt, cg_vel, cg_omega, time, dtime)
{
      real omega=M_PI;
      cg_omega[2]=omega;
}
DEFINE_CG_MOTION(movefluid, dt, cg_vel, cg_omega, time, dtime)
{
real v0;
      real v;
real omega=M_PI;
      real theta=omega*time;
      v0=omega*r_movefluid;
if (time>=0 && time<=M_PI/(180*omega))
         v=0;
else if (time<=8*M_PI/(180*omega))
         v=12.96*(theta-(M_PI/180))*omega/(7*pow(M_PI,2));
else if (time<=10*M_PI/(180*omega))
         v=0.072*omega/M_PI;
else if (time<=29*M_PI/(180*omega))
         v=(0.072*omega/M_PI)-12.96*(theta-M_PI/18)*omega/(19*M_PI*M_PI);
else if (time<=31*M_PI/(180*omega))
         v=0;

      cg_vel[0]=v*cos(theta)-v0*sin(theta);
      cg_vel[1]=v*sin(theta)+v0*cos(theta);
      cg_vel[2]=0;
      cg_omega[0]=0;
      cg_omega[1]=0;
      cg_omega[2]=omega;
}

刚注册的，目前只有两个金币，呵呵，不好意思啊~~

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

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

1楼 2012-10-09 10:43:48

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

huanghuang36

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 138.6
帖子: 92
在线: 43小时
虫号: 2050668
注册: 2012-10-09
性别: MM
专业: 机械工程

呵呵，目前我也不是很会啊，共同学习

3楼2012-10-10 10:02:36

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

查看全部 10 个回答

Mr.Silence

新虫 (初入文坛)

应助: 4 (幼儿园)
金币: 105.3
帖子: 34
在线: 5.3小时
虫号: 2023898
注册: 2012-09-23
专业: 流体力学

★
感谢参与，应助指数 +1
xiegangmai: 金币-1, 应助指数-1, 谢谢参与，非应助请勿选择“应助回帖” 2012-10-10 21:19:51

同求并行计算的方法。我也在学习UDF，希望楼主不吝赐教。

2楼2012-10-09 11:17:09

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

sunbird04

银虫 (小有名气)

应助: 36 (小学生)
金币: 710.2
红花: 4
帖子: 143
在线: 98.6小时
虫号: 881246
注册: 2009-10-23
性别: GG
专业: 钢铁冶金

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢应助 2012-10-10 21:20:01
huanghuang36: 金币+2, ★★★很有帮助, 谢谢你的帮助哦~~ 2012-10-16 09:34:24

在每段程序最前头加上 #if !RP_HOST 最末尾加上#endif 即可

赞一下(1人)

» 本帖已获得的红花（最新10朵）

huanghuang36

4楼2012-10-10 18:14:27

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

huanghuang36

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 138.6
帖子: 92
在线: 43小时
虫号: 2050668
注册: 2012-10-09
性别: MM
专业: 机械工程

引用回帖:

4楼: Originally posted by sunbird04 at 2012-10-10 18:14:27
在每段程序最前头加上 #if !RP_HOST 最末尾加上#endif 即可

太感谢你的回答了，我已经照你说的改了程序，还有我想问一下，那个计算节点（也就是处理器）数量的设置是不是和我的子域的数量有关啊，我看书上说子域的数量应该是计算节点的整数倍，是要这么设是吧~~

5楼2012-10-10 21:40:46

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

查看全部 10 个回答