UDF编译出错，求各位路过的帮忙看看

hujilin1229

下面是我的一个UDF程序代码，还望各位大虾指正，这是一个筒形的电场分布，然后求带电电荷的电场力，公式为
F=Q*US/(r*ln(b/a)，不过compiled时出现错误
#include "udf.h"
#define Q 80*1.69e-15         /* particle electric charge      */
#define US 15000        /* z component of ic field */
#define TSTART 0.0   /* field applied at t = tstart   */
#define a 0.0005
#define b 0.25
DEFINE_DPM_BODY_FORCE(particle_body_force,p,i)
{
        real bforce;
                if(P_TIME(p)>TSTART)
                {
                   if(i==0)
                           bforce=Q*US/(P_POS(p)[0]*log(b/a));
                   else
                           bforce=0.0;
                }
                else
                        bforce=0.0;
        /* an acceleration should be returned */
        return (bforce/P_MASS(p));
}

[ 本帖最后由 hujilin1229 于 2011-5-9 22:00 编辑 ]

swimming02

把错误内容贴出来啊