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