assert_sigfigs_equal(np.array([1/2]*2), np.mean(samples, axis=0), 2)