Time dependent hamiltonians in 4op 2t correlator

Modified 4op 2t correlator for hamiltonians with time dependence; works with cython or function td format. Added relevant test functions.
This Hamiltonian transformation is not specific to the 4op_2t solver, is it? Perhaps the function should be called something like `_transform_H_t_shift` instead?