def METHOD_NAME(self, mo_coeff=None, s=None):
    if not self.use_ao_symmetry:
        raise RuntimeError("AO symmetry not initiated")
    from pyscf.scf.ghf_symm import METHOD_NAME
    if mo_coeff is None:
        mo_coeff = self.mo_coeff
    if s is None:
        s = self.get_ovlp()
    cell = self.cell
    symm_orb = cell.symm_orb
    irrep_id = cell.irrep_id
    orbsym = []
    for k in range(len(mo_coeff)):
        orbsym_k = np.asarray(METHOD_NAME(cell, mo_coeff[k], s=s[k],
                                         symm_orb=symm_orb[k], irrep_id=irrep_id[k]))
        orbsym.append(orbsym_k)
    return orbsym