test_profiles.py 文件源码

def test_hernquist_mass():
    p = hernquist_mass_profile()
    a = 350.*apu.kpc
    M_0 = 1.0e15*apu.Msun
    p.set_param_values(a=a, M_0=M_0)
    assert_allclose(p(r_ap).to("solMass").value, (M_0*r_ap**2/(r_ap+a)**2).value)