+#if DUMP_RAW
+ dump_raw(mu, num_players);
+#else
+ compute_mu_uncertainty(mu, num_players);
+ dump_scores(players, mu, mu_stddev, num_players);
+ //fprintf(stderr, "Optimal sigma: %f (two-player: %f)\n", sigma[0], sigma[0] * sqrt(2.0f));
+ printf("%f -2\n", global_sigma / sqrt(2.0f));
+ printf("%f -3\n", prior_sigma);
+
+ float total_logl = compute_total_logl(mu, num_players);
+ printf("%f -4\n", total_logl);
+
+// construct_hessian(mu, sigma, num_players);
+#endif