The propagate method is not time symmetric and therefore not unitary if internal and free propagators are combined.
We should change the propagate method from
freeprop(Delta t) -> internalprop(Delta t)
to
freeprop(Delta t/2) -> internalprop(Delta t) -> ifreeprop(Delta t/2)
https://github.com/bleykauf/aisim/blob/3e186208d59be49e3dcbe25c954299cb3477ddf6/aisim/prop.py#L35
The propagate method is not time symmetric and therefore not unitary if internal and free propagators are combined. We should change the propagate method from freeprop(Delta t) -> internalprop(Delta t) to freeprop(Delta t/2) -> internalprop(Delta t) -> ifreeprop(Delta t/2)