HomeostaticSynapse.h

Go to the documentation of this file.

#ifndef HOMEOSTATICSYNAPSE_H_
#define HOMEOSTATICSYNAPSE_H_

#include <string>
using std::string;
#include <boost/format.hpp>

#include "PCSIMException.h"
#include "SpikeBuffer.h"
#include "SimObject.h"
#include "SimNetwork.h"

#include "GenericStaticSpikingSynapse.h"
#include "GenericDynamicSpikingSynapse.h"
#include "GenericConductanceBasedSpikingSynapse.h"
#include "GenericCurrentBasedSpikingSynapse.h"
#include "GenericHomeostaticSynapse.h"
#include "GenericGlutamateSynapse.h"

#include "ExponentialDecaySpikeResponse.h"

// CONDUCTANCE BASED

class DynamicHomeostaticSynapse : public GenericHomeostaticSynapse<
                                         GenericDynamicSpikingSynapse< 
                                         GenericConductanceBasedSpikingSynapse< 
                                         ExponentialDecaySpikeResponse > > >
{
    SIMOBJECT(DynamicHomeostaticSynapse, AdvancePhase::SpikeDriven)

public:
    DynamicHomeostaticSynapse(
        const float Winit          = 1e-9,
        const float tau            = 5e-3,
        const float delay          = 2e-4,
        const float Erev           = 0e-3,
        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 double alpha         = 0.01,
        const double activity_tau  = 10.0,
        const double activity_goal = 1.5,
        const double activity_init = 1.0,
        const double Wex           = 10e-9,
        const double back_delay    = 0.0,
        const bool active          = true
        )
    {
        this->Winit = Winit;
        this->tau = tau;
        this->delay = delay;
        this->Erev = Erev;

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

        this->alpha = alpha;
        this->Wex=Wex;
        this->active=active;
        this->inhibitory=(Erev < -20e-3);

        this->activity_tau = activity_tau;
        this->activity_goal = activity_goal;
        this->activity_init = activity_init;
        this->back_delay = back_delay;
    };
};


class StaticHomeostaticSynapse : public GenericHomeostaticSynapse<
                                        GenericStaticSpikingSynapse< 
                                        GenericConductanceBasedSpikingSynapse<
                                        ExponentialDecaySpikeResponse > > >
{
    SIMOBJECT(StaticHomeostaticSynapse, AdvancePhase::SpikeDriven)

public:
    StaticHomeostaticSynapse(
        const float Winit          = 1e-9,
        const float tau            = 5e-3,
        const float delay          = 2e-4,
        const float Erev           = 0e-3,
        const double alpha         = 0.01,
        const double activity_tau  = 10.0,
        const double activity_goal = 1.5,
        const double activity_init = 1.0,
        const double Wex           = 10e-9,
        const double back_delay    = 0.0,
        const bool active          = true
        )
    {
        this->Winit = Winit;
        this->tau = tau;
        this->delay = delay;
        this->Erev = Erev;

        this->alpha = alpha;
        this->Wex=Wex;
        this->active=active;
        this->inhibitory=(Erev < -20e-3);

        this->activity_tau = activity_tau;
        this->activity_goal = activity_goal;
        this->activity_init = activity_init;
        this->back_delay = back_delay;
    };
};




// CURRENT BASED

class DynamicHomeostaticCurrentSynapse : public GenericHomeostaticSynapse<
                                                GenericDynamicSpikingSynapse< 
                                                GenericCurrentBasedSpikingSynapse< 
                                                ExponentialDecaySpikeResponse > > >
{
    SIMOBJECT(DynamicHomeostaticCurrentSynapse, AdvancePhase::SpikeDriven)

public:
    DynamicHomeostaticCurrentSynapse(
        const float Winit          = 1e-9,
        const float tau            = 5e-3,
        const float delay          = 2e-4,
        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 double alpha         = 0.01,
        const double activity_tau  = 10.0,
        const double activity_goal = 1.5,
        const double activity_init = 1.0,
        const double Wex           = 10e-9,
        const double back_delay    = 0.0,
        const bool active          = true
        )
    {
        this->Winit = Winit;
        this->tau = tau;
        this->delay = delay;

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

        this->alpha = alpha;
        this->Wex=Wex;
        this->active=active;
        this->inhibitory=(Winit < 0.0);

        this->activity_tau = activity_tau;
        this->activity_goal = activity_goal;
        this->activity_init = activity_init;
        this->back_delay = back_delay;
    };
};


class StaticHomeostaticCurrentSynapse : public GenericHomeostaticSynapse<
                                               GenericStaticSpikingSynapse< 
                                               GenericCurrentBasedSpikingSynapse<
                                               ExponentialDecaySpikeResponse > > >
{
    SIMOBJECT(StaticHomeostaticCurrentSynapse, AdvancePhase::SpikeDriven)

public:
    StaticHomeostaticCurrentSynapse(
        const float Winit          = 1e-9,
        const float tau            = 5e-3,
        const float delay          = 2e-4,
        const double alpha         = 0.01,
        const double activity_tau  = 10.0,
        const double activity_goal = 1.5,
        const double activity_init = 1.0,
        const double Wex           = 10e-9,
        const double back_delay    = 0.0,
        const bool active          = true
        )
    {
        this->Winit = Winit;
        this->tau = tau;
        this->delay = delay;

        this->alpha = alpha;
        this->Wex=Wex;
        this->active=active;
        this->inhibitory=(Winit < 0.0);

        this->activity_tau = activity_tau;
        this->activity_goal = activity_goal;
        this->activity_init = activity_init;
        this->back_delay = back_delay;
    };
};



typedef GenericHomeostaticSynapse<
        GenericDynamicSpikingSynapse< 
        GenericGlutamateSynapse< 
        ExponentialDecaySpikeResponse > > > HomDynGlutamateExpBase;

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

public:
    HomDynGlutamateExpSynapse(
        const float Winit     = 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 double alpha         = 0.01,
        const double activity_tau  = 10.0,
        const double activity_goal = 1.5,
        const double activity_init = 1.0,
        const double Wex           = 10e-9
       )
    {
        this->Winit = Winit;
        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->fract_NMDA = fract_NMDA;

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

        this->rescale = false;
        this->inhibitory=false;
        this->active=true;
        this->back_delay = 0.0;

        this->alpha = alpha;
        this->Wex=Wex;

        this->activity_tau = activity_tau;
        this->activity_goal = activity_goal;
        this->activity_init = activity_init;
    };

    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 HomDynGlutamateExpBase::reset(dt);
    }
};

#endif /*GENERICHOMEOSTATICSYNAPSE_H_*/

---