源项DEFINE_SOURCE怎样加入速度和面积的内积-想编写MHD

ajene

我这段时间想做磁场的UDF，源项中矢量的内积，即在源项中加入速度和表面积矢量的点积值。我将UDF写好了，能通过编译，但是一旦加载就出现读取错误（ACCESS VIOLATION），不知是什么原因，我现在被卡了很长时间了，麻烦您帮忙看看，希望得到你的指点！谢谢了！ UDF如下：
DEFINE_SOURCE(C_a, cell, t, dS, eqn)
{
real source;
real rho=6000;
int n;
face_t f;
Thread *tf;
real V[ND_ND],A[ND_ND];

     c_face_loop(cell,t,n)

     {
                 f=C_FACE(cell,t,n);
                 tf=C_FACE_THREAD(cell,t,n);           
                 NV_DS(V, =, C_U(cell,t), C_V(cell,t),C_W(cell,t), *,rho);     
               
                 F_AREA(A,f,tf);
                 source=NV_DOT(V,A);
           }
    dS[eqn]=0;
   return source;

}

yanbowen89

DEFINE_SOURCE 这个宏本身是基于cells，不需要对cell循环，你可以试下去掉

transient

朋友，我在做这方面的，不知道你进展如何？

ajene

原帖由 yanbowen89 于 2013-10-24 15:15 发表

                               
登录/注册后可看大图

DEFINE_SOURCE 这个宏本身是基于cells，不需要对cell循环，你可以试下去掉

谢谢您的帮助，当我去掉循环后仍然是不能读取！我是想在cell内的face上循环，不知有这样的宏没有