def METHOD_NAME(self):
    numpy.random.seed(1)
    mf = scf.RHF(mol).density_fit(auxbasis='weigend')
    nao = mol.nao_nr()
    dms = numpy.random.random((2,nao,nao))
    vj, vk = mf.get_jk(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vj), -194.15910890730066, 9)
    self.assertAlmostEqual(lib.fp(vk), -46.365071587653517, 9)
    vj = mf.get_j(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vj), -194.15910890730066, 9)
    vk = mf.get_k(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vk), -46.365071587653517, 9)
    mf.with_df = None
    vj, vk = mf.get_jk(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vj), -194.08878302990749, 9)
    self.assertAlmostEqual(lib.fp(vk), -46.530782983591152, 9)
    vj = mf.get_j(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vj), -194.08878302990749, 9)
    vk = mf.get_k(mol, dms, hermi=0)
    self.assertAlmostEqual(lib.fp(vk), -46.530782983591152, 9)