aEIFNeuron.cpp

Go to the documentation of this file.

#include "aEIFNeuron.h"

#include <cmath>
#include <iostream>
using std::cout;
using std::endl;


// aEIFNeuronBase

ThreadSpecificRandomDistribution< NormalDistribution > aEIFNeuronBase::white_noise;


aEIFNeuronBase::aEIFNeuronBase(double a, double b, double Vt, double Vr,  double El, double gl,
    double Cm, double tau_w, double slope, double Vpeak, double Vinit, double Inoise, double Iinject)
    : a(a), b(b), Vt(Vt), Vr(Vr), El(El), gl(gl), Cm(Cm), tau_w(tau_w), slope(slope)
    , Vpeak(Vpeak), Vinit(Vinit), Inoise(Inoise), Iinject(Iinject), w(0.0), Isyn(0.0)
{
}


int aEIFNeuronBase::advance(AdvanceInfo const &ai)
{
    bool register hasFired = false;

    if (Vm >= Vpeak)
        Vm = Vr;

    double V=Vm;

    // all synapses have added the contributions to Isyn, we add Iinject    
    double Itot = Isyn + Iinject;

    if (Inoise > 0.0) {
        Itot += (white_noise() * Inoise);
//      else if(noiseType == 1)
//          Itot += (pink_noise() * Inoise);
    }

    Vm = V + _dt*1/Cm*(-gl*(V-El) + gl*slope*exp((V-Vt)/slope) - w + Itot);
    w = w + _dt*1/tau_w*(a*(V-El) - w);

    if (Vm >= Vpeak) {         // Note that the neuron has fired!
        Vm = 50e-3;
        w = w + b;
        hasFired = true;
    }

    // clear synaptic input for next time step
    clearSynapticInput();

    if (hasFired) {
        out_port.setSpike(ai);
        return ADVANCEFLAG_HASSPIKED;
    }

    return 0;
}


inline void aEIFNeuronBase::clearSynapticInput(void)
{
    Isyn = 0.0;
}


int aEIFNeuronBase::adjust(double dt)
{
    _dt = dt;

//    pink_noise.adjust(dt);
    clearSynapticInput();

    return 0;
}


int aEIFNeuronBase::reset(double dt)
{
    SingleOutputSpikeSender::reset();

    Vm = Vinit;     // set membrane voltage to its initial value
    w = 0.0;
    clearSynapticInput();
    adjust(dt);

    return 0;
}



// aEIFNeuron

aEIFNeuron::aEIFNeuron(double a, double b, double Vt, double Vr, double El, double gl, 
    double Cm, double tau_w, double slope, double Vpeak, double Vinit, double Inoise, double Iinject)
    : aEIFNeuronBase(a, b, Vt, Vr, El, gl, Cm, tau_w, slope, Vpeak, Vinit, Inoise, Iinject)
{
}


inline void aEIFNeuron::currentInput(double i)
{
    Isyn += i;
}



// CbaEIFNeuron

CbaEIFNeuron::CbaEIFNeuron(double a, double b, double Vt, double Vr, double El, double gl, 
    double Cm, double tau_w, double slope, double Vpeak, double Vinit, double Inoise, double Iinject)
    : aEIFNeuronBase(a, b, Vt, Vr, El, gl, Cm, tau_w, slope, Vpeak, Vinit, Inoise, Iinject)
{
}

inline void CbaEIFNeuron::currentInput(double i)
{
    Isyn += i;
}


inline void CbaEIFNeuron::conductanceInput(double g, double Erev)
{
    Isyn += g*(Erev-Vm);
}

---