Hello all,
I have a system of two differential equations that I am able to solve using Rosenbrock23()
.
I now would like to model this system with one modification, where at time t
I manually change the value of one of the differential equation results, causing a discontinuity.
I thought that simply passing t
to d_discontinuties would handle this, but I run out of iterations, or the solver produces very unusual results. I currently have it manually set up to solve from 0
to t-1
then from t
to f
starting with values from t-1
, but that does not seem like an elegant solution and I get the following error:
┌ Warning: dt(-1.6e-9) <= dtmin(1.862645149230957e-9) at t=1.059225e7. Aborting. There is either an error in your model specification or the true solution is unstable.
Any tips for the proper way to handle this? My equations are stiff.