sim_poisson.cpp File Reference

Example simulation which simulates a group of Poisson neurons. More...

#include "auryn.h"

Include dependency graph for sim_poisson.cpp:

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

Detailed Description

Example simulation which simulates a group of Poisson neurons.

Function Documentation

main()

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

{

        string dir = "./";
        string file_prefix = "poisson";

        char strbuf [255];
        string msg;

        NeuronID size = 1000;
        unsigned int seed = 1;
        double kappa = 5.;
        double simtime = 1.;

        int errcode = 0;

    try {

        po::options_description desc("Allowed options");
        desc.add_options()
            ("help", "produce help message")
            ("simtime", po::value<double>(), "simulation time")
            ("kappa", po::value<double>(), "poisson group rate")
            ("dir", po::value<string>(), "output directory")
            ("size", po::value<int>(), "poisson group size")
            ("seed", po::value<int>(), "random seed")
        ;

        po::variables_map vm;        
        po::store(po::parse_command_line(ac, av, desc), vm);
        po::notify(vm);    

        if (vm.count("help")) {
            std::cout << desc << "\n";
            return 1;
        }

        if (vm.count("kappa")) {
            std::cout << "kappa set to " 
                 << vm["kappa"].as<double>() << ".\n";
                        kappa = vm["kappa"].as<double>();
        } 

        if (vm.count("dir")) {
            std::cout << "dir set to " 
                 << vm["dir"].as<string>() << ".\n";
                        dir = vm["dir"].as<string>();
        } 

        if (vm.count("simtime")) {
            std::cout << "simtime set to " 
                 << vm["simtime"].as<double>() << ".\n";
                        simtime = vm["simtime"].as<double>();
        } 

        if (vm.count("size")) {
            std::cout << "size set to " 
                 << vm["size"].as<int>() << ".\n";
                        size = vm["size"].as<int>();
        } 

        if (vm.count("seed")) {
            std::cout << "seed set to " 
                 << vm["seed"].as<int>() << ".\n";
                        seed = vm["seed"].as<int>();
        } 
    }
    catch(std::exception& e) {
        std::cerr << "error: " << e.what() << "\n";
        return 1;
    }
    catch(...) {
        std::cerr << "Exception of unknown type!\n";
    }

        auryn_init(ac, av);
        sys->set_master_seed(seed);

        PoissonGroup * poisson = new PoissonGroup(size,kappa);

        sprintf(strbuf, "%s/%s.%d.ras", dir.c_str(), file_prefix.c_str(), sys->mpi_rank() );
        SpikeMonitor * smon_e = new SpikeMonitor( poisson, strbuf, size);

        sprintf(strbuf, "%s/%s.%d.prate", dir.c_str(), file_prefix.c_str(), sys->mpi_rank() );
        PopulationRateMonitor * pmon_e = new PopulationRateMonitor( poisson, strbuf, 1.0 );

        RateChecker * chk = new RateChecker( poisson , -1 , 20.*kappa , 10);
        if (!sys->run(simtime,false)) 
                        errcode = 1;


        if (errcode)
                auryn_abort(errcode);
        logger->msg("Freeing ...",PROGRESS,true);
        auryn_free();
        return errcode;
}

Here is the call graph for this function: