erf(0.5 * math.sqrt(2) * (e - m_xi) / std_xi)) * m_la #=========================================================================== # Lab #=========================================================================== slab = SPIRRIDLAB(s = s, save_output = True, show_output = True, exact_arr = mu_q_ex(e_arr, m_xi, std_xi, m_la)) #=========================================================================== # Compare efficiency of sampling types #=========================================================================== powers = np.linspace(1, math.log(1000, 10), 80) n_int_range = np.array(np.power(10, powers), dtype = int) # slab.sampling_efficiency(n_int_range = n_int_range) #=========================================================================== # Compare the structure of sampling #=========================================================================== # slab.sampling_structure(ylim = 18.0, xlim = 1.2,) #=========================================================================== # Compare the code efficiency #==============================8============================================= s.set(e_arr = np.linspace(0, 2.0, 40), n_int = 400) #slab.codegen_language_efficiency(extra_compiler_args = False) slab.codegen_language_efficiency(extra_compiler_args = True)