24小时热门版块排行榜    

查看: 1539  |  回复: 7

fanerye2015

新虫 (初入文坛)

[求助] udf程序出错,大家帮忙看下 已有5人参与

大家好,这是我写的一个在z向动量方程上,添加一个力,编译出错,大家帮忙看下是哪里出了错
DEFINE_SOURCE(zforce,c,t,ds,eqn)
{
        real zrc,x,y,z,zfen,frc,T;
        Thread *sec_th;       
        sec_th=THREAD_SUB_THREAD(t,1);
        T=C_T(c,t);

        x=C_T_G(c,t)[0];  //三个方向的温度梯度
        y=C_T_G(c,t)[1];
        z=C_T_G(c,t)[2];

        frc=101300*exp(12.37*(T-3200)/T)*2.5*0.000000001; //计算力的大小
        zfen=sqrt((z*z)/(x*x+y*y+z*z));  //z方向力的分量

        if(0.1<C_VOF(c,t,1)&&C_VOF(c,t,1)<0.8&&T>=3200)
                zrc=frc*zfen;
        else
                zrc=0;
        return zrc;
}
编译之后,提示错误
==============================================================================

Node 0: Process 47276: Received signal SIGSEGV.

==============================================================================

==============================================================================

Node 1: Process 42852: Received signal SIGSEGV.

==============================================================================

==============================================================================

Node 2: Process 36820: Received signal SIGSEGV.

==============================================================================

==============================================================================

Node 3: Process 45180: Received signal SIGSEGV.

==============================================================================
MPI Application rank 0 exited before MPI_Finalize() with status 2
The fl process could not be started.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (文坛精英)

感谢参与,应助指数 +1
本帖仅楼主可见
2楼2016-01-24 03:55:49
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

fanerye2015

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 硫酸亚铜 at 2016-01-24 03:55:49
你不用并行计算试试呢

也不行呜呜呜呜

发自小木虫Android客户端
3楼2016-01-25 22:22:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxh_qddx

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
x,y,z都是0怎么办……分母啊

发自小木虫Android客户端
4楼2016-01-25 22:35:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

summerho

铜虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
看程序没啥问题。你重新启动一下电脑。试试非编译的那个看能不能加载?
5楼2016-01-26 10:44:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微流咨询

新虫 (正式写手)


【答案】应助回帖

这个错误一般是计算的时候才出现,而你现在编译就出现这个错误,问题应该不是出在udf上。
6楼2016-01-26 14:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

死天

专家顾问 (著名写手)

【答案】应助回帖

ds[epn]的值呢?看帮助的源项UDF怎么写的。。。。。即使没有,也要赋值为0
7楼2016-01-27 14:08:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3_1415926

金虫 (小有名气)

可能会有问题的是:
1. 这句:sec_th=THREAD_SUB_THREAD(t,1); 如果加载在第二相上的话,DEFINE_SOURCE(zforce,c,t,ds,eqn)中的t就是第二相的thread
2. zfen=sqrt((z*z)/(x*x+y*y+z*z));这句为0情况
8楼2016-01-27 19:05:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanerye2015 的主题更新
信息提示
请填处理意见