int maass()
{
/*
зависимость выходной частоты от тока в первом эксперименте
зависимость порога от характера сигнала
зависимость порога от реализации
зависимость частоты выходного сигнала от порога
использовать сигнал кости
*/
srand(time(NULL));
double timeMax = 50.0;
double dt = 0.01;
double dumpPeriod = 10.0;
Node::dt = dt;
std::vector<Node *> generators;
std::vector<Synapse *> synapses1;
std::vector<Synapse *> synapses2;
std::cout<<"Creating neurons..."<<std::endl;
Node node1 = Node("leaky iaf");
node1.I_stim = genran::uniform(14.5, 15.5);
node1.addPsWaveType(PsWave(3.0), "current", 1.0); // exc input
node1.addPsWaveType(PsWave(6.0), "current", -1.0); // inh input
node1.initialiseSpikeBuffer(dt);
Node node2 = Node("leaky iaf");
node2.I_stim = node1.I_stim;
node2.addPsWaveType(PsWave(3.0), "current", 1.0); // exc input
node2.addPsWaveType(PsWave(6.0), "current", -1.0); // inh input
node2.initialiseSpikeBuffer(dt);
std::cout<<"Creating generators..."<<std::endl;
for (int i=0; i < 90; i++){
generators.push_back(new Node("poisson generator"));
Synapse* synapse1 = new Synapse("stdp and tm", &node1, "hard boundaries");
synapse1->setPreset(0); // excitatory constants
synapse1->control(0); // disable stdp and shuffle tm
synapse1->reset();
topology::connectN2N(generators.back(), &node1, synapse1);
synapses1.push_back(synapse1);
Synapse* synapse2 = synapse1->duplicate(&node2); //copy everything
synapse2->control(1); // enable stdp with gamma, don't touch tm
synapse2->reset();
topology::connectN2N(generators.back(), &node2, synapse2);
synapses2.push_back(synapse2);
std::cout<<"\t\t generators created: "<< i << "%\r";
}
for (int i=0; i < 10; i++){
generators.push_back(new Node("poisson generator"));
Synapse* synapse1 = new Synapse("stdp and tm", &node1, "hard boundaries");
synapse1->setPreset(2); // inhibitory constants
synapse1->control(0); // disable stdp and shuffle tm
synapse1->reset();
synapse1->waveType = 1; // inhibitory
topology::connectN2N(generators.back(), &node1, synapse1);
// synapses1.push_back(synapse1);
Synapse* synapse2 = synapse1->duplicate(&node2); //copy everything
topology::connectN2N(generators.back(), &node2, synapse2);
// synapses2.push_back(synapse2);
std::cout<<"\t\t generators created: "<< i+90 << "%\r";
}
std::cout<<std::endl;
std::cout<<"Start output..."<<std::endl;
Output output;
output.openSpikeFile();
output.openSynapticWeightFile("data/target.txt");
output.synapses = synapses1;
output.printSynapticWeights(0.0);
output.closeSynapticWeightFile();
output.openSynapticWeightFile("data/start.txt");
output.synapses = synapses2;
output.printSynapticWeights(0.0);
output.closeSynapticWeightFile();
output.printEssentialData(timeMax);
std::string s;
std::stringbuf dumpTimesBuffer;
std::ostream dumpTimes (&dumpTimesBuffer);
std::cout<<"Starting simulation..."<<std::endl;
for (double t = 0.0; t <= timeMax+dt; t += dt){
for (unsigned i=0; i<generators.size(); i++)
if (generators.at(i)->step(t))
output.push(t, i); //generator numbers range from 0 to 99
if (node2.step(t))
output.push(t, 199);
if (node1.step(t) and not node2.isSpiking(t)){
node2.forceSpike(t);
output.push(t, 151);
}
if (isPeriodNow(t, dumpPeriod, dt) and t>0.0){
s = std::string("data/dump_") + std::to_string(t)+ ".txt";
dumpTimes << s <<std::endl;
output.openSynapticWeightFile(s);
output.print(t);
output.closeSynapticWeightFile();
}
std::cout<<"Finished "<< t / timeMax <<" of the simulation.\t\t\r";
}
std::cout<<std::endl;
output.close();
std::cout<<"Finish output..."<<std::endl;
output.openSynapticWeightFile("data/finish.txt");
output.printSynapticWeights(timeMax);
output.closeSynapticWeightFile();
std::ofstream dumpTimesFile;
dumpTimesFile.open("data/dump_times.txt", std::ofstream::out);
dumpTimesFile << dumpTimesBuffer.str();
dumpTimesFile.close();
return 0;
}
