I ref this “Handling Instability When Solving ODE Problems - #5 by ChrisRackauckas”,Thanks.
But I found in solver option ’ * d_discontinuities:
Denotes locations of discontinuities in low order derivatives. This will force FSAL algorithms which assume derivative continuity to re-evaluate the derivatives at the point of discontinuity. The default is []
.',can you give some example how to use this d_discontinuities
add these code
function condition(u,t,integrator)
t-t_discontinue
end
function affect!(integrator)
println("u1=",integrator.u[1],"t=",integrator.t)
end
cb = ContinuousCallback(condition,affect!,save_positions=(false,false))
'''