LifNeuron.h

Go to the documentation of this file.

#ifndef _LIFNEURON_H_
#define _LIFNEURON_H_

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



class LifNeuronBase : public SimObject, public SingleOutputSpikeSender 
{
public:
    LifNeuronBase( float Rm       = 1e6,
                   float Cm       = 3e-8,
                   float Vresting =-0.06,
                   float Vthresh  =-0.045,
                   float Vreset   =-0.06,
                   float Vinit    =-0.06,
                   float Trefract = 3e-3,
                   float Inoise   = 0.0,
                   float Iinject  = 0.0);

    virtual ~LifNeuronBase()
    {
        /* NOOP */
    }

    virtual int reset( double dt );

    inline bool refractoryQ() { return (nStepsInRefr > 0); };

    virtual void clearSynapticInput(void)=0;

    float Rm;

    float Cm;

    float Vresting;

    float Vthresh;

    float Vreset;

    float Vinit;

    float Trefract;

    float Inoise;

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

    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;
    };
    
    
    virtual double currentNoiseInput()
    {
        return 0;
    }   

protected:
    double Vm;

    int nStepsInRefr;

    double C1,C2,I0;

        double _dt;
        
    static ThreadSpecificRandomDistribution< NormalDistribution > noise_gen;

};


class LifNeuron : public LifNeuronBase, public CurrentInputTarget
{
    SIMOBJECT( LifNeuron, AdvancePhase::One )

public:
    LifNeuron( float Rm       = 1e6,
               float Cm       = 3e-8,
               float Vresting =-0.06,
               float Vthresh  =-0.045,
               float Vreset   =-0.06,
               float Vinit    =-0.06,
               float Trefract = 3e-3,
               float Inoise   = 0.0,
               float Iinject  = 0.0);

    virtual ~LifNeuron()
    {
        /* NOOP */
    }

    virtual int adjust( double dt );

    virtual int advance(AdvanceInfo const &);

    virtual void clearSynapticInput(void);

    virtual void currentInput( double Isyn );

private:
    double Isyn;

};

inline void LifNeuron::clearSynapticInput(void)
{
    Isyn = 0;
}

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


class CbLifNeuron : public LifNeuronBase, public ConductanceInputTarget
{

    SIMOBJECT( CbLifNeuron, AdvancePhase::One )

public:
    CbLifNeuron( float Rm       = 1e6,
                 float Cm       = 3e-8,
                 float Vresting =-0.06,
                 float Vthresh  =-0.045,
                 float Vreset   =-0.06,
                 float Vinit    =-0.06,
                 float Trefract = 3e-3,
                 float Inoise   = 0.0,
                 float Iinject  = 0.0);

    virtual ~CbLifNeuron()
    {
        /* NOOP */
    }

    virtual int adjust( double dt );

    virtual int advance(AdvanceInfo const &);

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

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

private:
    double Isyn;
    double Gsyn;

};

inline void CbLifNeuron::clearSynapticInput(void)
{
    Isyn = Gsyn = 0;
}

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

inline void CbLifNeuron::conductanceInput( double g, double Erev )
{
    Gsyn += g;
    Isyn += (g*Erev);
}

#endif

---