DEFAULT_TRACE_MODEL * SpikingGroup::get_post_trace( AurynFloat x )
{
for ( NeuronID i = 0 ; i < posttraces.size() ; i++ ) {
if ( posttraces[i]->get_tau() == x ) {
stringstream oss;
oss << get_name() << ":: Sharing post trace with " << x << "s timeconstant." ;
logger->msg(oss.str(),DEBUG);
return posttraces[i];
}
}
DEFAULT_TRACE_MODEL * tmp = new DEFAULT_TRACE_MODEL(get_post_size(),x);
posttraces.push_back(tmp);
return tmp;
}
void IF2Group::integrate_nonlinear_nmda_synapses()
{
// decay of ampa and gaba channel, i.e. multiply by exp(-dt/tau)
auryn_vector_float_scale(scale_ampa,g_ampa);
auryn_vector_float_scale(scale_gaba,g_gaba);
// compute dg_nmda = (g_ampa-g_nmda)*dt/tau_nmda and add to g_nmda
AurynFloat mul_nmda = dt/tau_nmda;
auryn_vector_float_saxpy(mul_nmda,g_ampa,g_nmda);
auryn_vector_float_saxpy(-mul_nmda,g_nmda,g_nmda);
// BEGIN implement NMDA voltage dependence
auryn_vector_float_copy( mem, nmda_opening);
auryn_vector_float_add_constant( nmda_opening , -e_nmda_onset );
auryn_vector_float_scale( nmda_slope, nmda_opening );
for ( AurynState * ptr = auryn_vector_float_ptr( nmda_opening , 0 ) ; ptr != auryn_vector_float_ptr( nmda_opening , get_post_size()-1 )+1 ; ++ptr ) {
AurynFloat x = *ptr;
AurynFloat x2 = x*x;
AurynFloat r = x2/(1.0+x2);
if (x>0) *ptr = r; // rectification
else *ptr = 0;
// cout << *ptr << endl;
}
// END implement NMDA voltage dependence
// excitatory
auryn_vector_float_copy(g_nmda,t_exc);
auryn_vector_float_scale(-A_nmda,t_exc);
auryn_vector_float_mul(t_exc,nmda_opening);
auryn_vector_float_saxpy(-A_ampa,g_ampa,t_exc);
auryn_vector_float_mul(t_exc,mem);
// inhibitory
auryn_vector_float_copy(mem,t_inh);
auryn_vector_float_add_constant(t_inh,-e_rev);
auryn_vector_float_mul(t_inh,g_gaba);
}
