DifferentialEquations: Enforce that variables are inside a domain (isoutofdomain, PositiveDomain)

[ChrisRackauckas]

@LolianSh

I was thinking of simply clamping the array u after the integration step, in order to ensure my bounds. This way I would also avoid defining the differential equations outside the domain of my variables. Do you think this would be a reasonable solution ?

For a fixed time step that is reasonable. All of the other strategies are made for adaptive integration.

@LolianSh

So I would like to ask how I can manually change the state of integrator.u right after it has been updated, before going to the next integration step and before saving the result to the output object

DiscreteCallback. Make condition(u,t,integrator) = true so it happens after every step, and then affect!(integ) = (integ.u[...] = ...). You may then want to set save_positions = (false,false) in the construction of the callback so it doesn’t add extra save points.