def METHOD_NAME(self):
    full = rand_herm(2**4)
    beta = 0.624
    rhoth1 = thermal_state(full, beta)
    func = thermal_state(full, None, precomp_func=True)
    rhoth2 = func(beta)
    assert_allclose(rhoth1, rhoth2)