yade-users team mailing list archive

[Christian Jakob <question207350@xxxxxxxxxxxxxxxxxxxxx>]

Question #207350 on Yade changed:
https://answers.launchpad.net/yade/+question/207350

Christian Jakob posted a new comment:
Sorry, the code will not work, because you only want sphere-sphere contacts for coordination number.
So you can use this one here:

for b in O.bodies:
  c_coord_tmp = 0
  if isinstance(b.shape,Sphere):
    for i in b.intrs():
      if (i.geom.penetrationDepth > 0):
        c_coord_tmp += 1
    c_coord_sum += c_coord_tmp
    c_balls += 1

coord_num = c_coord_sum/c_balls

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.