Matt Fishman

Could you be more explicit with some code (for example pointing out problematic lines of code or showing a minimal example)? Is this an issue within the `tdvp` function, or...

For example you should be able to do: ```julia tdvp(H, ψ0, -0.1im; nsweeps=10) ``` in order to evolve to time `1.0`, does that address the problem you are seeing?

Currently the time that is input is the timestep per sweep, which was an arbitrary choice. We can change it so that the input time `t` is the total time,...

Right, it sounds like a simple way to go is to have a high level interface `tdvp` that does multiple sweeps of TDVP up to a desired number of sweeps...

We have a package: https://github.com/GTorlai/Observers.jl which I think is pretty simple to use. Basically you just create an `Observer` from a list of functions you want called at each iteration,...

The goal was to generalize the observer functionality in ITensor with just a single type, and ultimately the hope was to just use that `Observer` type in `dmrg`. So `tdvp`...

Yeah, it would be interesting to see if the noise term helps, I haven't seen that technique used in the context of TDVP.

What about at the beginning of the evolution, in cases that people use subspace expansion to avoid going along a bad trajectory (https://arxiv.org/abs/2005.06104)? I guess the noise may just add...

One reason I like the second one is that it is more similar to the equation `exp(Ht)|psi>` (in fact we could consider calling `tdvp` something like `exp` or `exponentiate`, in...

Yeah, I think for the sake of `ITensorsGPU.jl` my conclusion is that we can just overload the highest level `contract!!` function, which is the one that actually gets called by...