docs: def centre_of_atoms(atoms, mass_weighted=True):
      """ Returns centre point of any list of atoms.
  
      Parameters
      ----------
      atoms : list
          List of AMPAL atom objects.
      mass_weighted : bool, optional
          If True returns centre of mass, otherwise just geometric centre of points.
  
      Returns
      -------
      centre_of_mass : numpy.array
          3D coordinate for the centre of mass.
      """
      points = [x._vector for x in atoms]
      if mass_weighted:
          masses = [x.mass for x in atoms]
      else:
          masses = []
      return centre_of_mass(points=points, masses=masses)