def METHOD_NAME(self):
    td = tdscf.TDA(mf_lda).run(singlet=False, nstates=nstates)
    tdg = td.nuc_grad_method()
    g1 = tdg.kernel(state=3)
    self.assertAlmostEqual(g1[0,2], -0.3311324654, 6)
    td_solver = td.as_scanner()
    pmol = mol.copy()
    e1 = td_solver(pmol.set_geom_('H 0 0 1.805; F 0 0 0', unit='B'))
    e2 = td_solver(pmol.set_geom_('H 0 0 1.803; F 0 0 0', unit='B'))
    self.assertAlmostEqual(abs((e1[2]-e2[2])/.002 - g1[0,2]).max(), 0, 4)