+#if DUMP_RAW
+ dump_raw(mu, num_players);
+#else
+ construct_hessian(mu, num_players);
+ compute_mu_uncertainty(mu, 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("aux_param score_stddev %f\n", global_sigma / sqrt(2.0f));
+ printf("aux_param rating_prior_stddev %f\n", prior_sigma);
+
+ float total_logl = compute_total_logl(mu, num_players);
+ printf("aux_param total_log_likelihood %f\n", total_logl);
+#endif