24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者一苇以航航将赠送您 10 个金币

一苇以航航

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 2
在线: 6.4小时
虫号: 15713088
注册: 2019-07-24
专业: 可再生与替代能源利用中的

[求助] UDF分段函数求助

想写一个二维模型的分段函数，C1和C2关于x，y的分段函数。但是一直编译不过，摆脱大佬们给看看哪里出错了

DEFINE_ADJUST(dp_adjust,d)
{
Thread *t;
cell_t c;
real x[ND_ND];
real por;
t = Lookup_Thread(d,porous_ID);
begin_c_loop(c,t)
{
     C_CENTROID(x,c,t);
if(xc[0]>=0.0&&xc[0]<=0.3&&xc[1]>=0.0&&xc[1]<0.15)
{
por = dp1;
C_UDMI(c,t,num_c1)=44.5/(fa*pow(por,2.0));
C_UDMI(c,t,num_c2)=2.0*0.55/(por*pow(fa,2.0));
C_UDMI(c,t,num_dpi)=dp1;

}
if(xc[0]>=0.0&&xc[0]<=0.3&&xc[1]>=0.15&&xc[1]<=0.3)
{
por = dp2;
C_UDMI(c,t,num_c1)=44.5/(fa*pow(por,2.0));
C_UDMI(c,t,num_c2)=2.0*0.55/(por*pow(fa,2.0));
C_UDMI(c,t,num_dpi)=dp2;

}
if(xc[0]>0.3&&xc[0]<=0.6&&xc[1]>=0.0&&xc[1]<=0.3)
{
por = dp3;
C_UDMI(c,t,num_c1)=44.5/(fa*pow(por,2.0));
C_UDMI(c,t,num_c2)=2.0*0.55/(por*pow(fa,2.0));

C_UDMI(c,t,num_dpi)=dp3;

}
end_c_loop(c,t);
}

回复此楼