def sample_theta(): """Sample N=20 theta evolution""" ctrl, p, cdict, pdict = dec_regimes() sim = SimWrap(ctrl, p, cdict, pdict) sim.ctrl.keep_intermediate_values = True sim.show_CFO = False sim.set_all_nodisp() sim.show_TO = False sim.simulate() ax, _ = sim.post_sim_plots(save_TO='', save_CFO='')
def init_simwrap(): sim = SimWrap(ctrl, p, cdict, pdict); sim.set_all_nodisp() sim.TO_show_clusters = False sim.ctrl.saveall = False return sim