/**
* \copydoc Simulation::simulate()
*/
void MetropolisHastingsSimulation::simulate() {
GetRNGstate();
lScores[0].fill(0);
lDerivative[0].fill(0);
LOGS(Priority::DEBUG)<<"Simulate with theta: "<<rParameters().transpose();
for (int m = 0; m < lpData->observationCount() - 1; ++m) {
printRNGTrace
simulatePeriod(m);
}
PutRNGstate();
lMeanScoresMinusTargets = lScores[0];
fireResult(lResult);
}
/**
* \copydoc Simulation::simulate(bool)
*/
void StatisticsSimulation::simulate(bool withSeeds) {
simulatePeriods(withSeeds);
LOGS(Priority::DEBUG)<<"simulation finished, distribute results";
fireResultModificators(lResult);
fireResult(lResult);
}
