enclosed_mass_betainc#

jeanspy.model_jax.ZhaoModel.enclosed_mass_betainc

Class and construction: ZhaoModel.

ZhaoModel.enclosed_mass_betainc(r_pc, *, params)[source]#

Enclosed mass from the Zhao incomplete-beta closed form.

For beta <= 3 (outside the beta domain), or a saturated beta argument, use the regularized numerical integral. Shape autodiff is unsupported on this explicit analytic path; use auto/numeric for inference.

The NFW-limit branch (alpha,beta,gamma)=(1,3,1) is handled analytically because the raw beta/betainc expression becomes indeterminate there even though the physical enclosed mass remains finite.

Parameters:
Return type:

jax.Array