| 查看: 3773 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
huanghuang36铜虫 (小有名气)
|
[求助]
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; } 刚注册的，目前只有两个金币，呵呵，不好意思啊~~ |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
b口会评
已经有6人回复
-
信息学部会评时间
已经有3人回复
-
祈祷青基必中
已经有9人回复
-
在职考研
已经有7人回复
-
MSER送审了还被拒稿
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP并行计算出错
已经有9人回复
- 100个金币：udf改成并行的，求高人
已经有22人回复
- 曳力模型udf要并行计算，该怎么修改呢
已经有5人回复
- 并行计算如何实现k点平均分配？
已经有5人回复
- 【求助】用vof模型计算出的界面，有木有udf语句可以找得到
已经有6人回复
- 【求助】单机多核fluent并行udf
已经有13人回复
- 【求助】单机多核UDF怎么并行计算
已经有15人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
- 【求助】linux版fluent使用udf并行问题请教
已经有4人回复
- 【求助】castep并行计算问题
已经有15人回复
9楼2013-07-13 18:54:30
2楼2012-10-09 11:17:09
huanghuang36
铜虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 138.6
- 帖子: 92
- 在线: 43小时
- 虫号: 2050668
- 注册: 2012-10-09
- 性别: MM
- 专业: 机械工程
3楼2012-10-10 10:02:36
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢应助 2012-10-10 21:20:01
huanghuang36: 金币+2, ★★★很有帮助, 谢谢你的帮助哦~~ 2012-10-16 09:34:24
感谢参与，应助指数 +1
xiegangmai: 金币+1, 谢谢应助 2012-10-10 21:20:01
huanghuang36: 金币+2, ★★★很有帮助, 谢谢你的帮助哦~~ 2012-10-16 09:34:24
| 在每段程序 最前头加上 #if !RP_HOST 最末尾加上#endif 即可 |
» 本帖已获得的红花（最新10朵）
huanghuang364楼2012-10-10 18:14:27
50