GenericFroemkeDanStdpSynapse.h

Go to the documentation of this file.

#ifndef GENERICFROEMKEDANSTDPSYNAPSE_H_
#define GENERICFROEMKEDANSTDPSYNAPSE_H_

#include "GenericEachPairStdpSynapse.h"
#include "GenericOnePairStdpSynapse.h" 
#include "GenericStaticSpikingSynapse.h"
#include "GenericDynamicSpikingSynapse.h"
#include "GenericCurrentBasedSpikingSynapse.h"
#include "ExponentialDecaySpikeResponse.h"

#include <cmath>

#include <boost/format.hpp>
using boost::format;



template<class BaseSyn>
class GenericFroemkeDanStdpSynapse: public BaseSyn
{
public:
    bool useFroemkeDanSTDP;

    float tauspost;

    float tauspre;

    float taupos;

    float tauneg;

    float STDPgap;

    float Wex;

    float Aneg;

    float Apos;

    float mupos;

    float muneg;

    virtual double stdpLearning(const double & delta, const double  & t_post, const double & t_pre, const double & t_prev_post, const double & t_prev_pre );

    virtual double maxRelevantSpikeTimeDiff()
    {
        return 6*(fabs(tauneg)+fabs(taupos));
    }

};

#ifdef _MSC_VER
# pragma warning(disable:4244) // warning C4244: '+=' : conversion from 'double' to 'float', possible loss of data
#endif

template<class BaseSyn>
double GenericFroemkeDanStdpSynapse<BaseSyn>::stdpLearning(const double & delta, const double  & t_post, const double & t_pre, const double & t_prev_post, const double & t_prev_pre )
{
    double dw = 0;
    // cerr << format("inside learning delta=%10f ") % delta << " t_post=" << t_post << "t_pre=" << t_pre << endl;
    if (delta < -STDPgap) {
        // Depression
        // dw = Aneg * exp(delta/(fabs(Tmin)/3));
        dw = pow(BaseSyn::W, muneg) * Aneg * exp(delta/tauneg);
    } else if (delta > STDPgap) {
        // Potentiation
        // dw = Apos * exp(-delta/(fabs(Tmax)/3));
        dw = pow(fabs(Wex-BaseSyn::W), mupos) * Apos * exp(-delta/taupos);
    } else {
        return 0.0;
    }

        double epost = 1.0;
        double epre  = 1.0;
    if ( useFroemkeDanSTDP ) {
        if( t_prev_post > -1 ) epost = 1.0 - exp( -(t_post-t_prev_post) / tauspost );        
        if( t_prev_pre > -1 ) epre = 1.0 - exp( -(t_pre-t_prev_pre) / tauspre );
        BaseSyn::W += epost * epre * dw;
        // cerr << format("inside learning dw = %10e epost=%10e epre=%10e -> %10e") % dw % epost % epre % BaseSyn::W << endl;
    } else {            
        BaseSyn::W += dw;
        // cerr << format("inside learning dw = %10e -> %10e") % dw % BaseSyn::W << endl;
    }
    
    // check the sign
    if ((Wex < 0 && BaseSyn::W >0) || (Wex >0 && BaseSyn::W<0)) BaseSyn::W = 0;

    // check for greater Wmax
    if (fabs(BaseSyn::W) > fabs(Wex)) BaseSyn::W = Wex;
    
    return dw * epost * epre;
}

#ifdef _MSC_VER
# pragma warning(default:4244)
#endif

#endif /*GENERICFROEMKEDANSTDPSYNAPSE_H_*/

---