IzhiNeuron.h

Go to the documentation of this file.

#ifndef _IZHINEURON_H_
#define _IZHINEURON_H_

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



class IzhiNeuronBase : public SimObject, public SingleOutputSpikeSender
{
public:
    IzhiNeuronBase(double a       = 0.02,
                   double b       = 0.2,
                   double c       = -65.0,
                   double d       = 8.0,
                   double Vpeak   = 30.0,
                   double Vinit   = -0.06,
                   double Inoise  = 0.0,
                   double Iinject = 0.0,
                   double factI   = 1.0);

    IzhiNeuronBase(string type,
                   double Vinit   = -0.06,
                   double Inoise  = 0.0,
                   double Iinject = 0.0,
                   double factI   = 1.0);

    virtual ~IzhiNeuronBase()
    {
        // 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 Vpeak;
    double Vinit;
    double Inoise;
    double Iinject;


    double a;
    double b;
    double c;
    double d;

    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 u;
    double dt_msec;


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

protected:
    double factI;
};




class IzhiNeuron : public IzhiNeuronBase, public CurrentInputTarget
{
    SIMOBJECT(IzhiNeuron, AdvancePhase::One)

public:
    IzhiNeuron(double a       = 0.02,
               double b       = 0.2,
               double c       = -65.0,
               double d       = 8.0,
               double Vpeak   = 30.0,
               double Vinit   = -0.06,
               double Inoise  = 0.0,
               double Iinject = 0.0,
               double factI   = 1.0);


    IzhiNeuron(string type,
               double Vinit   =-0.06,
               double Inoise  = 0.0,
               double Iinject = 0.0,
               double factI   = 1.0);

    virtual ~IzhiNeuron() {
        // NOOP
    }

    virtual void currentInput(double Isyn);
};




class CbIzhiNeuron : public IzhiNeuronBase, public ConductanceInputTarget
{
    SIMOBJECT(CbIzhiNeuron, AdvancePhase::One)

public:
    CbIzhiNeuron(double a       = 0.02,
                 double b       = 0.2,
                 double c       =-65.0,
                 double d       = 8.0,
                 double Vpeak   = 30.0,
                 double Vinit   =-0.06,
                 double Inoise  = 0.0,
                 double Iinject = 0.0,
                 double factI   = 1.0);
/*
                 int NoiseType  = 0,
                 double ge0     = 0.012e-6, 
                 double sig_e   = 0.003e-6, 
                 double tau_e   = 2.7e-3, 
                 double gi0     = 0.057e-6, 
                 double sig_i   = 0.0066e-6, 
                 double tau_i   = 10.5e-3);
*/

    CbIzhiNeuron(string type,
                 double Vinit   =-0.06,
                 double Inoise  = 0.0,
                 double Iinject = 0.0,
                 double factI   = 1.0);
/*
                 int NoiseType  = 0,
                 double ge0     = 0.012e-6, 
                 double sig_e   = 0.003e-6, 
                 double tau_e   = 2.7e-3, 
                 double gi0     = 0.057e-6, 
                 double sig_i   = 0.0066e-6, 
                 double tau_i   = 10.5e-3); 
*/

    virtual ~CbIzhiNeuron() {
        // NOOP
    }

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

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




class ExIzhiNeuronBase : public IzhiNeuronBase
{
public:
    ExIzhiNeuronBase(double a       = 0.02,
                     double b       = 0.2,
                     double c       = -65.0,
                     double d       = 8.0,
                     double k       = 0.7,
                     double Vt      = -60.0,
                     double Vr      = -45.0,
                     double Cm      = 100.0,
                     double Vpeak   = 30.0,
                     double Vinit   = -0.06,
                     double Inoise  = 0.0,
                     double Iinject = 0.0,
                     double factI   = 1.0);

    ExIzhiNeuronBase(string type,
                     double Vinit   = -0.06,
                     double Inoise  = 0.0,
                     double Iinject = 0.0,
                     double factI   = 1.0);

    virtual ~ExIzhiNeuronBase(){
        // NOOP
    }

    virtual int advance(AdvanceInfo const &); 

protected:
    double k;
    double Vr;
    double Vt;
    double Cm;
};



class ExIzhiNeuron : public ExIzhiNeuronBase, public CurrentInputTarget
{
    SIMOBJECT(ExIzhiNeuron, AdvancePhase::One)

public:
    ExIzhiNeuron(double a       = 0.02,
                 double b       = 0.2,
                 double c       = -65.0,
                 double d       = 8.0,
                 double k       = 0.7,
                 double Vt      = -60.0,
                 double Vr      = -45.0,
                 double Cm      = 100.0,
                 double Vpeak   = 30.0,
                 double Vinit   = -0.06,
                 double Inoise  = 0.0,
                 double Iinject = 0.0,
                 double factI   = 1.0);

    ExIzhiNeuron(string type,
                 double Vinit   =-0.06,
                 double Inoise  = 0.0,
                 double Iinject = 0.0,
                 double factI   = 1.0);

    virtual ~ExIzhiNeuron() {
        // NOOP
    }

    virtual void currentInput(double Isyn);
};




class CbExIzhiNeuron : public ExIzhiNeuronBase, public ConductanceInputTarget
{
    SIMOBJECT(CbExIzhiNeuron, AdvancePhase::One)

public:
    CbExIzhiNeuron(double a       = 0.02,
                   double b       = 0.2,
                   double c       =-65.0,
                   double d       = 8.0,
                   double k       = 0.7,
                   double Vt      = -60.0,
                   double Vr      = -45.0,
                   double Cm      = 100.0,
                   double Vpeak   = 30.0,
                   double Vinit   =-0.06,
                   double Inoise  = 0.0,
                   double Iinject = 0.0,
                   double factI   = 1.0);

    CbExIzhiNeuron(string type,
                   double Vinit   =-0.06,
                   double Inoise  = 0.0,
                   double Iinject = 0.0,
                   double factI   = 1.0);

    virtual ~CbExIzhiNeuron() {
        // NOOP
    }

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

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

#endif

---