GlutamateSynapse.h

Go to the documentation of this file.

#ifndef _GLUTAMATESYNAPSE_H_
#define _GLUTAMATESYNAPSE_H_

#include "GenericGlutamateSynapse.h"

#include "GenericDynamicSpikingSynapse.h"
#include "GenericStaticSpikingSynapse.h"

#include "ExponentialDecaySpikeResponse.h"
#include "DoubleExponentialSpikeResponse.h"
#include "AlphaFunctionSpikeResponse.h"

#include "GenericHomeostaticSynapse.h"


// NEEDS TO BE TESTED

// Dynamic Synapses

// Exponential Response

typedef GenericDynamicSpikingSynapse< GenericGlutamateSynapse< ExponentialDecaySpikeResponse > > DynamicGlutamateExpBase;

class DynamicGlutamateExpSynapse : public GenericDynamicSpikingSynapse< GenericGlutamateSynapse< ExponentialDecaySpikeResponse > >
{
    SIMOBJECT(DynamicGlutamateExpSynapse, AdvancePhase::SpikeDriven)

public:
    DynamicGlutamateExpSynapse(
        const float W          = 1e-9,
        const float tau_AMPA   = 5e-3,
        const float tau_NMDA   = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const float U          = 0.2,
        const float D          = 0.1,
        const float F          = 0.4,
        const float u0         = 0.2,
        const float r0         = 1.0,
        const float f0         = -1.0,
        const int  has_NMDA    = 1,
        const bool rescale     = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau_NMDA = tau_NMDA;
        this->tau_AMPA = tau_AMPA;

        this->tau = tau_NMDA;
        this->resp_AMPA.tau = tau_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->U  = U;
        this->D  = D;
        this->F  = F;
        this->u0 = u0;
        this->r0 = r0;
        this->f0 = f0;

        this->rescale = rescale;
    };

    double tau_AMPA;

    double tau_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau = tau_NMDA;
        resp_AMPA.tau = tau_AMPA;

        if(rescale)
            scale_NMDA = tau_AMPA/tau_NMDA;
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return DynamicGlutamateExpBase::reset(dt);
    }
};


typedef GenericDynamicSpikingSynapse < GenericGlutamateSynapse < AlphaFunctionSpikeResponse > >
        DynamicGlutamateAlphaBase;

class DynamicGlutamateAlphaSynapse : public GenericDynamicSpikingSynapse < GenericGlutamateSynapse < AlphaFunctionSpikeResponse > >
{
    SIMOBJECT(DynamicGlutamateAlphaSynapse, AdvancePhase::SpikeDriven)

public:
    DynamicGlutamateAlphaSynapse(
        const float W          = 1e-9,
        const float tau_AMPA   = 5e-3,
        const float tau_NMDA   = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const float U          = 0.2,
        const float D          = 0.1,
        const float F          = 0.4,
        const float u0         = 0.2,
        const float r0         = 1.0,
        const float f0         = -1.0,
        const int  has_NMDA    = 1,
        const bool rescale     = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau_NMDA = tau_NMDA;
        this->tau_AMPA = tau_AMPA;

        this->tau = tau_NMDA;
        this->resp_AMPA.tau = tau_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->U  = U;
        this->D  = D;
        this->F  = F;
        this->u0 = u0;
        this->r0 = r0;
        this->f0 = f0;

        this->rescale = rescale;
    };

    double tau_AMPA;

    double tau_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau = tau_NMDA;
        resp_AMPA.tau = tau_AMPA;

        if(rescale)
            scale_NMDA = tau_AMPA/tau_NMDA;
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return DynamicGlutamateAlphaBase::reset(dt);
    }
};


// Double Exponential Response

typedef GenericDynamicSpikingSynapse< GenericGlutamateSynapse< DoubleExponentialSpikeResponse > >
        DynamicGlutamateDoubleExpBase;

class DynamicGlutamateDoubleExpSynapse : public GenericDynamicSpikingSynapse< GenericGlutamateSynapse< DoubleExponentialSpikeResponse > >
{
    SIMOBJECT(DynamicGlutamateDoubleExpSynapse, AdvancePhase::SpikeDriven)

public:
    DynamicGlutamateDoubleExpSynapse(
        const float W          = 1e-9,
        const float tau1_AMPA  = 0.1e-3,
        const float tau2_AMPA  = 5e-3,
        const float tau1_NMDA  = 2e-3,
        const float tau2_NMDA  = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const float U          = 0.2,
        const float D          = 0.1,
        const float F          = 0.4,
        const float u0         = 0.2,
        const float r0         = 1.0,
        const float f0         = -1.0,
        const int  has_NMDA    = 1,
        const bool rescale     = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau1_NMDA = tau1_NMDA;
        this->tau2_NMDA = tau2_NMDA;
        this->tau1_AMPA = tau1_AMPA;
        this->tau2_AMPA = tau2_AMPA;

        this->tau1 = tau1_NMDA;
        this->tau2 = tau2_NMDA;
        this->resp_AMPA.tau1 = tau1_AMPA;
        this->resp_AMPA.tau2 = tau2_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->U  = U;
        this->D  = D;
        this->F  = F;
        this->u0 = u0;
        this->r0 = r0;
        this->f0 = f0;

        this->rescale = rescale;
    };

    double tau1_AMPA;

    double tau2_AMPA;

    double tau1_NMDA;

    double tau2_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau1 = tau1_NMDA;
        tau2 = tau2_NMDA;
        resp_AMPA.tau1 = tau1_AMPA;
        resp_AMPA.tau2 = tau2_AMPA;

        if(rescale)
            scale_NMDA = (tau1_AMPA+tau2_AMPA)/(tau1_NMDA+tau2_NMDA);
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return DynamicGlutamateDoubleExpBase::reset(dt);
    }
        
};



// Static Synapses

// Exponential Response

typedef GenericStaticSpikingSynapse< GenericGlutamateSynapse< ExponentialDecaySpikeResponse > > StaticGlutamateExpBase;

class StaticGlutamateExpSynapse : public GenericStaticSpikingSynapse< GenericGlutamateSynapse< ExponentialDecaySpikeResponse > >
{
    SIMOBJECT(StaticGlutamateExpSynapse, AdvancePhase::SpikeDriven)

public:
    StaticGlutamateExpSynapse(
        const float W          = 1e-9,
        const float tau_AMPA   = 5e-3,
        const float tau_NMDA   = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const int  has_NMDA    = 1,
        const bool rescale     = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau_NMDA = tau_NMDA;
        this->tau_AMPA = tau_AMPA;

        this->tau = tau_NMDA;
        this->resp_AMPA.tau = tau_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->rescale = rescale;
    };

    double tau_AMPA;

    double tau_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau = tau_NMDA;
        resp_AMPA.tau = tau_AMPA;

        if(rescale)
            scale_NMDA = tau_AMPA/tau_NMDA;
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return StaticGlutamateExpBase::reset(dt);
    }
};


typedef GenericStaticSpikingSynapse < GenericGlutamateSynapse < AlphaFunctionSpikeResponse > >
        StaticGlutamateAlphaBase;

class StaticGlutamateAlphaSynapse : public GenericStaticSpikingSynapse < GenericGlutamateSynapse < AlphaFunctionSpikeResponse > >
{
    SIMOBJECT(StaticGlutamateAlphaSynapse, AdvancePhase::SpikeDriven)

public:
    StaticGlutamateAlphaSynapse(
        const float W          = 1e-9,
        const float tau_AMPA   = 5e-3,
        const float tau_NMDA   = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const int   has_NMDA   = 1,
        const bool  rescale    = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau_NMDA = tau_NMDA;
        this->tau_AMPA = tau_AMPA;

        this->tau = tau_NMDA;
        this->resp_AMPA.tau = tau_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->rescale = rescale;
    };

    double tau_AMPA;

    double tau_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau = tau_NMDA;
        resp_AMPA.tau = tau_AMPA;

        if(rescale)
            scale_NMDA = tau_AMPA/tau_NMDA;
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return StaticGlutamateAlphaBase::reset(dt);
    }
};


// Double Exponential Response

typedef GenericStaticSpikingSynapse< GenericGlutamateSynapse< DoubleExponentialSpikeResponse > >
        StaticGlutamateDoubleExpBase;

class StaticGlutamateDoubleExpSynapse : public GenericStaticSpikingSynapse< GenericGlutamateSynapse< DoubleExponentialSpikeResponse > >
{
    SIMOBJECT(StaticGlutamateDoubleExpSynapse, AdvancePhase::SpikeDriven)

public:
    StaticGlutamateDoubleExpSynapse(
        const float W          = 1e-9,
        const float tau1_AMPA  = 0.1e-3,
        const float tau2_AMPA  = 5e-3,
        const float tau1_NMDA  = 2e-3,
        const float tau2_NMDA  = 150e-3,
        const float delay      = 2e-4,
        const float Erev_AMPA  = 0e-3,
        const float Erev_NMDA  = 0e-3,
        const float fract_NMDA = 0.5,
        const int  has_NMDA    = 1,
        const bool rescale     = true)
    {
        this->W = W;
        this->Mg_conc = 1e-3;
        this->delay = delay;

        this->tau1_NMDA = tau1_NMDA;
        this->tau2_NMDA = tau2_NMDA;
        this->tau1_AMPA = tau1_AMPA;
        this->tau2_AMPA = tau2_AMPA;

        this->tau1 = tau1_NMDA;
        this->tau2 = tau2_NMDA;
        this->resp_AMPA.tau1 = tau1_AMPA;
        this->resp_AMPA.tau2 = tau2_AMPA;

        this->Erev_NMDA = Erev_NMDA;
        this->Erev_AMPA = Erev_AMPA;

        this->has_NMDA = has_NMDA;
        this->fract_NMDA = fract_NMDA;

        this->rescale = rescale;
        this->has_NMDA = has_NMDA;
    };

    double tau1_AMPA;

    double tau2_AMPA;

    double tau1_NMDA;

    double tau2_NMDA;

    bool rescale;

    virtual int reset(double dt)
    {
        tau1 = tau1_NMDA;
        tau2 = tau2_NMDA;
        resp_AMPA.tau1 = tau1_AMPA;
        resp_AMPA.tau2 = tau2_AMPA;

        if(rescale)
            scale_NMDA = (tau1_AMPA+tau2_AMPA)/(tau1_NMDA+tau2_NMDA);
        else
            scale_NMDA = 1.0;

        resp_AMPA.reset(dt);
        return StaticGlutamateDoubleExpBase::reset(dt);
    }

};


#endif //_GLUTAMATESYNAPSE_H_

---