Luca Gravina

icon indicating a website link https://inspirehep.net/authors/2134031

icon company EPFL icon location Lausanne, Switzerland icon location Doctoral Researcher at EPFL, Lausanne (Switzerland) Quantum Optics, Neural Quantum States, Condensed Matter Physics, Cavity QED

Results 3 comments of Luca Gravina

Rewinding integrator

Indeed `reinit!` serves the purpose but it has one problem that might be an interesting issue. Consider the following simple problem: ``` u0 = 1 tl = LinRange(0,1,1001) tspan =...

Rewinding integrator

I tried solving by using ``` function affect_1!(integrator) integrator.u = sol.u[select(t0, tl)] integrator.t = t0 reinit!(integrator, integrator.u; t0=t0, erase_sol=false) idx = findlast(integrator.sol.t .

Rewinding integrator

This is perfect. Indeed it works. The only problem now is that if one has a PresetTimeCallback this interferes with the reset of the integrator. Specifically, if I take again...