How to add Callbacks for FODEProblem solvers?

[junaidaliop]

I am looking to find fractional-order solutions of Izhikevich neuron?
Are callbacks such as the one implemented here:

https://docs.sciml.ai/SciMLTutorialsOutput/html/models/08-spiking_neural_systems.html

function thr(u,t,integrator)
    integrator.u[1] >= 30
end

function reset!(integrator)
    integrator.u[1] = integrator.p[3]
    integrator.u[2] += integrator.p[4]
end

threshold = DiscreteCallback(thr,reset!)
current_step= PresetTimeCallback(50,integrator -> integrator.p[5] += 10)
cb = CallbackSet(current_step,threshold)

be imported with

prob = FODEProblem(izh!, alpha, u0, tspan, p, callback=cb)

I think that the callback doesn't intregate with FODEProblem.
What would be a way around this?

Also, which solver would be recommended to use with this?

[ErikQQY]

Sorry, the current state of FractionalDiffEq doesn't support discrete or continuous callbacks, I think you might have to work on the rhs part to implement similar functionalities in FractionalDiffEq.

[junaidaliop]

I would appreciate if you could provide me with a starting point or guidance to get started on this.