请看看我的代码（为什么一运行就卡）

ecjtu2

想把块体是否接触用关系矩阵表达出来，块体已经生成del ball range color 0         
def re_head      ；重新定义表头
bp=ball_head
loop while bp # null
   bnext=b_next(bp) then
     if b_id(bp)=bid then
     ii=b_delete(bp)
     end_if
  bp=bnext
end_loop
end
re_head


def afill
cp=contact_head
array var1(3223,3223)
     loop while cp # null
       bp1=c_ball1(cp)
       bp2=c_ball2(cp)
       clp1=b_clump(bp1)
       clp2=b_clump(bp2)
       m=cl_id(clp1)
       n=cl_id(clp2)
       if c_nstrength(cp)=5.6e3
         var1(m,n)=1
       end_if
     end_loop
     if var1(m,n) # 1
        var1(m,n)=0
     end_if
end
def ashow
loop m (1,3223)
  hed=' '
  msg=' '+string(m)
  loop n (1,3223)
    hed=hed+'       '+string(n)
    msg=msg+' '+string(var1(m,n))
  end_loop
  if m=1
    dum=out(hed)
  end_if
  dum=out(msg)
end_loop
end
afill
set log on
ashow
set log off

ecjtu2

已经解决了