changed MPO and HeisenbergMPO to have complex prefactors

because this is needed for the new "AFMPO" class of hte auxiliary field operators.

I had to change the couplingmatrix to complex as well and this caused some changes in TrotterHeisenberg.
Everything ok now.