def run_simulation1()

in figs/fig_4.py [0:0]

• 21 lines of code
• 4 McCabe index (conditional complexity)

def run_simulation1():
    t1 = time.time()
    nsamp = 1000
    res = {'rembo': {}, 'hesbo': {}, 'unitsphere': {}}
    D = 100
    for d in [2, 6, 10]:
        for d_use in range(1, 21):
            if d_use < d:
                continue
            res['rembo'][(D, d, d_use)] = p_A_contains_optimizer(
                d=d, D=D, d_use=d_use, gen_A_fn=gen_A_rembo, nsamp=nsamp
            )
            res['hesbo'][(D, d, d_use)] = p_A_contains_optimizer(
                d=d, D=D, d_use=d_use, gen_A_fn=gen_A_hesbo, nsamp=nsamp
            )
            res['unitsphere'][(D, d, d_use)] = p_A_contains_optimizer(
                d=d, D=D, d_use=d_use, gen_A_fn=gen_A_unitsphere, nsamp=nsamp
            )
    print(time.time() - t1)
    with open('data/fig4_sim_output.pckl', 'wb') as fout:
        pickle.dump(res, fout)