code: def get_suggested_saxis(magmoms):
          
          # heuristic, will pick largest magmom as reference
          # useful for creating collinear approximations of
          # e.g. slightly canted magnetic structures
          # for fully collinear structures, will return expected
          # result
  
          magmoms = [Magmom(magmom) for magmom in magmoms]
          # filter only non-zero magmoms
          magmoms = [magmom for magmom in magmoms if abs(magmom)]
          magmoms.sort(reverse=True)
          if len(magmoms) > 0:
              return magmoms[0].get_00t_magmom_with_xyz_saxis().saxis
          else:
              return np.array([0, 0, 1], dtype="d")
docs: This method returns a suggested spin axis for a set of magmoms,
          taking the largest magnetic moment as the reference. For calculations
          with collinear spins, this would give a sensible saxis for a ncl
          calculation.
  
          :param magmoms: list of magmoms (Magmoms, scalars or vectors)
          :return: np.ndarray of length 3