Dear Charl,
Thanks for your suggestions!
We were able to get the log values of the synaptic trace we were interested in by running
the simulation in chunks.
However, we noticed that the instruction we wrote in the update block were never executed,
not even when receiving a presynaptic spike! Thus, we were only able to have a
"squared" exponential decay: the trace value was only determined by the equation
```syn_trace' = syn_trace / tau_syn ``` and was only updated in big steps when a spike
travelled through the individual synapse (note that we are currently using a toy network
with four neurons and very few spikes).
Here reported a figure (
https://imgur.com/a/AJfbcu1 ) to explain the concept: the four
lines are the exponentially-decaying synaptic traces for each connection in the network
and the markers are the spikes travelling from the presynaptic neurons; if we look at the
blue trace, for instance, we see that its value is updated only when the yellow-marked
spikes are present, causing an irregular shape of the exponential decay.
We might highlight the fact that we are using NESTML 5.0.0rc2 since NESTML 5.0.0 has
conflicts with the existing NEST installation on our workstation (version 3.1) that we
cannot change due to shared usage of the machine.
With this being said, what would be the best way to have a smooth decay of an exponential
function that depends the least on (or possibly is independent of) the actual spikes being
processed by the synapse model? (i.e. how could we achieve curves similar to the dopamine
traces shown in the tutorial you sent us?)
Best,
Léa D.