Open
Description
Describe the bug
When the time step of an integration is fixed, eg by choosing some small number of time steps for some total evolution time, the integration appears to fail by returning the wrong wavefunction at long times. Crucially, it does so silently without throwing any errors. A fix would be to check the integrator status (eg status(integrator)==True
) or to adaptively choose a smaller step if needed.
To Reproduce
In https://github.com/Happy-Diode/Bloqade.jl/blob/master/examples/ring/main.jl reduce nsteps by factors of two until integration obviously fails