def METHOD_NAME(self):
    numpy.random.seed(1)
    nao = mol.nao_nr()
    dms = numpy.random.random((2,nao,nao))
    mf = scf.RHF(mol).density_fit(auxbasis='weigend')
    vj0 = mf.get_j(mol, dms)
    vj1 = mf.get_jk(mol, dms)[0]
    self.assertAlmostEqual(abs(vj0-vj1).max(), 0, 12)
    self.assertAlmostEqual(lib.fp(vj0), -194.15910890730052, 9)