--- trunk/gadget/simulation.cc 2015/06/10 13:01:42 10 +++ 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++)