DoubleExponentialSpikeResponse.cpp

Go to the documentation of this file.

#include "DoubleExponentialSpikeResponse.h"

int DoubleExponentialSpikeResponse::adjust( double dt )
{
    if ( tau1 <= 0 )
        throw( PCSIM::Exception( "ExponentialDecaySpikeResponse::adjust", "Time constant tau1 not a positive vlaue!" ) );

    if ( tau2 <= 0 )
        throw( PCSIM::Exception( "ExponentialDecaySpikeResponse::adjust", "Time constant tau2 not a positive vlaue!" ) );

    if( tau1 == tau2 )
        throw( PCSIM::Exception( "ExponentialDecaySpikeResponse::adjust", "Time constants tau1 and tau2 must differ. Use AlphaSpikeResponse instead!" ) );

    c1 = C1( dt  );
    c2 = C2( dt  );
    c3 = C3( dt  );

    return 0;
}

int DoubleExponentialSpikeResponse::reset( double dt )
{
    FiniteSpikeResponse::reset( dt );
    adjust( dt );
    psr = x = 0.0;
    return 0;
}

int DoubleExponentialSpikeResponse::advance(AdvanceInfo const &)
{
    // Do the exponantial decay
        psr = c2 * x + c3 * psr;
        x *= c1;
    return advanceReturn();
}

int DoubleExponentialSpikeResponse::spikeHit( spikeport_t port, SpikeEvent const& spike )
{
        double _c1_ = C1(spike.delta);
        double _c2_ = C2(spike.delta);
        double _c3_ = C3(spike.delta);
        double psr_tgt = _c2_ * x + _c3_ * psr;
        x   += spike.weight / _c1_;
    psr = (( psr_tgt - _c2_ * x ) / _c3_);
     return spikeHitReturn( AdvanceInfo( spike.dt ) );
}

int DoubleExponentialSpikeResponse::psrLength(double dt) const
{
    return (int)( (double)PSR_MULTIPLE_TAU * 2 * (tau1+tau2) / dt + 0.5 );
}

---