/*! Take a single sample of time. If at least one previous sample has been taken, the difference
between the current time and the previous time is stored as the duration of that "lap" through the code. */
void CX_LapTimer::takeSample(void) {
_timePoints.push_back(_clock->now());
_durationRecalculationRequired = true;
if ((_samplesBetweenLogging != 0) && (_timePoints.size() == _samplesBetweenLogging)) {
CX::Instances::Log.notice("CX_LapTimer") << "Stats for last " << _samplesBetweenLogging << " samples." << getStatString();
restart();
}
}
void NEATRunner::runLoop()
{
//setting up signal handlers
(void)signal(SIGINT,signalhandler);
(void)signal(SIGTERM,signalhandler);
(void)signal(SIGKILL,signalhandler);
//
time_t startt,tmpt;long totaltime=0;
int countruns=0;
vector<double> beststate;
generations++;
Evaluator * bak = NULL;
if(nodes==0&&localFE==false){
localFE=true;
}else if(!localFE)
cerr << "running cluster code.." << endl;
writeRunfile(false,basefile,infoline,pid);
pop->fe = icb->fe;
stringstream sgfc; sgfc << sgf.str() << "-" << countruns << ".xml";
if(speciationGraph)
sg = new SpecGraph((int)pop->getMembers()->size(),generations,sgfc.str());
stringstream sCurrentGenomeFilesuffixless; sCurrentGenomeFilesuffixless << sCurrentGenomeFile.str();
stringstream sCurrentGenomeFilec; sCurrentGenomeFilec << sCurrentGenomeFile.str() << "-" << countruns;
stringstream sCurrentGenGenomeFilec; sCurrentGenGenomeFilec << sCurrentGenomeFilec.str();
stringstream sCurrentXMLGenomeFilec; sCurrentXMLGenomeFilec << sCurrentXMLGenomeFile.str() << "-" << countruns << ".xml";
stringstream sCurrentGraphFilec; sCurrentGraphFilec << sCurrentGraphFile.str() << "-" << countruns;
if(currentgraphf==NULL){
currentgraphf = new ofstream(sCurrentGraphFile.str().c_str());
}
ofstream ofs(sCurrentGenomeFilec.str().c_str());
ofstream ofs2(sCurrentGenomeFilesuffixless.str().c_str());
ofstream ofscurg(newestgenome.c_str());
bool slaveStop = false;
int * bestspecid = new int;
bool teststop = false;
int gc=0; double ftest = 0;
FitnessEvaluator * fe = ev->getFitnessEvaluator();
while(!stop){
if(coevo != NULL && pop->getGeneration() == (coevo->getStartGeneration()+1)){
cout << "doing coevo!?" << endl;
bak = ev;
ev = coevo;
}
teststop = false;
startt = time(0);
if(localFE){
ev->evaluate(pop->getMembers(),pop->getMembers()->size());
}else{
if(generations>0&&(pop->getGeneration()+2)==generations&&runs==(countruns+1))
slaveStop = true;
comm->outputPopulation(pop,nodes,coevo,mc,slaveStop); //stream the population out to nodes for evaluation
ev->evaluate(pop->getMembers(),mc);//sweet..
comm->readFitness(pop,mc); //read the corresponding returned fitness values
}
//checking and updating for the overall best phenotype.
//do population/species sorting and stat updating
pop->updateSpeciesStats();
pop->sortmembers();
pop->sortspecies();
avgf = pop->calcAvgFitness();
//keeping a copy of generation champ:
gbest = pop->getCopyOfCurrentBest();
setChamp(best,gbest,bestspecid);
if(ftest < pop->getMembers()->at(0)->getFitness()){
ftest = pop->getMembers()->at(0)->getFitness();
}
else if(ftest != 0 && ftest > pop->getMembers()->at(0)->getFitness()){
cout << "old genome:\n " << best->getGenome();
cout << "new genome:\n " << pop->getMembers()->at(0)->getGenome();
cout << "fitness went DOWN ftest: " << ftest << " pop->getMembers()->at(0)->getFitness(): " << pop->getMembers()->at(0)->getFitness() << endl;
Phenotype * p = pop->getMembers()->at(0);
vector<double> input = p->getNet()->getInput();
Phenotypes * testp = new Phenotypes();
testp->push_back(p);
testp->push_back(gbest);
for(unsigned int i=0;i<10;i++){
ev->evaluate(testp
, 1);
}
}
best->getGenome()->setSeed(rands);
sCurrentGenGenomeFilec.str("");
sCurrentGenGenomeFilec << sCurrentGenomeFilec.str();
sCurrentGenGenomeFilec << "-" << pop->getGeneration();
ofs.open(sCurrentGenomeFilec.str().c_str());
ofs << best->getGenome();
ofs.close();
ofs2.open(sCurrentGenomeFilesuffixless.str().c_str());
ofs2 << best->getGenome();
ofs2.close();
ofscurg.open(newestgenome.c_str());
ofscurg << best->getGenome();
ofscurg.close();
ofstream ofs3(sCurrentGenGenomeFilec.str().c_str());
ofs3 << best->getGenome();
ofs3.close();
// if(pop->getGeneration()%2==0){
// cerr << icb->fe->show(best);
// }
icb->best = best;
writenetwork(best,sCurrentXMLGenomeFilec.str());
if(speciationGraph)
sg->update(pop);
//writing stats to file
*currentgraphf << getStatString(pop,avgf);
currentgraphf->flush();
//updating smoothed graph data..
updateSmoothData(smoothdata,pop,avgf,countruns+1);
if(coevo!=NULL)
coevo->update(pop); // update the coevolution data..
tmpt = (time(0)-startt);
totaltime += tmpt;
cerr << (pop->getGeneration()+1) << ":"
<< " curmaxid: "<<pop->getMembers()->at(0)->getID()
<< "(" << pop->getMembers()->at(0)->getSpecies()->getID() << ")"
<< " curmax: " << pop->getMembers()->at(0)->getFitness()
<< " bestid: "<< best->getID()
<< "(" << *bestspecid << ")"
<< " bestfitness: "<< best->getFitness()
<< " maxfitness: " << pop->getHighestFitness()
<< " avgfitness: " << pop->calcAvgFitness()
<< " curmin: " << pop->getMembers()->at(pop->getMembers()->size()-1)->getFitness()
<< " species: " << pop->getSpecies()->size()
<< " (" << pop->spectarget << ") "
<< " size: " << pop->getMembers()->size()
<< " time: " << tmpt
<< " time/size: " << (double)tmpt/(double)pop->getMembers()->size() << endl;
//run the code that adjusts fitness according to species age and size..
//select the lucky ones for reprocicration..
sel->select(pop,sel->getySpeciesForElitism());
//do the mating
teststop = fe->stop(best);
rep->reproduce(pop);
if(((generations>0&&(pop->getGeneration()+1) ==generations) || teststop)
&& runs==(countruns+1)){ // stopconditions
if(teststop && comm!=NULL)
comm->outputPopulation(pop,nodes,coevo,mc,true);
gc += pop->getGeneration();
setChamp(sbest,best,bestspecid);
stop = true;
if(speciationGraph){
sg->writetofile();
delete sg;
}
}else{
if(generations>0&& ((pop->getGeneration()+1)==generations || teststop )){
gc += pop->getGeneration();
//generation run is over lets go on to the next run..
countruns++;
//keep superchamp across runs..
Phenotype * tmp = best;
best = NULL;
setChamp(sbest,tmp,bestspecid);
//reset population...
if(pop->spawn)
pop->resetSpawn();
else
pop->resetGenesis();
if(bak!=NULL){
ev = bak;
bak = NULL;
}
if(speciationGraph){
sg->writetofile();
delete sg;
sgfc.str(""); sgfc << sgf.str() << "-" << countruns << ".xml";
sg = new SpecGraph((int)pop->getMembers()->size(),generations,sgfc.str());
}
sCurrentGenomeFilec.str(""); sCurrentGenomeFilec << sCurrentGenomeFile.str() << "-" << countruns;
sCurrentXMLGenomeFilec.str(""); sCurrentXMLGenomeFilec << sCurrentXMLGenomeFile.str() << "-" << countruns << ".xml";
sCurrentGraphFilec.str(""); sCurrentGraphFilec << sCurrentGraphFile.str() << "-" << countruns;
ftest = 0;
delete currentgraphf; currentgraphf = new ofstream(sCurrentGraphFilec.str().c_str());
}
}
}
ofstream ofsuper(sFinalGenomeFile.c_str());
ofsuper << sbest->getGenome();
ofsuper.close();
cout << "final best fitness" << sbest->getFitness() << endl;
if(countruns>1){
cout << "avg gc: " << (double)gc/(double)(countruns+1) << endl;
}
delete sbest;
ofstream finalgraphf(finalgraphfile.c_str());
for(int i=0;i<generations-1;i++)
finalgraphf << smoothdata[i][0]/(double)runs << " "
<< smoothdata[i][1]/(double)runs << " "
<< smoothdata[i][2]/(double)runs << endl;
for(int i=0;i<generations-1;i++)
delete[] smoothdata[i];
delete[] smoothdata;
//close graph files..
finalgraphf.close();
currentgraphf->close();
delete currentgraphf;
delete bestspecid;
writeRunfile(true,basefile,infoline,pid);
}
/*! This function stores the difference between the current time and the time captured with t1().
If enough samples have been collected, equal to the value of `logSamples` during setup(), a
summary statistics string will be automatically logged. */
void CX_SegmentProfiler::t2(void) {
_durations.push_back((_clock->now() - _t1).value());
if ((_samplesBetweenLogging != 0) && (_durations.size() == _samplesBetweenLogging)) {
CX::Instances::Log.notice("CX_SegmentProfiler") << "Stats for last " << _samplesBetweenLogging << " samples." << getStatString();
restart();
}
}
