+}
+
+// Compute uncertainty (stddev) of mu estimates, which is 1/sqrt(H_ii),
+// where H is the Hessian (see construct_hessian()).
+void compute_mu_uncertainty(const float *mu, int num_players)
+{
+ memset(mu_stddev, 0, sizeof(mu_stddev));
+
+ for (unsigned i = 0; i < all_matches.size(); ++i) {
+ const match &m = all_matches[i];
+
+ int p1 = m.player;
+ int p2 = m.other_player;