density_3d_numerical#

jeanspy.model.SersicModel.density_3d_numerical

Class and construction: SersicModel.

SersicModel.density_3d_numerical(r_pc, *, epsrel=1e-06, epsabs=0.0, limit=200)[source]#

Deproject the Sérsic surface density by numerical Abel inversion.

Notes

Inputs and units. Nonnegative r_pc in pc; epsrel/epsabs are adaptive-integral tolerances; limit is the subdivision limit.

Returns and shape. pc^-3 with radius shape; infinity at r=0 for n>=1, a finite analytic center for n<1, zero at +infinity.

Validity. Positive finite re_pc and supported interpolation-table n. At r>0 integrate over theta in [0,pi/2].

Errors. Negative/NaN radii raise ValueError; SciPy convergence warnings can propagate. The returned value is not accompanied by a certified deprojection error.

Parameters: