#include #include #include int main() { int i; /* individual index */ int p; /* population index */ int o; /* observation index */ char buf[50]; FILE *f; Random *r = ran_random(); /* We repeat the process in order to construct 50 populations. */ for (p = 0; p < 50; p++) { sprintf(buf, "synthetic_populations_%2d.txt", p); f = fopen(buf, "w"); if (!f) { ran_random_free(r); return -1; } for (i = 0; i < 696; i++) { /* Group size */ fprintf(f, "%d; ", 1 + ran_intrandom(r, 3)); /* Week day */ fprintf(f, "%d; ", 1 + ran_intrandom(r, 7)); /* Hour. I assume we cannot leave after 7pm */ fprintf(f, "%d; ", 5 + ran_intrandom(r, 14)); /* Observations */ for (o = 0; o < 16; o++) { fprintf(f, "%g; ", ran_gumbel(r, -M_EULER, 1.0)); fprintf(f, "%g; ", ran_gumbel(r, -M_EULER, 1.0)); } fprintf(f, "\n"); ran_random_resetnextsubstream(r); } } ran_random_free(r); return 0; }