def METHOD_NAME(self):
    mf = dft.ROKS(h2o_p)
    mf.run(conv_tol=1e-14, xc='lda,vwn')
    g = roks.Grad(mf)
    self.assertAlmostEqual(lib.fp(g.grad()), -0.12051785975616186, 6)
    g.grid_response = True
    self.assertAlmostEqual(lib.fp(g.grad()), -0.12052121736985746, 6)