Added tests for interacting Hamiltonians in test_time_solver.cc
where_to_truncate() with tolerance=0.0 now eliminates the trailing zeros.
Added some debug statements and computation of errors in trotter_unitary.cc
MPS time evolution now defaults to truncating at DBL_EPSILON