Resolve "Exploit the property that hopping matrices are hermitian"

Closes #58 (closed) .The results are not spectacular since it seems that internally the libraries still call essentially zgemm. Since Op_Exp is only called twice for the whole runtime of the program I implemented Kahan's summation to get really accurate version of the hopping matrix.