int main(int argc, char* argv[]) {
MPI_par mp;
setup_mpi(mp, argc, argv);
if (argc != 2) {
cerr << "\n\tUsage: ./abc_mpi abc_config_file.json\n\n";
return 100;
}
const gsl_rng* RNG = gsl_rng_alloc (gsl_rng_taus2);
gsl_rng_set(RNG, time (NULL) * getpid()); // seed the rng using sys time and the process id
AbcSmc* abc = new AbcSmc(mp);
abc->set_simulator(simulator);
abc->parse_config(string(argv[1]));
abc->run(RNG);
MPI_Finalize();
return 0;
}
int main(int argc, char* argv[]) {
MPI_par mp;
setup_mpi(mp, argc, argv);
if (argc != 2) {
cerr << "\n\tUsage: ./abc_mpi abc_config_file.json\n\n";
return 100;
}
// const gsl_rng* RNG = gsl_rng_alloc (gsl_rng_taus2);
gsl_rng_set(RNG, time (NULL) * getpid()); // seed the rng using sys time and the process id
AbcSmc* abc = new AbcSmc(mp);
abc->set_simulator(simulator);
abc->parse_config(string(argv[1]));
time(&GLOBAL_START_TIME);
abc->run(RNG);
//caseFile << "vaccine retro catchup target year inf1 inf2 inf3 inf4 sym1 sym2 sym3 sym4" << endl;
MPI_Finalize();
return 0;
}