def METHOD_NAME(self):
    mf = scf.UHF(mol).density_fit ()
    mf.conv_tol = 1e-14
    e0 = mf.kernel()
    g = mf.nuc_grad_method ().kernel()
    mf_scanner = mf.as_scanner()
    e1 = mf_scanner('''O    0.   0.       0.
                    1    0.   -0.758   0.587
                    1    0.   0.757    0.587''')
    e2 = mf_scanner('''O    0.   0.       0.
                    1    0.   -0.756   0.587
                    1    0.   0.757    0.587''')
    self.assertAlmostEqual(g[1,1], (e2-e1)/2e-3*lib.param.BOHR, 5)
    e1 = mf_scanner('''O    0.   0.       0.
                    1    0.   -0.7571  0.587
                    1    0.   0.757    0.587''')
    e2 = mf_scanner('''O    0.   0.       0.
                    1    0.   -0.7569 0.587
                    1    0.   0.757    0.587''')
    self.assertAlmostEqual(g[1,1], (e2-e1)/2e-4*lib.param.BOHR, 7)
    mf = scf.UHF(mol1).density_fit ()
    mf.conv_tol = 1e-14
    e0 = mf.kernel()
    g = mf.nuc_grad_method ().kernel()
    mf_scanner = mf.as_scanner()
    e1 = mf_scanner('''