def METHOD_NAME():
  sites_cart = matrix.col_list([
    (0,0,0),
    (1,0,0),
    (0.5,0.5*3**0.5,0)])
  com = matrix.col((0,0,0))
  for xyz in sites_cart:
    com += xyz
  com /= len(sites_cart)
  sites_cart = [xyz-com for xyz in sites_cart]
  E = RBDA_Eq_4_12(matrix.col((0.81, -0.32, -0.42, -0.26)).normalize())
  sites_cart = [E*xyz for xyz in sites_cart]
  return sites_cart