def METHOD_NAME():
  atom_number = flex.int()
  coordinates = flex.vec3_double()
  for i in range (3):
    atom_number.append(i)
    coordinates.append( (i,i,i) )
  atom_types = multipolar.assign_atom_types( atom_number, coordinates )
  #print atom_types
  for atom_type in atom_types:
    print(atom_type)