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.

Thanks,

Pira