def METHOD_NAME(self):
    vsize = kci.vector_size()
    vec = np.random.rand(vsize) * 1j
    hc = kci.matvec(vec, 1, eris=eris)
    h = ci.kcis_rhf.cis_H(kci, 1, eris=eris)
    ref = h.dot(vec)
    self.assertAlmostEqual(abs(hc - ref).max(), 0, 7)