24小时热门版块排行榜    

查看: 1927  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:48 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 09:52 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
信息提示
请填处理意见