I have a problem where the difference of two relative entropies is finite, but they individually diverge. The function that calculates relative entropies in qutip is qt.entropy_relative(state1, state2), however there is no option where I can increase the precision of the operation.