Speed improvement for herm=1

I don't know why but np.real produces big overhead. For small matrices cy_expect_rho_vec with herm=1 was two times slower than herm=0. Now it is almost as fast.