sim_epsp_stp.cpp File Reference

Example simulation which simulates one Poisson input with short-term plasticity onto a single postsynaptic neuron and records the membrane potential. More...

#include "auryn.h"

Include dependency graph for sim_epsp_stp.cpp:

Macros
#define	N   1

Functions
int	main (int ac, char *av[])

Detailed Description

Example simulation which simulates one Poisson input with short-term plasticity onto a single postsynaptic neuron and records the membrane potential.

As opposed to sim_epsp.cpp this examples uses a connection with short-term plasticity (Tsodyks Markram model) which is implemented in STPConnection.

Macro Definition Documentation

N

#define N   1

Function Documentation

main()

int main	(	int	ac,
		char *	av[]
	)

{

        int errcode = 0;
        char strbuf [255];
        string outputfile = "out_epsp_stp";
        string tmpstr;
        AurynWeight w = 1.0;

        // BEGIN Global definitions
        auryn_init( ac, av );
        sys->set_simulation_name(outputfile);
        // END Global definitions

        // Sets up a single presynaptic Poisson neuron which fires at 1Hz
        PoissonGroup * poisson = new PoissonGroup(N,1.);

        // Sets up a single postsynaptic integrate-and-fire neuron
        IFGroup * neuron = new IFGroup(1);


        // Initializes the STP connection
        STPConnection * con = new STPConnection(poisson,neuron,w,1.0,GLUT);

        // Sets STP parameters (depression dominated)
        double U = 0.2;
        double taud = 200e-3; // s
        double tauf = 50e-3;  // s

        // Passes the parameters to the connection instance
        con->set_tau_d(taud);
        con->set_tau_f(tauf);
        con->set_ujump(U);


        // Sets up recording

        // Records the input spike train
        tmpstr = outputfile;
        tmpstr += ".ras";
        SpikeMonitor * smon = new SpikeMonitor( poisson, tmpstr.c_str() );

        // Records the postsynaptic membrane potential
        tmpstr = outputfile;
        tmpstr += ".mem";
        VoltageMonitor * vmon = new VoltageMonitor( neuron, 0, tmpstr.c_str() );

        // Records the postsynaptic AMPA conductance
        tmpstr = outputfile;
        tmpstr += ".ampa";
        StateMonitor * amon = new StateMonitor( neuron, 0, "g_ampa", tmpstr.c_str() );

        // Records the postsynaptic NMDA conductance
        tmpstr = outputfile;
        tmpstr += ".nmda";
        StateMonitor * nmon = new StateMonitor( neuron, 0, "g_nmda", tmpstr.c_str() );

        // simulate for 5s
        logger->msg("Running ...",PROGRESS);
        sys->run(5);

        // Changes firing rate of PoissonGroup
        poisson->set_rate(50.0);

        // simulate for 0.5s
        sys->run(0.5);

        // Changes firing rate of PoissonGroup
        poisson->set_rate(1.0);

        // simulate for 4.5s
        sys->run(4.5);


        if (errcode)
                auryn_abort(errcode);

        logger->msg("Freeing ...",PROGRESS,true);
        auryn_free();

        return errcode;
}

Here is the call graph for this function: