def run_simulation()

in figs/fig_S6.py [0:0]

• 15 lines of code
• 5 McCabe index (conditional complexity)

def run_simulation():
    t1 = time.time()
    nsamp = 1000
    res = {'unitsphere': {}}
    for D in [50, 100, 200]:
        for d in range(2, 19, 2):
            for d_use in range(2, 21, 2):
                if d_use < d:
                    continue
                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
                )
    with open('data/figS6_sim_output.pckl', 'wb') as fout:
        pickle.dump(res, fout)
    print(time.time() - t1)