def METHOD_NAME(self):
    mf1 = copy.copy(mf)
    nao, nmo = mf.mo_coeff[0].shape
    numpy.random.seed(10)
    mf1.mo_coeff = numpy.random.random((2,nao,nmo))
    numpy.random.seed(12)
    nocca, noccb = mol.nelec
    nmo = mf1.mo_occ[0].size
    nvira = nmo - nocca
    nvirb = nmo - noccb
    t1a  = .1 * numpy.random.random((nocca,nvira))
    t1b  = .1 * numpy.random.random((noccb,nvirb))
    t2aa = .1 * numpy.random.random((nocca,nocca,nvira,nvira))
    t2aa = t2aa - t2aa.transpose(0,1,3,2)
    t2aa = t2aa - t2aa.transpose(1,0,2,3)
    t2bb = .1 * numpy.random.random((noccb,noccb,nvirb,nvirb))
    t2bb = t2bb - t2bb.transpose(0,1,3,2)
    t2bb = t2bb - t2bb.transpose(1,0,2,3)
    t2ab = .1 * numpy.random.random((nocca,noccb,nvira,nvirb))
    mycc = cc.GCCSD(mf1)
    t1 = mycc.spatial2spin((t1a, t1b        ))
    t2 = mycc.spatial2spin((t2aa, t2ab, t2bb))
    eris = mycc.ao2mo()
    e3a = gccsd_t.kernel(mycc, eris, t1, t2)
    self.assertAlmostEqual(e3a, 9877.2780859693339, 5)