Implementation of Izhikevich neuron model

[adam]

Hi Iann.
The short answer is both reset conditions can be expressed within an equation string:

reset_eq = 'v = c; u += d;'
N = NeuronGroup(..., reset = reset_eq)

the long answer is there are some unclear details when it comes to implementing the Izhikevich model, you can see some recent discussion in this thread:

Realization the model of network in Brain 2

[Screen Shot 2022-01-13 at 11.01.10 AM] [Screen Shot 2022-01-13 at 11.01.24 AM] the above are from my run of the code presented so far, with some modified plotting. I added a 40Hz sine to compare to the gamma rhythm mentioned in the original paper: [Screen Shot 2022-01-13 at 11.03.07 AM] visually inspecting the rasters, the oscillation seems less pronounced (and less synchronous), but looking at the rates, there does appear to be a 40Hz oscillation in the inhibitory firing rate. here is m…