I’m trying to replicate some previous work (eg: Naud & Sprekeler 2018) that has a “square” signal in response to spiking. I had two thoughts on going about this. My first thought was to use
not_refractory, but the signal needs to be delayed (“on” if 0.5ms < delta < 3ms, 0 otherwise). My second idea was to put a conditional, compared to
lastspike in the derivative equation, eg:
eqsPyr = ''' ... dx/dt = (-(Cd/tau_x)*(x-El) + (gd / (1+exp(-(x-Ed)/Dd))) + gc*( ((t-lastspike)< (3*ms)) * ((t-lastspike)>0.5)) : '''
However, it seems like 1) ‘lastspike’ isn’t known within the update space, and 2) comparisons (<) can’t be made in this space either.
Does anyone have any thoughts on how I might implement this square signal?