Add the `return_eigenvectors` parameter to the call of `operator_spectrum`

Could be useful from user perspective, when calculating the Tr and Tr \ln terms, to optimize the CPU time vs precision ratio.

@parras @pfrank