def test_u_dt_vs_point_source(self): """test that in the stationary reference frame, for large enough distances, the energy density of the DT tends to the one of a point like source approximating it""" # point source with the same luminosity as the DT ps_dt = PointSourceBehindJet(dt_test.xi_dt * dt_test.L_disk, dt_test.epsilon_dt) # r >> R_dt r = np.logspace(21, 24, 10) * u.cm assert u.allclose( dt_test.u(r), ps_dt.u(r), atol=0 * u.Unit("erg cm-3"), rtol=1e-2 )
def test_u_blr_vs_point_source(self): """test that for large enough distances the energy density of the BLR tends to the one of a point like source approximating it""" # point source with the same luminosity as the BLR ps_blr = PointSourceBehindJet(BLR.xi_line * BLR.L_disk, BLR.epsilon_line) # r >> R_line r = np.logspace(19, 23, 10) * u.cm assert u.allclose(BLR.u(r), ps_blr.u(r), atol=0 * u.Unit("erg cm-3"), rtol=1e-2)
def test_u_dt_vs_point_source_comoving(self): """test that in the reference frame comoving with the Blob, for large enough distances, the energy density of the DT tends to the one of a point like source approximating it""" # point source with the same luminosity as the DT ps_dt = PointSourceBehindJet(DT.xi_dt * DT.L_disk, DT.epsilon_dt) # r >> R_line r = np.logspace(21, 24, 10) * u.cm assert u.allclose( DT.u(r, PWL_BLOB), ps_dt.u(r, PWL_BLOB), atol=0 * u.Unit("erg cm-3"), rtol=1e-1, )