How to implement double-exponential functions for synapses in C++
by Ravichandran-Schmidt, Pirathitha (SRH Hochschule Heidelberg)
Dear NEST Community,
is there a good way to implement a double-exponential function for synapses in C++ within a nest module?
As I need to implement my own neuron model in C++, I adapted the code from extraction the information from other neuron models (for example: iaf _cond_beta.cpp) and used the following lines
f[S::DG_EXC] = -y[ S::DG_EXC ] / node.P_.tau_synAMPA_on;
f[S::G_EXC ] = y[ S::DG_EXC ] - y[ S::G_EXC ] / node.P_.tau_synAMPA_off
At first, it seemed to work pretty good. But when I swap the values for rise und decay time, the evolution of g_ADP did not change at all (for example: rise_time = 1ms, decay_time = 100ms).
I would be happy to get any feedback on it.