--- trunk/gadget/simulation.cc 2014/02/10 17:09:07 1 +++ trunk/gadget/simulation.cc 2015/07/23 19:00:38 11 @@ -101,8 +101,11 @@ for (j = 0; j < Area->numAreas(); j++) this->updatePopulationOneArea(j); - for (j = 0; j < likevec.Size(); j++) + for (j = 0; j < likevec.Size(); j++){ likevec[j]->addLikelihood(TimeInfo); +// cout << "likelihood:" <getUnweightedLikelihood() << endl; + } +// cout << "###################################" << endl; if (print) for (j = 0; j < printvec.Size(); j++)