def METHOD_NAME(self):
    t1 = mcc.t1
    t2 = mcc.t2
    eris = mcc.ao2mo()
    nocc, nvir = t1.shape
    nmo = nocc + nvir
    mol = mcc.mol
    orbsym = symm.addons.label_orb_symm(mol, mol.irrep_id, mol.symm_orb,
                                        mcc.mo_coeff)
    orbsym = numpy.asarray(orbsym, dtype=numpy.int32)
    t2ref = t2.copy()
    mo_energy, t1T, t2T, vooo, foVT, restore_t2_inplace = \
            ccsd_t._sort_t2_vooo_(mcc, orbsym, t1, t2.copy(), eris)
    self.assertAlmostEqual(abs(t2ref-restore_t2_inplace(t2T.copy())).max(), 0, 12)
    o_sorted = numpy.hstack([numpy.where(orbsym[:nocc] == i)[0] for i in range(8)])
    v_sorted = numpy.hstack([numpy.where(orbsym[nocc:] == i)[0] for i in range(8)])
    o_sym = orbsym[o_sorted]
    oo_sym = (o_sym[:,None] ^ o_sym).ravel()
    oo_sorted = numpy.hstack([numpy.where(oo_sym == i)[0] for i in range(8)])
    ref_t2T = t2.transpose(2,3,1,0)
    ref_t2T = ref_t2T[v_sorted][:,v_sorted][:,:,o_sorted][:,:,:,o_sorted]
    ref_t2T = ref_t2T.reshape(nvir,nvir,-1)[:,:,oo_sorted].reshape(nvir,nvir,nocc,nocc)
    ref_vooo = numpy.asarray(eris.ovoo).transpose(1,0,2,3)
    ref_vooo = ref_vooo[v_sorted][:,o_sorted][:,:,o_sorted][:,:,:,o_sorted]
    ref_vooo = ref_vooo.reshape(nvir,-1,nocc)[:,oo_sorted].reshape(nvir,nocc,nocc,nocc)
    self.assertAlmostEqual(abs(ref_vooo-vooo).sum(), 0, 9)
    self.assertAlmostEqual(abs(ref_t2T-t2T).sum(), 0, 9)