code: def get_consistent_set_and_saxis(magmoms, saxis=None):
          
          magmoms = [Magmom(magmom) for magmom in magmoms]
          if saxis is None:
              saxis = Magmom.get_suggested_saxis(magmoms)
          else:
              saxis = saxis/np.linalg.norm(saxis)
          magmoms = [magmom.get_moment(saxis=saxis) for magmom in magmoms]
          return (magmoms, saxis)
docs: Method to ensure a list of magmoms use the same spin axis.
          Returns a tuple of a list of Magmoms and their global spin axis.
  
          :param magmoms: list of magmoms (Magmoms, scalars or vectors)
          :param saxis: can provide a specific global spin axis
          :return: (list of Magmoms, global spin axis) tuple