def connected_angles(VW,Z,L,connections, N_nodes, vertex_points=1, cut = 5):
Cube = [[[] for j in range(N_nodes[0])] for l in range(len(connections))]
for k,i in enumerate(connections):
for j in i:
#find the two connections
#we need to + Adivide by the number of sites to get the right
#major particle index
A=j[0]
B=j[1]+N_nodes[0]
Cube[k][A].append(B)
counter = np.zeros((len(Cube),4))
for k in range(len(Cube)):
count_face = 0
count_all = 0
count_face_less = 0
count_all_less = 0
for i in range(len(Cube[k])):
for j in points.unique(Cube[k][i]):
if Cube[k][i].count(j)>=cut:
count_all += 1
if angle.face_check(VW,Z,k,i,j,1,1,L):
count_face += 1
else:
count_all_less += 1
if angle.face_check(VW,Z,k,i,j,1,1,L):
count_face_less += 1
counter[k] = [count_face,count_all,count_face_less,count_all_less]
print 'finished counting'
return counter
def sides_faces_mod(VW,Z,L,connections, N_nodes, vertex_points=1, cut = 5):
A_points = vertex_points[0]*N_nodes[0]
Cube = [[[] for j in range(sum(N_nodes))] for l in range(len(connections))]
for k,i in enumerate(connections):
for j in i:
#find the two connections
#we need to + Adivide by the number of sites to get the right
#major particle index
A=j[0]
B=j[1]+N_nodes[0]
#This provides the position of sides which are connected
Cube[k][A].append(B)
#add the connections to a list
counter = np.zeros((len(Cube),3))
for k in range(len(Cube)):
count_face = 0
count_all = 0
for j in range(len(Cube[k])):
if len(Cube[k][j]) in cut:
for N in points.unique(Cube[k][j]):
count_all += 1
a = angle.face_check(VW,Z,k,j,N,L)
if a+b==2:
count_face += 1
if count_all == 0:
counter[k] = [ 0,0,0]
else:
counter[k] = [count_face,count_all,count_face/float(count_all)]
print 'finished counting'
return counter
def side_faces(VW, Z, connections, N_nodes, vertex_points = 1, cut = 5):
#Cube = np.zeros((len(connections),N_nodes,6))
Cube = [[[0.0 for i in range(6)] for
k in range(N_nodes)] for l in range(len(connections))]
for k,i in enumerate(connections):
for j in i:
#find the two connections
#we need to divide by the number of sites to get the right
#major particle index
A=j[0]/vertex_points
B=j[1]/vertex_points
print vertex_points
side1 = j[0]%vertex_points/(vertex_points/6)
side2 = j[1]%vertex_points/(vertex_points/6)
#This provides the position of sides which are connected
Cube[k][A][side1].append([side2,B])
count_face = 0
count_all = 0
counter = []
for k in range(len(Cube)):
for j in ge(len(Cube[i])):
for side1 in range(6):
if len(Cube[k][j][i]) >= 1:
for side2 in (Cube[k][j][i]):
count_all += 1
count_face += angle.face_check(VW,Z,k,j,side2[1],side1,side2[1],L)
counter.append([count_all,count_face])
print counter
return counter