def METHOD_NAME(mol, *args):
    if mol.nelectron == 1:
        if not mol.symmetry or mol.groupname == 'C1':
            return rohf.HF1e(mol)
        else:
            return hf_symm.HF1e(mol, *args)
    elif not mol.symmetry or mol.groupname == 'C1':
        return rohf.METHOD_NAME(mol, *args)
    else:
        return hf_symm.METHOD_NAME(mol, *args)