Hi,
I am new to Julia and my aim is to create a network of spiking neurons. I use simple leaky integrate and fire (LIF) model for now with ODEProblem solver. Whenever the neuron (u) reaches a threshold value, it should enter a refractory period of some fixed ms, where du would be 0 and u would be kept in a reset value. Also the time points at which u reaches the threshold must be recorded as a vector for each neuron. However I couldn`t really understand how I can implement this, I use callbacks for thresholds passing condition, I am not sure how to use callbacks for setting refractory time. In python I would use a counter inside the function but since ODE handles the time itself I am confused how to intervene with it.