[python]kaj je hitreje

resultFull = [ []*10 for i in range(len(db['cohesive'].keys())*3) ]
#data = [ [0] * 10 ] * len(db['cohesive'].keys())*3
num = 0
for el in db['cohesive'].keys():
  for es in db['elementSet']:
    check = set([i[0] for i in db['elementSet'][es]])
    if el in check:
      final = es
      break
  #for el in temp2:
  #print el
  pt1 = getCoords(db['cohesive'][el]['nodes'][0])
  pt2 = getCoords(db['cohesive'][el]['nodes'][1])
  pt3 = getCoords(db['cohesive'][el]['nodes'][2])
  a = [(pt2[0] - pt1[0]), (pt2[1] - pt1[1]), (pt2[2] - pt1[2])]
  b = [(pt3[0] - pt1[0]), (pt3[1] - pt1[1]), (pt3[2] - pt1[2])]
  normal = cross2(a,b)
  n1 = getFirstDirection(pt1, pt2, pt3)
  n2 = cross2(normal, n1)
  stress = [getStress(sM, normal, normal), getStress(sM, normal, n1), getStress(sM, normal, n2)]
  for i in range(3):
    temp = 'int' + str(i+1)
    resultFull[num] = [final[0], final[1], el, i+1, db['cohesive'][el][temp][0], db['cohesive'][el][temp][2], db['cohesive'][el][temp][4], stress[0], stress[1], stress[2]]
    num += 1
  #print num
  if num%10000 == 0:
    print num

def getStress(sM, normal, vec):
  p = matmult(sM, normal)
  return dot(p, vec)

def matmult(m, v):
  rows = len(m)
  w = [0]*rows
  irange = range(len(v))
  sum = 0
  for j in range(rows):
    r = m[j]
    for i in irange:
      sum += r[i]*v[i]
    w[j],sum = sum,0
  return w
  
  return [v[0] - dot(v,n)*n[0], v[1] - dot(v,n)*n[1], v[2] - dot(v,n)*n[2]]

def getFirstDirection(n1, n2, n3):
  a = [(n2[0] - n1[0]), (n2[1] - n1[1]), (n2[2] - n1[2])]
  b = [(n3[0] - n1[0]), (n3[1] - n1[1]), (n3[2] - n1[2])]
  n = cross2(a, b)
  angle = getAngle([1,0,0], n)
  
  if angle <= 0.1:
    temp = projectVectorToPlane(a, b, [0,0,1])
    norm = sqrt(temp[0]**2 + temp[1]**2 + temp[2]**2)
    return [temp[0]/norm, temp[1]/norm, temp[2]/norm]
  else:
    temp = projectVectorToPlane(a, b, [1,0,0])
    norm = sqrt(temp[0]**2 + temp[1]**2 + temp[2]**2)
    return [temp[0]/norm, temp[1]/norm, temp[2]/norm]