aEIFNeuron.h

Go to the documentation of this file.

#ifndef _AEIFNEURON_H_
#define _AEIFNEURON_H_

#include "InputTargetTypes.h"
#include "SpikeSender.h"
#include "SimObject.h"
#include "ThreadSpecificRandomDistribution.h"
//#include "PinkGenerator.h"


class aEIFNeuronBase : public SimObject, public SingleOutputSpikeSender
{
public:
    aEIFNeuronBase(double a       = 4e-9,       // S
                   double b       = 0.0805e-9,  // A
                   double Vt      = -50.6e-3,   // V
                   double Vr      = -70.4e-3,   // V
                   double El      = -70.4e-3,   // V
                   double gl      = 30e-9,      // S
                   double Cm      = 281e-12,    // F
                   double tau_w   = 144e-3,     // s
                   double slope   = 2e-3,       // V
                   double Vpeak   = 20.0e-3,    // V
                   double Vinit   = -60.0e-3,   // V
                   double Inoise  = 0.0,
                   double Iinject = 0.0);
    
    virtual ~aEIFNeuronBase()
    {
        // NOOP
    }

    virtual int adjust(double dt);

    virtual int advance(AdvanceInfo const &); 

    virtual int reset(double dt);

    inline bool refractoryQ()
    {
        return false;
    }

    virtual void clearSynapticInput(void);


    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;

    virtual double getManagedDelay() const { return 0; }
    virtual int nSpikeInputPorts() const { return 0; };
    virtual int nSpikeOutputPorts() const { return 1; };
    virtual int nAnalogInputPorts() const { return 0; };
    virtual int nAnalogOutputPorts() const { return 0; };

    virtual PortType outputPortType(port_t o) const
    {
        if( o==0) return spiking; else return undefined;
    };

    virtual PortType inputPortType(port_t i) const
    {
        return undefined;
    };

protected:
    double Vm;
    double w;
    double _dt;

    static ThreadSpecificRandomDistribution< NormalDistribution > white_noise;
/*
    PinkGenerator pink_noise;
*/
    int noiseType;
    double Isyn;

protected:
    double factI;
};



class aEIFNeuron : public aEIFNeuronBase, public CurrentInputTarget
{
    SIMOBJECT(aEIFNeuron, AdvancePhase::One)

public:
    aEIFNeuron(double a       = 4e-9,                   // S
               double b       = 0.0805e-9,              // A
               double Vt      = -50.6e-3,               // V
               double Vr      = -70.4e-3,               // V
               double El      = -70.4e-3,               // V
               double gl      = 30e-9,                  // S
               double Cm      = 281e-12,                // F
               double tau_w   = 144e-3,                 // s
               double slope   = 2e-3,                   // V
               double Vpeak   = 20.0e-3,                // V
               double Vinit   = -60.0e-3,               // V
               double Inoise  = 0.0,
               double Iinject = 0.0);

    virtual ~aEIFNeuron() {
        // NOOP
    }

    virtual void currentInput(double Isyn);
};




class CbaEIFNeuron : public aEIFNeuronBase, public ConductanceInputTarget
{
    SIMOBJECT(CbaEIFNeuron, AdvancePhase::One)

public:
    CbaEIFNeuron(double a       = 4e-9,         // S
                 double b       = 0.0805e-9,    // A
                 double Vt      = -50.6e-3,     // V
                 double Vr      = -70.4e-3,     // V
                 double El      = -70.4e-3,     // V
                 double gl      = 30e-9,        // S
                 double Cm      = 281e-12,      // F
                 double tau_w   = 144e-3,       // s
                 double slope   = 2e-3,         // V
                 double Vpeak   = 20.0e-3,      // V
                 double Vinit   = -60.0e-3,     // V
                 double Inoise  = 0.0,
                 double Iinject = 0.0);

    virtual ~CbaEIFNeuron() {
        // NOOP
    }

    virtual void currentInput(double Isyn);
    virtual void conductanceInput(double g, double Erev);

    virtual double getVm()
    {
        return Vm;
    }
};

#endif

---