求助，自定义积分后处理

curarpikt

我是个新手啊，现在对模拟完的结果进行后处理就编了简单的udf，式子是如图形式的。但是输出不出来，实在是困扰我好久了，万望高手们帮帮忙，谢谢啦！
编的程序如下：
/********************************************/
/*sn.c*/
/********************************************/

#include "udf.h"

DEFINE_EXECUTE_AT_END(execute_at_end)
{
real x[ND_ND];
real sa=0.0,sb=0.0,r,sn=0.0;
real sc=0.0,sd=0.0;
Domain *d;
Thread *t;
cell_t c;
d=Get_Domain(2);
thread_loop_c(t,d)
{

begin_c_loop(c,t)
{
C_CENTROID(x,c,t);
r=sqrt(x[0]*x[0]+x[1]*x[1]);
sc=C_W(c,t)*C_V(c,t)*r;
sd=C_W(c,t)*C_W(c,t)*r;
sa+=sc*C_VOLUME(c,t);
sb+=sd*C_VOLUME(c,t);
}
end_c_loop(c,t)
sn=sa/sb;
}

printf("sn=\n",sn);
}

gearboy78

Message("sn=\n",sn);

curarpikt

非常感谢，但是出现一个另外的问题：
Error: Chip: internal error: invalid builtin -1: pc=335

Error Object: ()

curarpikt

高手们帮帮忙吧，年底毕业了。做不完我就死定了

lml3504

学习中ing，谢谢分享