[交流] 邻域网格 udf

#include "udf.h"
#define Ncell 100000 /*网格数;*/
#define NeighC 26 /*邻域网格;*/
cell_t c1,c2,neighbourcell[Ncell][NeighC];
Thread *t; int ii,jj,kk;
DEFINE_INIT(init,domain)
{
for(ii=0;ii for(jj=0;jj neighbourcell[ii][jj]=-1;
thread_loop_c(t,domain)
{
begin_c_loop(c1,t)
{
kk=0; begin_c_loop(c2,t)
{
for(ii=0;ii {
for(jj=0;jj if(C_NODE(c1,t,ii)==C_NODE(c2,t,jj))/*共享Node;*/
{ neighbourcell[c1][kk]=c2;
kk++;
break;
}
if(C_NODE(c1,t,ii)==C_NODE(c2,t,jj)) break;/*退出c1的node循环;*/
}
}
end_c_loop(c2,t)
}
end_c_loop(c1,t)
}
}