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

登录注册

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

24小时热门版块排行榜

返回列表

苦逼锅炉工

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 0.3
帖子: 41
在线: 10.2小时
虫号: 5883918

[交流] 关于输出捕获的颗粒物的坐标问题

最近在做颗粒物的捕获。即当颗粒物碰到壁面后判定颗粒物速度，然后进行捕获或者反射，并将捕获的颗粒物坐标输出到文件，UDF如下
#include "udf.h"
#include "dpm.h"
FILE *fp;
DEFINE_DPM_BC(JIEZHA,p,t,f,f_normal,dim)
{
fp=fopen("d:\\out.txt","a+"

;
if(NV_MAG(P_VEL(p))<8)
{
   fprintf(fp," %g %g %g\n",P_POS(p)[0],P_POS(p)[1],P_POS(p)[2]);
   return PATH_ABORT;
   fclose(fp);
}
else
   return PATH_ACTIVE;
}

然而在进行非稳态计算后，计算总是会卡在Advancing DPM injections ....，out.txt文件会被创建，但是里面并没有数据，请问是什么原因造成的，我的UDF有什么问题吗

回复此楼

» 猜你喜欢

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

三、关于“天然动力能量场核”构造推想已经有0人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2017-04-13 10:36:23

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

匿名

应助: (幼儿园)
在线:
虫号: 0

★
苦逼锅炉工(金币+1): 谢谢参与

本帖仅楼主可见

赞一下

2楼2017-04-13 11:32:25

已阅申请仿真EPI 回复此楼编辑查看我的主页

苦逼锅炉工

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 0.3
帖子: 41
在线: 10.2小时
虫号: 5883918

自挽

回复此楼

4楼2017-04-14 13:33:22

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

zhengjimi

铁虫 (小有名气)

应助: 2 (幼儿园)
金币: 3.2
帖子: 291
在线: 231.5小时
虫号: 593980

★
苦逼锅炉工(金币+1): 谢谢参与

/* Message(" This FACE thread no. %d..\n", THREAD_ID(t));*/

赞一下

回复此楼

5楼2017-04-27 14:16:36

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

db0320

铁虫 (初入文坛)

应助: 4 (幼儿园)
金币: 47.7
帖子: 42
在线: 14.3小时
虫号: 3608837

★
苦逼锅炉工(金币+1): 谢谢参与

要我说吧，这个问题应该这么看:
1. 能生成out.txt文件，那么这个UDF能够正常调用并至少能够运行至 fp=xxxxx 这一行；
2.能运行fp=xxxxx这一行照理来说没理由不能继续运行if这一行；
3.然而运行了if这一行之后没有记录结果，说明没有任何一次if条件判断是True ；
4.于是乎得出结论，你的判断条件有问题，没有任何一句 NV_MAG(P_VEL(p))<8；
5.怎么验证？把NV_MAG(P_VEL(p))<8中的这个8设定成88888888888888
6.如果这样修改后out.txt中依旧没有数据，那么，我上面说的都是废话(′;ω;｀)
7.再次ps：fclos(fp）这一句永远不会被执行；只要运行到Return就表示本次对该程序的调用结束，并返回Return中的值，Return之后的语句不会被执行；
8.在另一个帖子里回复你关于输出好多点的问题。
9.对8的补充：首先还是确认一下你的case是只在某一个时间点射入制定数量的颗粒，有可能你的case一直在入射颗粒所以才会一直trap。

赞一下

回复此楼

7楼2017-04-27 16:41:04

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

WILLTONG94

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 4119.2
帖子: 476
在线: 44小时
虫号: 5614882

★
苦逼锅炉工(金币+1): 谢谢参与

引用回帖:

7楼: Originally posted by db0320 at 2017-04-27 16:41:04
要我说吧，这个问题应该这么看:
1. 能生成out.txt文件，那么这个UDF能够正常调用并至少能够运行至 fp=xxxxx 这一行；
2.能运行fp=xxxxx这一行照理来说没理由不能继续运行if这一行；
3.然而运行了if这一行之后没 ...

您好！
能否请问一下关于dpm中udf编写的问题？
非常感谢！

赞一下

回复此楼

8楼2017-06-08 09:45:26

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