make res of trace_rho_op real if op + rho are hermitian

[the-hampel]

• add function to check if a triqs many body operator is hermitian
• check if input density matrix and operator to trace_rho_op are hermitian and return real part of averaged value

[the-hampel]

We could expose the function is_op_hermitian to python?