I’d like to use TerminateSteadyState from DiffEqCallbacks.jl on a DiscreteProblem, but the default test function is based on derivatives (which my model doesn’t have). I see that one can define their own test function based on the integrator - any tips on how this should be structured in order to correspond to the ‘allDerivPass’ function for continuous systems?
Yup - I just need to use the difference in states uprev and u - I didn’t know whether there was a canned function for this, or how to build the callback appropriately.
Looking closer, allDerivPass appears to be able to handle DiscreteProblem, but I get an error that the FunctionMap stepper doesn’t support derivatives, and the callback fails.