figs/fig_4.py [107:118]:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
        x = [d_use for d_use in range(21) if d_use >= d]
        y1 = np.array([res['rembo'][(100, d, d_use)] for d_use in x])
        y2 = np.array([res['hesbo'][(100, d, d_use)] for d_use in x])
        y3 = np.array([res['unitsphere'][(100, d, d_use)] for d_use in x])
        y1err = 2 * np.sqrt(y1 * (1 - y1) / nsamp)
        y2err = 2 * np.sqrt(y2 * (1 - y2) / nsamp)
        y3err = 2 * np.sqrt(y3 * (1 - y3) / nsamp)
        ax.errorbar(x, y1, yerr=y1err, color=plt.cm.tab10(0), marker='')
        ax.errorbar(x, y2, yerr=y2err, color=plt.cm.tab10(1), marker='')
        ax.errorbar(x, y3, yerr=y3err, color=plt.cm.tab10(2), marker='')
        ax.set_title(r'$d={d}$'.format(d=d))
        if i == 0:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -



figs/fig_S5.py [19:30]:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
        x = [d_use for d_use in range(21) if d_use >= d]
        y1 = np.array([res['rembo'][(100, d, d_use)] for d_use in x])
        y2 = np.array([res['hesbo'][(100, d, d_use)] for d_use in x])
        y3 = np.array([res['unitsphere'][(100, d, d_use)] for d_use in x])
        y1err = 2 * np.sqrt(y1 * (1 - y1) / nsamp)
        y2err = 2 * np.sqrt(y2 * (1 - y2) / nsamp)
        y3err = 2 * np.sqrt(y3 * (1 - y3) / nsamp)
        ax.errorbar(x, y1, yerr=y1err, color=plt.cm.tab10(0), marker='')
        ax.errorbar(x, y2, yerr=y2err, color=plt.cm.tab10(1), marker='')
        ax.errorbar(x, y3, yerr=y3err, color=plt.cm.tab10(2), marker='')
        ax.set_title(r'$d={d}$'.format(d=d))
        if i == 0:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -