GenericDynamicSpikingSynapse.h

Go to the documentation of this file.

#ifndef GENERICDYNAMICSPIKINGSYNAPSE_H_
#define GENERICDYNAMICSPIKINGSYNAPSE_H_

#include "SimObject.h"

#include "ExponentialDecaySpikeResponse.h"
#include "AlphaFunctionSpikeResponse.h"
#include "DoubleExponentialSpikeResponse.h"
#include "GABA_BResponse.h"
#include "GABA_BResponseODE.h"
#include "GenericNonLinGABASynapse.h"
#include "AlphaInfiniteSpikeResponse.h"

#include "GenericCurrentBasedSpikingSynapse.h"
#include "GenericConductanceBasedSpikingSynapse.h"

#include "GenericGlutamateSynapse.h"
#include "GenericGABASynapse.h"
#include "GenericNMDASynapse.h"

template<class SpikeResponse>
class GenericDynamicSpikingSynapse : public SpikeResponse
{
public:
    GenericDynamicSpikingSynapse() : W(1e-9f), delay(2e-4f), U(0.2f), D(0.1f), F(0.4f), u0(0.2f), r0(1.0f)
    { /* NOOP */ };

    virtual ~GenericDynamicSpikingSynapse()
    { /* NOOP */ };

    float W;


    float delay;
    virtual double getManagedDelay() const { return delay; };

    float U;
    float D;
    float F;
    float u0;
    float r0;
    float f0;


    virtual int reset( double dt )
    {
        if (f0 >= 0.0) {
            u0 = U / (1.0 - (1.0-U) * exp(-1.0/(f0*F)) );
            r0 = (1.0 - exp(-1.0/(f0*D))) / (1.0 - (1.0-u0) * exp(-1.0/(f0*D)) );
        }

        u         = u0;
        r         = r0;
        lastSpike = -1;
        SpikeResponse::reset( dt );

        return 0;
    }

    virtual int spikeHit( spikeport_t port, SpikeEvent const& spike )
    {
        if ( lastSpike > 0 ) {
            float isi = (float)( spike.time() - lastSpike );
            r = 1.0 + (r*(1.0 - u) - 1.0)*exp(-isi/D);
            u = U + u*(1.0 - U)*exp(-isi/F);
        }

        lastSpike = (float)spike.time();

        return SpikeResponse::spikeHit( port, SpikeEvent( spike, W * u * r ) );
    }

private:
    float u;

    float r;

    float lastSpike;
};

#endif /*GENERICDYNAMICSPIKINGSYNAPSE_H_*/

---