def interd_compare_r12():
graphs.change_fontsize(graphs.FONTSIZE + 2)
ctrl, p, cdict, pdict = dec_dpll()
sim = SimWrap(ctrl, p, cdict, pdict)
sim.TO_show_clusters = False
sim.ctrl.saveall = False
sim.simulate()
sim.post_sim_plots(save_TO='dpll_theta_evol')
ctrl, p, cdict, pdict = dec_r12()
sim = SimWrap(ctrl, p, cdict, pdict)
sim.TO_show_clusters = False
sim.ctrl.saveall = False
sim.simulate()
sim.post_sim_plots(save_TO='r12_theta_evol')
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='')
