def test_iba_vs_rayleigh_active_m2():
em_iba, em_ray = setup_func_rayleigh()
mu = setup_mu(1. / 64, bypass_exception=True)
ok_((abs(em_iba.ft_even_phase(2, mu, npol=3) / em_iba.ks - em_ray.ft_even_phase(2, mu, npol=3) / em_ray.ks) < tolerance_pc).all())
评论列表
文章目录