def METHOD_NAME(self):
    eri4 = ao2mo.restore(4, mf._eri, nmo)
    dm = mf.make_rdm1()
    vj0, vk0 = _vhf.incore(eri4, dm, hermi=1)
    vj1, vk1 = scf.hf.get_jk(mol, dm, hermi=1)
    self.assertTrue(numpy.allclose(vj0,vj1))
    self.assertTrue(numpy.allclose(vk0,vk1))