I was trying to use the one center TDVP algorithm to study a system with a time-dependent term that is not part of the Hamiltonian. The overall system evolution is described by
$$\dot y = H y + f(t)$$
where y is the state, H the Hamiltonian, and f(t) is a time dependent vector. How should I modify the original TDVP code to treat this case?