def METHOD_NAME(self):
    td = mf_lda.TDA()
    es = td.kernel(nstates=5)[0]
    a,b = td.get_ab()
    nocc, nvir = a.shape[:2]
    nov = nocc * nvir
    e_ref = numpy.linalg.eigh(a.reshape(nov,nov))[0]
    self.assertAlmostEqual(abs(es[:3]-e_ref[:3]).max(), 0, 5)
    self.assertAlmostEqual(lib.fp(es[:3] * 27.2114), 3.220469947746697, 5)