DoubleExponentialSpikeResponse.h

Go to the documentation of this file.

#ifndef DOUBLEEXPONENTIALSPIKERESPONSE_H_
#define DOUBLEEXPONENTIALSPIKERESPONSE_H_

#include "FiniteSpikeResponse.h"

class DoubleExponentialSpikeResponse : public FiniteSpikeResponse
{

        SIMOBJECT( DoubleExponentialSpikeResponse, AdvancePhase::SpikeDriven )

public:
    DoubleExponentialSpikeResponse(
        const float tau1=1e-3,
        const float tau2=3e-3
    ): FiniteSpikeResponse(), tau1(tau1), tau2(tau2)
    { /* NOOP */ };

    virtual ~DoubleExponentialSpikeResponse()
    { /* NOOP */ };

    float tau1;

    float tau2;

    // Update internal variables
    virtual int adjust( double dt );

    // Reset to initial condition.
    virtual int reset( double dt );

    virtual int advance(AdvanceInfo const &);

    virtual int spikeHit( spikeport_t port, SpikeEvent const& spike );


    virtual int psrLength(double dt) const;

private:

        inline double C1( double const& delta ) {
                return exp( -delta/tau1 );
        };

        inline double C2( double const& dt ) {
                return exp( - dt*(1.0/tau1+1.0/tau2) ) * ( exp(dt/tau2) - exp(dt/tau1) ) * pow( tau1 / tau2, (tau1+tau2)/(tau1-tau2) ) /
                       ( pow( tau1/tau2, tau1/( tau1-tau2 ) ) - pow( tau1/tau2 , tau2/(tau1-tau2) )  );
        };

        inline double C3( double const& delta ) {
                return exp( - delta/tau2 );
        };

    double c1,c2,c3;

    double x;
};


#endif /*DOUBLEEXPONENTIALSPIKERESPONSE_H_*/

---