AlphaFunctionSpikeResponse.cpp

Go to the documentation of this file.

#include "AlphaFunctionSpikeResponse.h"

#include <iostream>
using std::cerr;
using std::endl;

int AlphaFunctionSpikeResponse::adjust( double dt )
{
    if ( tau > 0 ) {
        c1 = C1( dt );
        c2 = C2( dt );
    } else {
       throw( PCSIM::Exception( "AlphaFunctionSpikeResponse::adjust", "Time constant not a positive vlaue!" ) );
    }
    return 0;
}

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

int AlphaFunctionSpikeResponse::advance(AdvanceInfo const &)
{
    psr  = c2 * x + c1 * psr;
        x   *= c1;
    return advanceReturn();
}

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

int AlphaFunctionSpikeResponse::psrLength(double dt) const
{
    return (int)( (double)PSR_MULTIPLE_TAU * 3 * tau / dt + 0.5 );
}

---