def METHOD_NAME(self):
    gref = mol.RKS(xc='wb97').run().nuc_grad_method().kernel()
    g1 = mol.RKS(xc='wb97').density_fit().run().nuc_grad_method().kernel()
    self.assertAlmostEqual(abs(gref - g1).max(), 0, 4)