test_profiles.py 文件源码

def test_beta_model():
    p = beta_model_profile()
    r_c = 100.*apu.kpc
    rho_c = 1.0e-25*apu.g/apu.cm**3
    beta = 1.0
    p.set_param_values(r_c=r_c, rho_c=rho_c, beta=beta)
    assert_allclose(p(r_ap).value, (rho_c*(1.+(r_ap/r_c)**2)**(-1.5*beta)).value)